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16. Abstract When the average scattering properties of a non-coherent extensive scene with homogeneous 
statistics are measured, It is essential that an appropriate radar-scene interaction relation be 
employed to properly interpret the results. The interaction relationship hopefully will be sufficient 
ly general to isolate various polarized scattering parameters and to predict the average return for an 
arbitrarily polarized antenna. The complete non-coherent radar equation is derived within this work i 
and a measurement-inversion technique is developed to isolate all the scattering paramerers for a 
scene satisfying reciprocity. VJhen the measurement technique Is specialized to scenes having a 
scattering characteristic similar to the sea, antenna polarization specifications for accurate retriev 
al of the scattering parameters are established through simulations. The extrapolation of these 
results directly or indirectly through the use of the non-coherent radar equation can yield antenna 
polarization specifications for other scenes. The results serve as a guideline for designing meaning- 
ful and accurate scatterometer experiments. 

It is specifically shown that for scenes satisfying reciprocity, one must admit three complex 
valued scattering coefficients in addition to the three well known real valued scattering coefficients 
The complex valued coefficients are associated with the cross-correlation scattering properties of the 
scene. As a result of the spatial integration to acquire an average return, the scattering coeffici- 
ents must satisfy Schwarts' inequality. An immediate consequence of strict inequality for the complex- 
ed valued coefficients is that radar returns from non-coherent scenes must be partially polarized. 

The characteristics of the complex valued coefficients are demonstrated from scattering theories 
applicable to the sea. 

The measurement-inversion technique is proposed with and without regard to the difference between 
antenna and surface polarizations. It is demonstrated that the distinction between polarizations in 
negligible for narrow beam antenna at all but small view angles. For small view angles one must, in 
general, employ the inversion based on surface polarizations if he wishes to compare his measurements 
with theory or that of other experimenters. Inversions based on antenna polarizations can; however, 
be performed at small incident angle if very small beamwidths are employed. In this case, nadir can 
only be probed In en asymptotic sense. This is also a necessary procedure when an anisotropic charac- 
ter is to be measured at small incident angles. (continued) 
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16. Abstract (continued) 


Within the simulations it is shown that there is little difficulty in recovering the dominant 
scattering coefficients with modest realizations of the polarization specifications. Retrieval of the 
weaker scattering parameters requires more careful observation of the polarization requirements. In 
the latter case, it is shown that more relaxed realizations of the polarization specifications can be 
tolerated for many of the measurements if the phase of the cross polarized leakage can be adjusted to an 
□ptiniun value. 

In general it is indicated that three real and three complex valued scattering coefficients can 
interact with the scatterometer antenna in an undesirable fashion when attempting recovery of any one 
coefficient. The measurement error arises either as a result of inadequate realizations of the speci- 
fied antenna polarization or as a result of the inherent mis-match between antenna and surface 
polarizations for small angles. 
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This report was prepared by the Remote Sensing Laboratory of the University 
of Kansas Space Technology Laboratories under Contract NAS1''10048, Under this 
contract the principal investigator is Dr, R. K. Moore and the project engineer is 
Dr, Ad K. Fung. 

This document covers a particular task In an on-going effort between NASA 
Langley Research Center and the University of Kansas to demonstrate the value of 
the microwave scatterometer as a remote sea wind sensor. Specifically the inter- 
action between an arbitrarily polarized scatterometer antenno and a non“coherent 
distributive target is derived and applied to develop a measuring technique to recover 
all the scattering parameters. The results are helpful for specifying antenna polar- 
ization properties for accurate retrieval of the parameters not only for the sea but 
also for other distributive scenes . 
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LO SUMMARY 


The non-cohereni- radar equation is derived within the framework of a general- 
ized reception theory. For scenes satisfying reciprocity, the resulting equation confirms 
a previously derived theory [6]; this result, however, was extended to occount for the 
difference between antenna and surface polarizations. The present theory permits one 
to interpret the radar return and its reception within the context of scattering and co- 
herence theories (see Section 5.2), Under the reciprocity assumption it is shown that 
in addition to the three commonly known real valued scattering coefficients there are 
three complex valued coefficients (without reciprocity there are four real and six 
complex valued coefficients). As a result of the new coefficients, the definition of 
a scattering coefficient had to be extended. Specifically a descriptive definition was 
suggested, viz,. 
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where, 

9 

-^ij^jt = scattered field component 


= linear polarized scattering operator 
= incident field component 

R = range to the illuminated area 

A = incremental area of illumination 

0 = incident angle 

o 

The subscripts denote the polarization states of the incident and scattered fields, either 
vertical v or horizontal h. The above definition encompasses the old as well as the 
new scattering coefficients. In the new notation and denote 

the polarized coefficients, < is the cross polarized coefficient, and 

> , and are the new complexed valued coefficients. 

Other scattering coefficients participate in the scattering process as implied by Equation (I- 
i)‘ above; however, when reciprocity is satisfied the above set Is sufficient (See 
Equations (4-26) and (4-29).) 

The complex valued coefficients account for the relative phase induced between 
the vertically and horizontally polarized components by the scattering surface. The 
phase characteristic of these scattering coefficients interact with the relative phase pro- 


perttes of fhe transmjssion and recepHon antennas to contribute on observed power com** 
posed of real and complex valued scattering coefficients. This interaction occurs 
within the coherent radar equation also; however the interaction must be Intel- 
preted differently for the non-coherent case„ As a result of the spatial integration 
to acquire an average return, the complex valued coefficient must, in genera! satisfy 
Schwartz* inequality 

For the coherent case, equality is always assured q However, for the non~coherent 
case strict inequality can occur. As a result of the strict inequality, one can attribute 
a partially polarized character to non-coherent radar returns (See Section 5.4), Also 
as a result of the inequality, techniques for measuring the scattering matrix for coherent 
targets cannot be employed for non-coherent targets. To illustrate the character of 
these scattering coefficients, several scattering theories applicable to sea returns were 
examined (See Section 5,3) . 

On the basis of the above theory a measurement and inversion technique was 
developed to measure all six coefficients (nine parameters when the real and imaginary 
parts are considered). The technique is based on intensity measurements by narrow 
beam radar scatterometers (See Section 6.4), Inversions are proposed with and without 
regard to the distinction between antenna and surface polarizations (See Section 6,4 
and Section 7.2, respectively). It is demonstrated that the distinction between polar- 
izations is negligible for narrow beam antennas at all but small view angles (See Sections 
5,5 and 7,3), For small view angles, inversions based on surface polarizations are more 
accurate, in general, if the measurements are to b>j compared with theory or with other 
experimenters. For example, a 50% error occurs at nadir in inverting for > 

(defined with respect to the surface polarizations) when ine inversion technique is based 
on antenna polarizations. Comparison of the inversions with and without regard to the 
distinction are shown in Figures (7.3) through (7.8), Inversions based on antenna polar- 
izations can, however, be performed at small view angles if very small beamwidths 
are employed. In this case nadir can only be probed in an asymptotic sense. The degree 
to which one can approach nadir and yet meet the constraint that the antenna polar- 
izations across the main beam approximately match those of the surface is dependent on 
beamwidth. Figure (7.2) parametrically shows the beamwidth requirement as a function of 
view angle to minimize unwanted orthogonally polarized content in the measurement. 
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This laHer technique is preferred in as much as the measurements may be restricted to a 
partial set of coefficients whereas when inversions are performed with respect to the sur- 
face polarizations the entire set of measurements must be performed. It is also advantag- 
eous to use inversions based on antenna polarizations and small beamwidth antennas when 
an anisotropic characteristic is to be measured at small view angles. 

Computer simulations were conducted to determine the effect of deviations 
from the ideal antenna polarizations (required by the measurement technique) on the 
accurate recovering of all nine scattering parameters. The deviations, for example, 
can be introduced by the mis-match between surface and antenna polarizations pre- 
suming the scattering parameters are to be reported with respect to the surface polar- 
izations. Also, deviations obviously occur because ideal antenna polarization sped" 
ifications cannot be realized by practical antennas. Within these simulations a scatter- 
ing characteristic similar to that of the sea was employed as iHustrated in Figure (7,1)* 

All simulations were conducted with the assumption that the relative phase between the 
cross-polarizations was stationary across the main beam. 

The simulations indicated that there is little difficulty in rfscovering the three 
dominant scattering coefficients with off-the shelf antennas as illustrated by Figures 

(7.9) through (7.12) and Figure (7.15), Some difficulty can be anticipated when <)S^^^[^> 
is more than 10 dB beneath <|S^j^> as illustrated by Figures (7,11) and (7.12), In this 
case the cross polarized level must be better than 20 dB below the dominant (h) polar- 
ization. On the other hand, the antenna polarization requirement must be more care- 
fully observed when retrieving the six weak scattering parameters as illustrated by 
Figures (7.13) through (7.16). In some cases an adjustment in the relative phase of the 
cross polarization (if possible) can relax the antenna requirement. When the relative 
phase cannot be controlled in the case of cross polarized measurements, a rule of thumb 

for the quality of the antenna was established. If the measurement is to be performed 

2 2 
with a 0.5 dB accuracy and > is X dB beneath the geometric mean of <1S^[ > 

and<.|S^I^(^>, then the one-way cross polarized pattern must be X + 16 dB beneath 

the dominont. 

When the dynamic range of the scattering coefficients is large, it is 
clear from the simulation studies that the experimenter must carefully design his 
antenna to accurately retrieve the weaker coefficients. Certain types of antennas 
which have potential in achieving the Ideal polarization states are suggested In 
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chapter 6. The antenna specifications for observations over the sea or scenes 
having similar scattering characteristics can be based on the results reported in 
Chapter 7. However, for scenes having an entirely different characteristic, it 
is advisable to conduct simulations similar to those reported here. The simulation 
projjram, documented in Chapter 7 and Appendix D, may be easily modified 
for this purpose. These observations as we!) as others serve as a guide for designing 
meaningful and accurate scatterometer* experiments. 


The term scatterometer was introduced by R, K, Moore of the University of Kansas, 
A scatterometer is a radar designed to accurately measure the scattering properties 
of non“coherent scenes. The term scatterometer and non-coherent radar will be 
used interchangeably , 
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r. 


2.0 INTRODUCTION 


Various research programs have been proposed or are in progress to demonstrate 
the potential of monitoring, on a global basis, important geological, environmental, 
hydrological, oceanographic, meteoroloical, and agrarian parameters. The usefulness 
of remotely sensing certain parameters has been repeatedly demonstrated with optical 
and infrared sensors. More recently however, satellite and microwave technologies 
have developed to a point where microwave sensors are also suitable candidates as 
remote sensing devices. The microwave radiometer and radar scatterometer are prime 
candidate sensors. 

In remote sensing technology it is common knowledge that the retrieval of the 
remotely sensed parameters often entails compensation of the measurements for sensor 
and dtmosph'eric effects. The antenna is one element of the sensor system that requires 
special consideration. An understanding of the antenna-scene interaction is essential 
to designing meaningful experiments and for specifying the antenna with which the 
experiments are to be conducted. 

The radio astronomers, for example, have developed a rigorous theory involving 
the complex visiblity function to describe the interaction of o radiometer antenna with a 
small celestial scene, [I]. , Measurement techniques were based on the theory to derive 
complete emission properties of the scene. Recently Claassen and Fung [2l and Peake 
[3] have reported radiometer interaction relationships for nominally flat scenes having 
a simple partially polarized emission property, A measurement technique based on fhe 
relationship was developed by Claassen and Fung. Grody [4] has illustrated how the 
difference between antenna and surface polarizations imp ct radiometer experiments. 

To date little has been done to develop and use o comprehensive radar scatterometer 
antenna" scene interaction relationship for non” coherent targets. Most efforts have 
treated Dniy the spatial extent of the antenna pattern and have avoided general 
antenno and scene polarization properties [5], An exception occurs in the theory 
developed by Williams, et al . [6] . Their characterization of the scene parameters was 
based on the coherent radar equation and no rresurement technique was reported. 

In this study a complete non-coherent radar equation is derived and interpreted. 
The resulting expressions are valid for an arfsitrary antenna. The resuH is also extended 
to consider the differences between antenna and surface polarizations. The distinction 
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Is imporfanf when measuremeni'S are fo be compared with theoretical predicftonso It 
IS shown that six differential scattering coefficients are required to describe the antenna- 
scene interaction when reciprocity applies. Three of the six coefficients are complex 
valued. The coefficients are interpreted within the context of scattering and coherence 
theories, A retrieval technique based on intensity measurements is proposed to measure 
all the scattering coefficients. Computer simulations, based on the technique and a 
scattering characteristic similar to that of the sea, were conducted. The results of the 
simulation were employed (1) to validate an approximation used in the inversion, (2) 
to demonstrate antenna requirements for accurate retrieval of the scattering coefficients, 
and (3) to evaluate whether the distinction between ai;fenna and surface polarizations is 
important , 

The development of material in the subsequent chapters is accumulative. Chapter 
3 develops the background theoiy relevant to the derivation and understanding of the 
complete non-coherent radar equation. An adequate number of references are cited so 
that the reader can fill in background more deeply if he so desires. The derivation of 
the non*'coherenf radar equation is presented in Chapter 4, !n the latter section of this 
chapter the equation is altered to account for the difference between antenna and surface 
polarizations. Chapter 4 is strongly supported by the contents of Appendix A and B. 
Chapter 5 Is devoted to developing an understanding of the non~coherent radar equation 
and the polarization properties of radar returns. Certain scattering theories described 
in Appendices A and C are visited to illustrate the behavior of the complex valued co- 
efficients, The difference between antenna and surface polarizations is also illustrated. 
The measurement and inversion technique is presented in Chapter 6, The mathematical 
aspects of the inversion are treated in general and then specialized to the radar problem. 
Certain antenna properties which simplify the inversion are described. Antenna types 
capable of realizing these properties are suggested. The measurement and Inversion 
technique is evaluated within Chapter 7, A computer program which simulates the 
measurement and retrieval of all nine scattering parameters is described briefly. Full 
documentation of the scatterometer simulation program is provided in Appendix D, The 
results of the simulation are employed to illustrate antenna polarization requirements 
to measure all nine scattering parameters. Other practical aspects in making radar scatter* 
ometer measurements are also discussed. The measurement of the pattern amplitudes is 
specifically treated. Appendix E describes a computer program which specifies the points 


6 


at which a pattern must be measured. The conclusions and recommendations are presented 
in Chapter Bp A summary of all significant results is presented in Chapter 1 , It is 
advisable to read the summary before entering the technical chapters. 
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3.0 BACKGROUND 


3*1 In^^roduc^ion 

The theory and measurement of radar cross sections have been well developed 
for discrete coherent targets. An excellent review on the measurement of radar cross 
sections is found in a special issue of the Proceedings of the IEEE [7]. The theory of 
measuring non“coherent radar cross sections, except for the isolated works of Williams, 
et al. [6] and to some degree Hagfors [8], is largely lacking. Williams, et al, simply 
extended the theory for coherent targets to a non“coherent scene. In doing so, they 
over-looked some subtle distinctions between coherent and non“coherent theories as 
shown in Chapter 5. No measurement technique was presented. Hagfors, on the other 
hand, related Stoke's parameters for the incident wave to Stoke's parameters for the 
scattered wave in terms of the Mueller matrix [9]. In general, there are sixteen 
parameters in the Mueller matrix. However, as shown by Hagfors, targets exhibiting 
reciprocity and circular symmetry can be characterized by five independent entries 
in the Mueller matrix. Hagfors related his measurements to some of the five inde- 
pendent entries but no attempt was made to isolate all five entries. By using 
“Gedanken Experimented as Hagfors did, one can show that for a flat scene there 
are nine independent entries. At nadir there can conceivably be only five if the 
scene is cylindrical symmetric (isotropic). The fact that there ate nine independent 
entries in the Mueller matrix for flat scenes implies that there should be nine 
scattering parameters. To date only three scattering coefficients have been reported 
by the earth resources community [10] [11] [12], 

In preparing the background for this effort the author chooses to avoid the use 
of Stoke’s parameters and Mueller matrices since the earth resource community is, for 
the large part, unfamiliar with them. Instead polarization coherency matrices, an 
entirely equivalent representation for the polarization state of the transverse wave, are 
employed. The relationship between the entries in the coherency matrix and the standard 
differential scatterin;^ coefficients are clearer. To properly introduce the more general 
reception theory in terms of coherency matrices, the background for the reception of 
(polarized) monochromatic waves is first established. It is then employed to derive the 
coherent radar equation. In doing so the importance of reception theory in understanding 
the radar equotion is clarified. 
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3,2 The RecepHon of Monochromatic Waves ond the Radar EguaHon 


3.2 Transmitted Ffelds 

Schelkonoff has shown that the far field of any antenna has a dipoie field 
characteristic 113], viz.. 




2Xr 


(3-1) 


where 

N = radiation vector 

Z^“ intrinsic impedance of the medium 

X = wave length 

r = distance to the far field point 
k “ propagation constant 

and where the time factor is suppressed. The radiation vector in general, 

has complex components and induces a relative phase between the far field exponents. 
As a consequence, the far field has an arbitrary elliptical polarization. Now since 
N is proportional to the antenna input current I, Sinclair [14] proposed that a complex 
effective height vector L be introduced so that 


L(9,(>) = S(e,^.)/i 


(3-2) 


The far field can therefore be expressed as 

-dcoy 


t = 


(3-3) 


4'ht 


where 


w = radian frequency 

permeability of free space 

In general, L may have both 9 and J3 components in o spherical coordinate system and 
both may be complex. Specifically to emphosize this property, we may write t in 
normalized form 


t = 


|t| (cosfiT -t sin6e^*^i . ) 


(3-4) 


9 


at l•(m 1 pact'ly as 


c 


ifi 


(’v'f 




(3-5) 


iMlIm [s file orientation of a linear polarization when 0=0 and when 0?^ 0, 
(1 h llin lelative phase between the components* With a little effort 6 and 0 
c‘iM 1 mi inlated to the axial ratio and orientation of a polarization ellipse [15], 




Itru -pivinp Polarized Waves 

Suppose the above antenna is used to receive a plane wave described by 


t 


(Me 


+ E, 






(3-6) 


I 

whotp k Is the propagation vector* and E|^ are the vertically and horizontally 
polnrl/ml amplitudes, respectively. Then it can be shown by the reciprocity theorem 
[l6l tliut Iho open circuit vdtage induced into the terminals of an antenna having 
efloi.tlvo height^ is given by [14] 


V - t-t 


(3-7) 


Tlio |H»wnr available at the antenna terminals under matched conditions is given by 

P = /8R^ (3-8) 

wlmiei \{^ Is the radiation resistance of the antenna* 

M onochromatic Reception and the Radar Equation 

It has been shown by Sinclair [17] and by Kennaugh[I8] that a radar target 
- oon nut ns a polarization transformer. Sinclair expresses the transformation by a 
scnilui Inp matrix which can be incorporated in the radar equation. The scattering 
mnti hi h defined by 

S 






j? 

jP 


vv 


hv 




jpvh 


i? 


hh 


(3-9) 
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where 


o' 

pq 


= radar cross section for a q linearly polarized incident wave and a p 
reflected wave 


P 


pq 


= phase center for each component of the reflected wave 


p^q = V or h 

If the incident field is denoted as 


or 



4tt r 


-j kr 


(3-iO) 

(3"ii) 


then the scattered field ^ in component form is given by* 

s 




= l/(Vte r) S ) 

.®hs- 

L ®ht. 


(3-12) 


If1l^ denotes the complex effective height vector of the receiving antenna coHocated 
with the transmit antenna, then the power received under matched conditions is given 
by 



P = 

\% • 1 

t^|V8R^ . 

(3-13) 

or 

P = 

L. 


(3-14) 

where 

8(4tt) 

r"R^ 



^t,r 

= Itt.rl 


(3-15) 




■\t1 

.^htj 

I 



*Matrix notation for transverse wave components is frequently used throughout. 

II 


(3-15) 


’vr'' 


1 


hr 


Now It IS well known [18] that the antenna gain is given by 


G, 


I 9><|>) I 


t,r 




‘dJ2 


and the radiation resistance by 


Rt,r = 


L 


t ,r 


•dS2 


As a result, the received power can be written in more familiar form 


(3-!6) 


(3-17) 


W. 


x 2G.(0.(}>} G^(e,0)W.a 


(4¥ )"r‘‘ 

where the radar cross-section has been identified as 


(3-18) 


a = 



(3-19) 


The above expression for the radar cross section reduces to the linear polarized 
cases when both 1^ and contain a single non-zero component. For the coherent 
target the above formulation completely describes the interaction between three 
apertures, the transmitting and receiving antennas and the target*. 

Methods for measuring the elements of the scattering matrix have been re- 
viewed by Huynen [21], Methods of measuring radar cross sections ef have been 
reviewed by Blacksmith, et al, [22] and by Kell and Ross [23], 


3.3 The Non-Coherent Radar Equation 

Radar returns for a non "coherent scene have been defined in terms of a different 
ial scattering coefficient o-° rather than a scattering cross section o' . is unitless 

*It has been shown that the target actually acts like two coupled apertures [20], 
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and eKpresses ihe equivalent average rodar cross section per unit area. Moore [5] has 
shown from elementary considerations that when cr® is employed the average return 
is given by 


- W^/(4tt}^ j6^Gj,cr7r‘‘ dA (3-20) 

where the integration is performed over the illuminated area. For linear polarizations 
the differential scattering coefficients have been defined in analogy to ^ for the 
coherent case [5] [24] 


4Trr^< E (3“2I) 

^0 = ^ 

qp 

where 

A ’= illuminated area 
r = distance between the illuminated area and the point of observation 
= scattered field intensity 
■ E. = incident field intensity 

p,q = V or h 

Williams et al, [6] have shown that radar returns cannot be characterized by linear 
polarized scattering coefficients for an arbitrary antenna polarization and an arbitrary 
scene. They offer on expression for the differentia! power contribution by a small 
patch of the scene. Their formulation, however, is entirely identical to the radar 
equation for coherent targets, i,e,, the effects of spatial averaging have not been 
considered. 



3,4 The Reception of Quasi -Monochromatic, Portially Polarized Waves 
3,4,1 General 

The treatment of radar returns for non"coherenf scenes to date has relied 
on intuitive extensions of (polarized) monochromatic theory. Yet when one con- 
templates how the measurement of the non“coherent scattering coefficients is actually 
performed, one is acutely aware that the measurement involves estimating the mean of 
a fading signal having a certain doppler bandwidth. The resulting returns are, as a 
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consequence, quest -monochromatic rather than monochromatic. Furthermore, it is 
presumptive to anticipate that average returns from a randomly rough target are com- 
pletely polarized*. Indeed one should anticipate that the actual return will be a 
mixture of randomly polarized and polarized waves, i.e., will be partially polarized. 
Hence, any derivation of the non-coherent radar equation should include this possibility. 

Ko 115] has developed a comprehensive reception theory for quasi-monochromatic 
partially polarized waves. This theory is reviewed below and will be employed in the 
succeeding chapter. Important to this theory are the notions of an analytic signal as 
defined by Gabor [26] and the polarization coherence matrix as originated by Wiener 
[27] and Perrin [28] and later developed by Wolf [29], An excellent discussion of 
both topics appears in a text by Born and Wolf [30], 

3r4,2 Quasi -Monochromatic Partially Polorized Waves 

Wolf [29] has shown that a quasi-monochromatic wave whose bandwidth is 
small in comparison to the mean angular frequency ^ can be represented in analytic 
signal form of the type 

r(r,0,tj)st) = ey(r,e,(|),t)TQ + (3-22) 

where 

6v = (3-?3) 

The actual signal may be isolated by taking the real port of the above expression. 

The elements of this analytic signal have properties such that a L(r,t) >0 and 

V/ n 

ol u(r,t) is real . The correlation of the 0 and 0 components determines the state of 
n 

polarization of the wave. Wolf [29] defines the correlation by the complex factor 

< e > 

V h 

Mvh ° (3-24) 

/ < I e„ I ■‘x I a. I '> 

V n ’ 


’''Monochromatic waves are completely polarized. 


14 


Where the angular bracket > represents a time average. By Schwartz’ inequality, 

I i , The absolute value of P is a measure of the degree of correlation between 

V and h components while the phase angle of u^j^reflects the relative phase between 
the two components. If (y =1, the wave is said to be completely polarized. If 
^ ^ vh I ®v ®h wove is randomly polarized. The wave 

is said to be partially polarized when { y^j^ [ is between zero and one. The state of 
polarization may be completely characterized by a coherency matrix 


C ^ 


<e 


<e 


e*> 

V V 


p*> 

h®v 



(3-25) 


as shown by Wolf [29] (See also Born and Wolf [30]), 

Following Ko [25] and Collin [9] we may now suppose that a quasi-monochromatic 
partiolly polarizejl wave with coherency matrix C is incident on an antenna with 
effective height L, If the bandwidth of the wave or receiver is sufficiently narrow, 
then the open circuit voltage in analytic signal form is given by 


V = t*t(e,<i.,w) 


(3-26) 


where w is the mean frequency of the wave. If a coherency matrix is introduced 
for the antenna 


1 1* 
V V 


1 1 * 
'h'v 


1 1 * 

‘v‘h 


1 1 * 

‘h‘h 


(3-27) 


2 2 

where [1^1 + |lj^| = 1, then as shown by Ko [25], the power observed at the 

antenna terminals under matched conditions is givc.i by 


W(0.4>) = 


X^G(e,<J>) 

4tvZ. 


tr C^C 


(3-28) 


where tr is the trace operator andt is the transpose operator. The coherency matrix 
for the impinging wave is the transpose of that defined by Ko, All coherency matrices 
employed within this work are defined with respect to a coordinate system located at 
the observing antenna. Further interpretation of this expression is deferred until Chapter 
5 where a similar expression is discussed In the context of the scatteiom eter equation. 
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4.0 derivation of THE SCATTEROMETER EQUATION 


4.1 Introduction 

A generalized reception theory [91 [25] and notions from scattering theory are 
combined to derive the complete scatterometer equation for a scatterometer antenna 
having a specified but otherwise arbitrary transmit and receive property. The radar 
return is treated as a quasi “monochromatic “partially "polarized wave. The quasi" 
monochromatic character is induced into the return signal as the antenna linearily 
scans the scene. The scan is, of course, important in achieving a spatial average. 

The partially polarized assumption as well as the quasi “monochromatic characters per“ 
mits one to derive the scatterometer equation elegantly within the framework of the 
generalized reception theory. Intuitively, it is reasonable to assume that scatterometer 
returns are partially polarized since a spatial average constitutes the return. This inter" 
pretation will be illustrated in Section 5.4, 

The scatterometer equation is initially derived assuming that the scatterometer 
antenna transmission and reception properties are defined in terms of the surface polar" 
izations. In the last section of the chapter the distinction between antenna polarizations 
and surfoce polarizations is introduced and the impact of this distinction on the scatter" 
ometer equation Is shown. 


4.2 Derivation 


To determine the average power return from a homogeneous randomly extensive 
target, we suppose that a narrow beam scatterometer linearily scans across the scene 
with its antenna pointed in direction £2^ = (0^, If the scene has an anisotropic 
character it is important that j^be maintained constant during the scan {see Figure 4.1), 
The incident (transmitted) field may be related to the antenna complex effective 
height vector t^, a reception property. In the standard way [25]* 





j (ti3t-kr( t) ) 


(4-1) 


4‘irr 


* Reception and scattering relationships in the far field adapt well to the matrix 

notation. Capital letters will denote matrices and lower case letters will denote their 
elements. 
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X 



where is the vertically polarized component, is the horizontally polarized 


ht 


component ond 


1 


Vt 

ht 


(4-2) 


The subscript v and h are employed to denote the vector components aligning with 
the spherical polarized unit vectors iQond i^, respectively, associated with the sur- 
face coordinate system of Figure 4.1 , The backscattered field arriving with direction 
(df0) from a differential patch of the surface will be denoted 


Eg (e,(|),t) 



(4-3) 


Only transverse components for each line of sight (6,^ are admitted in the matrix. The 
field has thp units of volts/meter per steradian . Each component of must be regarded 
as an analytic signal since the relative motion between the antenna and the rough scene 
induces a time varying response for each line of sight. 

Now the antenna does not respond to the resultant field at the point of observation. 
Rather, if Lj^ denotes the complex effective height vector during reception, the antenna 
integrates the field components arriving with different directions so that the open circuit 
voltage appearing at the antenna terminals Is given by 






(4-4) 


where as the reader will recall, denotes the look direction and where the symbol 
T denotes the transpose operator. For narrow beam scatterometers the integration may 
be limited to the main beam and under worst circumstances to the first side lobes. The 
overage power observed at the terminals of the antenna under matched conditions is given 

<1 V(S2 ) I 

W(iIJ = 


8R. 


where R^ is the radiation resistance during reception (r) and < > denotes a time average 
or equivalently a spatial average since the scatterometer is scanning across the scene. 
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Expanded# the received power is given by 


W(%) 


tJf 


(4-6) 


Define a mutual (polarization) coherence matrix for the scattered fields as 


<ehs(n)e5,(fi')> <e|,s(!2)e^^(s2')> 


(4-7) 


Similarly a mutual coherence matrix^can be defined for the receiving antenna 




Then the average return can be written in compact form 


tr ' o' 


t-// 


tr M M dS2 d£2‘ 
r s 


(4-8) 


(4-9) 


where tr denotes the trace operator. 

For a rondom scene it Is reasonable to assume that the scattered fields are 
angularly non“coherent, i.e., 

<e.JSi) e1ja'}> = <e,. j£!)e*. fri)> 6(SJ-n') (4-10) 

1 S J S 13 J a 

The pragmatic aspect of this assumption is established in Appendix A. There it is shown 
that for a finitely conducting -smoothly undulating surface the degree of coherency 
(correlation) defined by 


“ij = < 6 . 5 ( 0 ) e* (a')> / < 6 ^ 5 ( 0 ) el^la) > (4-11) 
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is given by 


= 2 exp {-k^sin^6 a^A^d/Z) J1nc(k cosSR^Ae) 


where 


(4-12) 


0"^ = surface height variance 

k = 2 tfA 

0 = incident angle 
A0= small angular deviation from 0 
= radius of the illuminated area 
t, j = V or h 

The delta function type character of the angular coherency D.. is illustrated in Figure 
4„2 for a patch of rough surface having a radius of one meter and illuminated at 13,9 
GHz. A close examination of D.. reveals that, in general, the size of the illuminated 
area rather than the surface roughness dominates the correlation property at all angles 
of incidence except for the very large incident angles. The above result is based on 
plane wove illumination. The degree of coherency is thought to have a stronger delta 
function character in the case of spherical wave illumination since returns arriving 
from different directions arise from different patches of the scene whose statistical 
characteristics are poorly correlated, A discussion of this latter point within the context 
of a scattering theory appears in Appendix B. 

Under the above assumption the return power reduces to 






c 

r s 


(4-13) 


where 




(4-14) 


are the coherency matrices for the receiving antenna and the scattered fields, respectively, 

• 2 

As a result of the integration the units of the elements within become (v/m) per 

steradian. The change in units is clarified in Appendix B. 

Now also under the non~coherent assumption, it is permissible to introduce the 
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CORRELATION FACTOR 


1.0 
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ANGULAR SEPARATION (RADIANS x IG?) 

FIGURE A. 2 ANGULAR COHERENCY OF BACKSCATTER FOR 

A CIRCULAR PATCH NITH RADIUS OF ONE IlETER 
AilD FOR A FREODEICY OF 13.9 Ghz ' 
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notion of a matrix of differential scattering operators so that for each arrival direction 
the backscattered field (coming from a differential patch of the surface) Is related to 
the incident field in the following way 


% 


3, 


vv ■'vh 
'^hh 


(4-15) 


The second subscript indicates the polarization of the incident field and the first sub- 
script denotes the polarization of the resulting backscattered field. The objective for 
introducing this operator is that it identifies the scattered field components for each 
component of the incident field. With the introduction of this matrix the coherence 
matrix associated with the scattered field may be written as 


[Csl vv 





^^vh^vh^ht®ht^ 
'^'^vv^hv^vt^vt ^ 

<q Q-k a p* > 

^hh^hh^ht^ht 

[g:h 


2Re<2vvnh®vt®ht" ^ 

^'^vv^hh^vt^rit^ '** ^^vh'^hv®ht®vt^ 

(4-16) 


The action of the scattering operators on the incident fields is clearly evident in the 
above expression. 

If the incident wave were a plane wave, it is natural to define a scattering 
coefficient as 


?kl ^tt^Tt (4-17) 
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If the above definition Is employed for a spherical wave, the scattering operator will 
have to contend with a quadratic phase factor In the incident wave and with a varying 
intensity across the surface. The resulting scattering coefficient would depend on the 
geometry of the antenna pattern. However, under the non~coherent assumption the 
incident wave may be considered locally plane on each patch of the surface and the 
scattering action is then interpreted In accordance with the plane wave definition for 
the scattering coefficient. In particular, the expectations in may be written as 




^0 t 


4iTr 


a * 

jt^lt (4-18) 


where Equation (4"1) has been employed. The scattering coefficient is now allowed 
to vary with (6,J0) across the illuminated area. The integrand of the equation can now 
equation can now be written as ^ 


tr Cr Cg = 




(4TTr) 


o't^ 

0^1 


1 




vr 




+ 2 Re 


Ivrlhrhvv^hv" ^ 




2Re<S>,..S£u> 


^hv'"hh 'vt ht 


(4-19) 


As a result of the non~coherent assumption and the introduction of the scattering co- 
efficients, the transmit antenna pattern parameters have been dtvorsed from the composite 
scattering operators. The reader will observe that the scattering coefficient employed 
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2 2 

here has the units of m /m per steradlan. This definition is natural to this derivation 
and is a direct consequence of the integrating action of the antenna about its observe- 
tion point (Equation (4-4)). Further discussion of the scattering coefficients Is deferred 
until Section 5.2. The above steps in the derivation are clarified in the context of a 
simple scattering theory in Appendix B. As illustrated there, the above theory can be 
expressed as a continuum limit of an incremental theory which treats the backscatter on 
a patch by patch basis. Each patch is associated with an arrival direction. 

Now the following Identifications are helpful in re-formulating the results in 
more common terminology: 

o 

n 

(4-20) 




21 . 


max 

0,(j) 






= tan 


dm \il|;5/Re (4-21) 




tMyjl + 


(4-22) 




max [g‘. (0,(j>)} 


(4-23) 


-2 


- 2W,/R, 


(4-24) 




(Zq/ 4X ^) 


(^“25) 
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where ! = t (transmit) or r (receive) and p = v or h polarization. Descriptively, during 
transmission g is the normalized gain of the vertically polarized pattern whereas g^^ 
is the accompanying horizontally polarized pottern. The relative phase between these 
two polarizations is denoted os 0^. In general oil three are functions of the pattern 
coordinates. G‘ is the gain under a matched polarization condition and is the 
maximum gain (presumably on boreslght) . is the transmitted power and is the 
radiation resistance when the ontenna is transmitting. Similar explanations apply to 
the reception parameters. They are Identified with a subscript r. With the introduction 
of the above pattern parameters, the scatterometer equation can be written os. 


w(n^) = (x/'iit) 


dn 


( 4 - 26 ) 


where 


. Hr 


•f 2/ g„z.ght ^ 


^ht 


vt^ht hV 

2 . 


(4-27) 


•U is interesting to note at this point that there are ten scattering coefficients. Addition' 
al simplification occurs when reciprocity applies. Under this assumption 

\h (4-28) 
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since field reciprocity implies that the operators must be identical. When the above 
property is applied to the definition of a scattering coefficient 


= g g x<|S \^> -t 2Re(g / g + 

tr ^vr^vt * vv’ '^^vr ^vt-'ht 


2 


Svr^hr^ ^ ^vt^ht' 


Sht'' 3vr9hr®^^'’>^^h5hh^ g^rShr’ 

^vt^ht ^ t r) + + ^hr^vt * 

2/ g„^gh,.3„t9u,. Re <lS„.|^> + 


'vr^hr^vt^ht 
® hr^ht*^ ^ ^hh ^ ^ 


(4-29) 


When reciprocity applies the number of coefficients reduces to six. 

The above result is the complete non-coherent radar equation under the re- 
ciprocity assumption. Although the equation was derived from the viewpoint of polar- 
izations ascribable to the surface, the same equation would have resulted had the 
antenna and surface polarization states been defined with respect to the antenna. In 
the latter case the scattering coefficients would not be comparable with those defined 
by the theorist who derives scattering coefficients with respect to the surface polari- 
zations. In addition, the scattering coefficients for an arbitrarily line of sight would, 
in general, be a function of antenna view angle also. Pragmatical ly,the antenna polar- 
izations are referenced to a coordinate system rigidly bound to the physical antenna. 

The antenna polarization vectors, consequently, move with the antenna as it changes 
view angle. The surface polarization vectors on the otherhand, remain rigidly oriented 
with respect to the surface. The transformation between the two polarizations description 
is derived in. the succeeding section. The distinction between antenna and surface polar- 
izations on the scatterometer equation is treated simply by transforming the transmission 
and reception coherency matrices, and C^, from the antenna coordinate system in 
which they were measured to ihe surface coordinate system in which the surface polar- 
izations are naturally defined. 
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4.3 The Scatterometer Equation lnciudin<:? the Distinction Between Antenno and Surface 
PolanzQfrions 

Suppose that the antenna patterns, both polarized and cross-polarized patterns, 
are measured with the scatterometer antenna mounted on an azimuth -over-elevation 
positioner. To describe the antenna polarizations measured from such an antenna 
positioner, afix a primed coordinate system rigidly to the antenna. Let the x' axis 
denote the boresight axis and let the z' axis be oriented In a direction coinciding with 
the vertical polarization sense (with respect to the antenna) for an observer on the bore- 
sight axis. Then the antenna polarizations, vertical and horizontal, will coincide with 
the spherical polar unit vectors ig, and , respectively, of the afixed coordinate 
system. The antenna coordinate system is illustrated with respect to the pattern measuring 
antennas in Figure 4.3. Patterns are "cut" by incrementing the positioner in elevation 
when the y* and y" axis coincide and then rotating the positioner about the z‘ axis. The 
measuring antennas are located on the x" axis of the range coordinates (x",y",z"). 

Within the antenna coordinate system the transmitted fields will be denoted by 
e 0 i| and e^,^ and the complex effective reception heights by Ig,^ and Both 

pairs of parameters are, in general, complex (to convey the relative phase between 
members within the pairs) and vary with 0 ' and 0 \ 

Now locate the antenna (primed) coordinate system so that its origin coincides 
with the origin of the surface coordinate system (Figure 4,1). Without loss in generality 
it is assumed that the antenna scans linearily in the x direction of the surface coordinate 
system and that observations are conducted in the x - z plane. The antenna is so oriented 
that Its vertical polarization sense coincides with the surface vertical polarization sense 
at the intersection of the boresight point with the surface. Within the xz plane the antenna 
is pointed at an angle 9^ with respect to the local vertical (z axis). The geometry of 
the two coordinate systems relative to one another is shown in Figure 4.4, 

To develop the relationship between the antenna and surface coordinates consider 
any line of sight vector which emanates from the common origin and whose extension 
intersects the surface (See Figure 4.4). By definition, the antenna polarization pair 
(Tg*, 1^’) and the surface polarization pair (Tg/l^) are both perpendicular to T^, It 
follows that the polarization pairs at every line of sight are related by a simple rotation, 

* Not normalized as in Chapter 3. 
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FIGURE COMPARISON OF ANTENNA AND SURFACE 
COORDINATE FRAMES WITH THE SURFACE POLARIZATIONS 
DEFINED WITH RESPECT TO A GENERAL LINE OF SIGHT 
VECTOR 


say 4 >. Define so that* 




r , = cos 

<P 


and 


(4-30) 


^9 • • ^ 6 ' = 


(4-31) 


By noting the transformation between the coordinate systems, the reader can easily 
show that 

cos 1 |J = cos (j> cos 4>' + sin d sin 4)' sin 0^ (4-32) 


and " 

sin 'f' = cos 0 ( sin ^ cos <j)' - cos c|) sin 4i' sin 9^) + sin 0 cos 0^ sin ^ 

(4-33) 

where 0 ' is the spherical azimuthal angle in the primed coordinate system. Now 0 ' 
can be eliminated by observing that 


tan If' = { • iy. ) / ( • i,,, ) 


(4-34) 


to get 


4)' = tan 


-1 


sin 9 sin 4i 

cos 0 cos 0Q + sin 0 sin cos 4> J 


(4-35) 


Finally from the above we have established the transform T between the antenna and 
surface polarizations, viz.. 


. ^ Ui 


(4-36) 


where 


T = 


cos i|; sin tj; 

-sin 4^ cos t|j 


(4-37) 


* Note: An alternate method of mounting the antenna could have resulted in defining 
4* so that cos^t' = ig etc. The difference between the two is discussed in Chapter 
5* 30 


The entries in T are provided b/ Equations {4-32) and (4-33) with the assist of 
Equation (4-35). When the antenna pattern Is finally introduced the following 
relationship 


COS0 


1 . 


'z' 


(4-38) 


or 


cos0‘ == cos0sin9^ - si n0cos0^cos(|) ( 4 - 39 ) 

in addition to Equation (4-35) will be helpful in identifying the pattern coordinates 
when the surface coordinates are given. 

Now from the preceding derivation (Equation (4-13)), we had 


where 


k fk 


'^tr ^ 8^ I 4 


(4-40) 
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= > 


(4-41). 
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(4-42) 


'vt 


^ht ^vt 


^ht 


'ht 


(4-43) 


The above coherency matrices are written in terms of the surface polarizations 
since the scattering operators are defined on the basis of these polarizations. When 
the antenna transmission and reception properties are, however, defined within 
another coordinate system (the primed coordinate system), these properties must be 
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appropriately transformed. It is easily shown that if C'^ and C'^. are the coherency 
matrices in the primed coordinate system, then in the surface coordinate system 


and 


C = T C T 

r r 


t A 

C* = T C, f 


(4-44) 


(4-45) 


't ' “t 

From the above expressions the following identities can be established 

hvrl^ = cos^ + sin 2^ Re l’^ + sin^ if. 


* 

vr hr 


n - iCi^ 


vr' 


sin ^ cos t|j + cps ijj 1 ^^1 j!|* 


2 I* t 

- sin I- ’vr 


* 

^hr ^vh 


-[ 


1 1 
vr 


* ~ 

hr_ 


(4-46) 


1 


hr 


= sin^ ^ Pvr|^ “ ^vr V Phrj^ 


I V ■ I T * ... I I 

A similar expression can be established for the elements of C^., It is noted that th 
coherency matrices reduce to those In the surface coordinate system when ^ = 0, 
Now let g0t|.i> 9j^i-|.r describe the antenna during transmission and ggij.j 
gj^ij. ond /3^, during reception. When the transformed coherency matrices are in- 
corporated into the scatterometer equation and relationships of the type as shown I 
Equations (4“20) through (4-25) are noted in the antenna coordinate system, the 

scatterometer equation can be written as 

9 


»tr < ) = 




( 4it )‘ 


/ 


tr 




(4-47) 


V / 

provided that the following identities are understood • 

rv 


Syp " 3g.p cos^ + + •'Sgip g^,.p sin 2<i, cos e'p + g^,p sin^ * 


9hp 


= 9e'p 1' - ■'9 0'p %'p ^ ^'p %'p 


COS^ l(> 
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S, = tan 
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V3^‘0 
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where p = i- or r. The latter identities indicate how the common antenna parameters 
transform. It is noted that /3 , g or gj^p is each dependent on all three antenna 
porametett, , g'g °nd g-^p. To appreciate the additionol complexity in the 
scatteromei-er equation resulting from the transformation expand the integrand in the 


form 
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Then it will be noted that 
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When accurate measurements of the scattering coefficients, say the complete 
set of nine parameters is desired, one must contend with inverting a system of integral 
equations of the type derived above „ The scattering coefficients are rigorously defined 
in terms of the surface polarization, a definition universally employed by the scattering 
theorist. If comparisons with theory are necessary then the antenna properties must be 
transformed to conform with this definition. To date measurements have been reported 
without the recognition that Equations (4-47) through (4-50) governor the interaction 
between the scatterometer antenna and the scene. Yet reasonable agreement between 
measurements from targets with known statistics and theory have been reported flO] 

[31] for the polarized scattering coefficients. This indicates that the complexity of 
through may be avoidable under some circumstances. To resolve this problem and 
related ones,the polarization coordinate systems will be compared and the character of 
the scatterometer equation will also be examined in depth in succeeding chapters. 

Once the character of the scatterometer equation is established, a measurement technique 
to recover all six scatting coefficients is specified. Computer simulations based on 
the specified technique are then conducted to determine antenna requirements for 
accurate measurements. 
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5.0 DISCUSSION OF THE 5CATTEROMETER EQUATION 


5,1 IntroducHon 

This chapter is devoted to developing on understanding of the scatterometer 
equation. The character of the scattering coefficients Ts established by reference to 
previous definitions, both coherent and non-coherent. It is shown that the non" 
coherent definition appearing in the literature must be extended to include new kinds 
of coefficients. The composition of the average return is examined from the standpoint 
of coherence theory and the complete set of scattering coefficients. The importance 
of the phase characteristic of the wave and the receiving antenna in governing the 
observed power is described. It is also shown that certain properties of the coherent 
scattering coefficients cannot be extrapolated to the non~ coherent case. Well known 
theories applicable to the sea are also employed to illustrate the behavior of the 
scattering coefficients having a cross" correlation property. Other possibilities for the 
cross-correlation coefficients are also treated intuitively. 

Within this chapter it is also shown that this formulation of the scatterometer 
equation admits partially polarized returns. A previous formulation |6] failed in this 
respect. The degree of polarization of the average sea return is specifically illustrated 
using a simple scattering theory o 

Finally the distinction between surface and antenna polarizations is illustrated. 
Certain aspects of this distinction ore qualitaively applied to specifying antenna 
requirements. 


5,2 General 


5.2ol The Scattering Coefficient 

The scattering coefficients within the scatterometer equation may be partially 
identified with the differential scattering coefficients defined by Peake [241. As the 
reader may recall,, Peake defines 

47rr2 <|e. 
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A cos^e.^1 
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where |e.^| is the polarized incident intensity>< (e.g) > is the i polarized backscatter 

intensity in volts^/m^, 6 is the incident angle and A is the illuminated area. The 

scattering coefficient employed in this formulation is simply related to Y.j in the 
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following way 


c|s,,f> = i.yiV (5-2) 

The difference by 4it occurs since the scattered intensities were defined in terms of 
inverse steradians. It is clear that, in view of the three addition coefficients, it is 
more appropriate to define the coefficients in terms of the scattering operators 

(5~3a) 

where Ej^ yields a scattered field with units volts/meter. The operators in the 
derivation are related to those in the definition in the following way 

■ . =■ AA cos (S-3b) 


To understand the function of these cross"corre lotion coefficients one must 
examine the coherence matrix for the scattered wave. Under the non“coherence 
and reciprocity assumptions the elements of are given by 
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The power in the scattered wave is carried in the trace of whereas the relative 
phase between the orthogonal components in an average sense is carried in the off 
diagonal elements. From the structure of the coherence matrix It Is evident that the 
cross-correlation coefficients can be complex valued. The cross-correlation co- 
efficients therefore alter the phase property of the scattered wave. Some of the re- 
lative phase is attributable to the incident wave and some to the surface, eg,, e^^ 

®ht* °"^<^v\?hh*> ' 
diagonal terms and consequently contribute to the total power available in the scattered 

wave when both polarizations appear in the incident wave. 

During reception the coherency matrix of the antenna interacts with the co- 
herency matrix for the wave. The interaction is completely described by taking the trace 
of C C ‘ , The trace is given by 

r $ yj t ^ 


respectively. The cross- correlation terms also appear in the 
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(5-5) 


(The expanded version of the trace is given in Equation (4-19) of Chapter 4). The phase 
interaction between the scattered wave and the antenna is described by the middle terms 
in the trace expression. These terms are complex conjugate pairs and consequently make 
a real contribution to the observed power. 

To show that the phase properties of the antenna and the wave are important to 
the observed return it must be recalled that in the case of polarized waves, the antenna 
polarization states must be matched to the polarization state of the wave to observe 
maximum power [14] , This requirement in terms of coherency matrices implies that 


C^/trC^ 


= C^/tr 


(5-6) 


Under the polarized assumption takes the form 


C^/trC^ 


(5-7) 


where a + b = 1 , c = and a is the relative phase between the v and h 

2 

components. The observed power will be proportional to (a■^ b) trC^ tr or 
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frC frC . If the reception matrix had been given by 

fb -cl 


no power would be observed at the antenna terminals as can be easily demonstrated. 

In this case the antenna polarization state is said to be orthogonal to the polarization 
state of the arriving wave . 

If the wave is partially polarized Ko [25] has shown that the observed power 
may vary from a minimum of ”P)trC^/8 to a maximum of +P) trC Z^. 

where P is the degree of polarization (See Section 4 of this chapter). To understand 
this result it must be noted that a partially polarized wave can be uniquely decomposed 
into a sum of a randomly polarized wave and a completely polarized wave [30] , As a 
consequence for an arbitrarily polarized backscattered wave,J’he decomposition can 
be written as 
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where P is defined in Equation (5- 30) and 
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The first term is the randomly polarized component and the second is the completely 
polarized component. If the receiving antenna is orthogonal to the completely polarized 
part, then only the randomly polarized component is observed at the antenna terminals. 

If the antenna is matched to the completely polarized part maximum power is observed. 
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The extremes in the observable pov/er are a positive indication of the importance of 
the phase interaction of the antenna and the wave. The cross“correiatlon scattering 
-coefficients and the cross-polarized scattering coefficient can be effective in 
altering the phase property of the return. 

The cross-correlation terms have their analogues in scattering theory for coherent 
targets [21], In the theory for discrete targets the complex scattering matrix is commonly 
employed to define scattering properties. The elements of this matrix have the 
property that 




(5-11) 


However for a statistical target this property is not necessarily true. Since the scatter- 
ing coefficients can be considered as a inner product of the form 






where it is crincluded by Schwartz* inequality that 
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(5-12) 


(5-13) 


As a consequence the magnitudes of the scattering coefficients may not be simply re- 
lated as suggested by William^ et al . [6] . The inequality is an admission that the 
amplitudes or phase centers between scattered field component^ can be correlated. 

One can identify two scattering parameters with each complex valued scattering 
coefficient, viz,, its real and imaginary parts. As a result one may attribute nine 
scattering parameters to Equation (4-29) where reciprocity has been applied. Similarly 
from Equation (4-27) where reciprocity has not been applied, sixteen scattering para- 
meters can be identified. These observations are in agreement with the "Gedanken 
Experimented' cited in Chapter 3. 

5.2 2 Special Cases 

An examination of the scatterometer equation indicates that the equation under 
appropriate conditions reduces to the classical cases. For example, when vertically 
polarized measurements are conducted, i,e,, g^^ = g^^ = 0, the integrand factor of 
the scatterometer equation becomes 


®vr®vt ' ^vv' 


(5”14) 
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Similarly when horizontally polarized measurements are conducted, i.e., g . = g = 0, 


I 




vt '^vr 

(5-15) 


and when cross“polarized measurements are conducted, i.e,, = 0, 
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tr 


^vt^hr 




(5-16) 


It should be noted that the reductions result from highly idealized representations of 
practical antennas. Invariably antennas have cross polarized leakage; and when leakage 
is present other scattering coefficients, both auto-correlation and cross-correlation 
types, will be excited. As shown in the last section of Chapter 4, even if the leakage is 
not present, the difference between antenna and surface polarizations can introduce, 
in effect, cross-polarized components in the incident wave and in the reception antenna. 
The impact of undesirable antenna properties and polarization mis“match on the 
measurement of isolated scattering parameters will be treated in Chapter 7. 

An understanding of the cross -correlation coefficients from theories applicable 
to sea returns is developed in the succeeding section. 

5.3 Characteristics of the Correlation Terms 

Several scattering theories are examined to disclose the. character of the cross- 
correlation terms in the scatterometer equation. Approximate backscatter solutions to 
the small perturbation theory [32] [33] [34] [35] and the Kirchhoff theory [36] [37] are 
specifically examined. These theories with some slight alterations are thought to apply 
to ocean backscatter and have shown reasonable agreement with measured results. The 
selection of these theories by no means exhausts the possibilities. Physically intuitive 
arguments ore given at the end of this section to further enhance our understanding. 

When returns are considered from a surface having a small roughness, satisfying 

2 2 2 

k m cos 0 £ 1 , where o' is the rms surface height and k is the wave number, it can be 
shown that (see Appendix C)* 


<S S. u * > 

vv hh 
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- COS 2k sine, 0 ) (5"17) 


* The reader should be aware that in constructing the scattering coefficients from 
scattering theory and identifying them with measured coefficients involves an 
ergodtc assumption, i.e., an ensemble average is equated with a spatial average, 
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In the above equation 
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cose + ^ sin^e 

Comparison of the magnitude of this term with the polarized scattering coefficients 
diows that 

“ V (5-20) 

The equality is true at least to the order to which these solutions are valid. The 
magnitude of this term is illustrated in Figure 5,1 wherein it Is also compared with the 
polarized coefficients. The computations were based on a slightly rough sea. The 
phase of the cross correlation, defined by 

*= tan-l (lm<S^„S*|,>/ Rs<S^^S*|^>) (5-21) 

was computed and is shown in Figure 5,2 for three different water temperatures. The 
sea water temperature alters the complex dielectric constant of the surface and con- 
sequently influences Rj^. It is observed from the graph that the imaginary part of 

<5^^ ^hh*^ small in comparison to its real part and only tends to become significant 
at the larger angles when compared with the real part. 

An examination of the integrand of the scatterometer equation (4-29) indicates 
that the above cross-correlation term con make a significant contribution to a radar return 
when like and cross antenna polarizations are present during transmission and reception. 
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FIGURE 5.2 CROSS-CORRELATION PHASE PROPERTY BASED ON SMALL PERTURBATION 





The conlribuHon can be positive or negative depending upon the value of (^^, + and 
can be comparable to the sum of the contributions arising from the polarized scattering 
coefficients. To illustrate the above statement it is sufficient to observe that the polarized 
contributions are proportional to > and ^hh^ ^ other- 


hand, if the imaginary part of <S is smai 
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< I ^ contribution by the cross-correlation coefficient is given by 
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If right circular polarization* is transmitted and received, then cos = “I and 

g^j. = gj^^ = g^P = g|^^* !<■ is apparent that v\/hen <| , the magnitude 

of the cross-correlation contribution is Identical to the sum of the polarized terms. 

The sign of the contribution is, in this case, negative. However, had the wave been 

received with a LC polarized antenna, the sign of the contribution would have been 

positive. The contribution by this scattering coefficient can also be very effective 

2 

when attempting measurement of a weak scattering coefficient such as with 

a "tinearily" polarized antenna having some cross polarized leakage. This will be 
illustrated in Chapter 7, 

When the cross -correlation terms of the type and ^hh*^ 

are examined In the context of small perturbation theory, it Is easily shown that these 

2 2 2 

coefficients vanish at the lowest order where <IS 1 > , <■! S, ■ I > and ^IS i 1 >■ 

vv ' hh vh' 

are non -zero (see Appendix C)*t The lack of correlation is physically reasonable since 
it is believed that the cross-polarized fields result. from multiple scatter. When higher 
order solutions are included these cross-correlction terms will not vanish; however, 
their magnitudes will be extremely smoll. 

The above theory is thought to apply with some modification to the sea for angles 
of observation between 30 and 80 degrees [37] . At smaller angles Kirchhoff theory [371 
has predicted sea returns reasonably well. When the theory reported by Fung [36] is 
employed to explain near vertical returns it can be shown that (see Appendix A) for 
an isotropic stationary gaussian surface 
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Circular polarization is defined with respect to the antenna; but for narrow beams and 
sufficiently large angles, the circular polarization state transforms to the surface with- 
out significant alteration. This will be clarified at the end of this chapter and within 
Chapter 7. 


** Recall 
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order solution. 
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and 


are first order solutions and <■/ S^|^l > is a second 
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where 
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Comparison of the magnitude of this cross-correlation terms with the magnitudes of 
the polarized scattering coefficients again shows that (Appendix A) 

l<Svv=hh*>l = V (5-28) 

The equality is valid to at least first order in corrections to the reflection coefficient for 
the local slope. The magnitude and phase of the above result is illustrated in Figures 
5,3 and 5,4, respectively, for an isotropically rough sea surface having a moderate 
rms surface slope. It is noted that the phase property may be attributed to the linear 
corrections of the reflection coefficients for the local slope. The resulting reflection 
coefficients’ compare favorably with that for normal incidence. 

When cross-correlations involving or are considered within the Kirchhoff 
approximation little can be said regarding their character. A typical cross-correlation 
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between field componenl-s is given by 
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(5-29) 


where integration by parts has simplified and the reflection coefficient has 
been linearily approximated. Since the expectation involves higher order slope 
terms it is anticipated that the correlation will be weak. The stationary phase tech- 
nique for solving the integral, for example, would cause the integrand to vanish. 

On the basis of the above simple scattering theories, it is clear that the cross- 

correlation coefficient <S Si i *> can contribute to a radar return when both like 

vv hh 

and cross polarizations are present during transmission and reception. The theory for 
the slightly rough surface Indicates that the phase of the correlation product is depend- 
ent on the relative phase between the so-called Rice reflection coefficients. The 
phase factor Is somewhat significant at the larger angles. Correlation products contain- 
ing or vanish for the slightly rough surface and appear to be negligible for 
Kirchhoff type surfaces also. These two theories by no means exhaust the possibilities. 
Consider, for example, radar returns from a strongly de-polarizing scene such 
as a tenuous vegetated terrain in which the depolarization is attributable to linear re- 
radiation, Intuitively, one would expect sizeable contributions from the cross-cor- 
relotion products containing or The correlation contributions will likely arise 

from a single scatter process, particularly at the canopy. The contributions, as an 
examination of C shows, will arise if like and cross antenna polarizations are present 
during transmission op reception or both. When / S i I > is somewhat less than 

2 9 vn 2 

<1S^^^1> and and one attempts retrieval of with a "linearily” 

polarized antenna possessing a cross leakage pattern one can anticipate contamination 
notonlyby <jSy^l> , and but also by ^ Sj^/> and 

^ ^vh ^hh*> * 
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5.4 The Degree of PoiorizaHon of Radar Returns 


Another important aspect of the scatterometer equotion derived here is that the 
backscattered fields can be considered partially polarized. The partially polarized 
character is induced by measuring average returns from a non**coherent scene. This 
permits us to consider a statistical coherency matrix as a suitable representation for 
the return. It is well known that the degree of polarization of a wave is given by [30] 
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1 - 


(tr O' 


(5-30) 


where Cgis the coherency matrix of the wave and IlCgildenotes its determinant. If 
P = 1 the wave is said to be completely polarized. This occurs If and only if II Cj/= 0. 
To show that the present formulation admits partially polarized waves {possibly randomly 
polarized, P = 0), it is sufficient to show that HC^I! ^ 0, Now is given in Equation 
(4“16) and with a little tedious effort one can show 
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where the subscript t has been dropped. If the determinant is to vanish independent of 
the transmitted fields then each difference term in the above expression must vanish. 
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A non’^taHstical target having a scattering matrix 
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(5-32) 


will obviously meet the requirement. However,, the general result Indicates that the 
backscattered wave is partially polarized as the examples below illustrate. 

For the case where terms of the type ^vh*^ assumed negligibly small, 
as we suspect they are over the sea, we have 




(5-33) 


At moderate to large incident angles over the ocean it is anticipated that the above 
term can be significantly different from zero* Specifically iPa horizontally polarized 
wave is transmitted, we have 



(5-34) 


The corresponding degree of polarization is given by 


h. 


(5-35) 


<l^hhl <ISvhl > 

Similarly when a vertically polarized wave is transmitted, the degree of polarization 
is given by 




(5-36) 
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If a circularly polarized wave is transmitted, then 



<|SJ^>[<|Svvl^> ^-<ISKhl^>* 


2R„<S S. *> 
e vv hh 





(5-37) 


and the degree of polarization is given by Equation (5-30), 

The above cases were evaluated as a function of incident angle for the scatter- 
ing characteristic of Figures 5,1 and 5.2 , The results are shown In Figure 5,5 , 

It is apparent that the partially polarized character is an Important factor when the 
non"coherent scatterometer equation is appropriately Interpreted. 


5.5 Visualization of the Polarization Properties of the Antenna and Scene 

Within the latter section of Chapter 4 the scatterometer equation accounting 
for the difference between antenna and surface polarizations was derived. It was shown 
that the polarization mis-alignment could be characterized by a simple rotation of either 
orthogonal polarization pair through an angle 1 ^. To show this mis-alignment character, 
rather than study the functional behavior of on (0, 0 ), it is more convenient to fall 
back on the properties of the spherical polar vectors i0 and y. 

Regardless of whether one considers the antenna or surface coordinate system, 
the projection of the polar vector i0 and the azimuthal vector on any sphere whose 
center is located at the point of observation can be depicted,respectively, by longitudes 
and latitudes on that sphere. For any line of sight emanating from the origin of the 
sphere the longitude and latitude lines intersecting the line of sight on the sphere will 
correspond to the orientation of vertical and horizontal polarization, respectively, for 
that line of sight. We can therefore employ spheres marked with longitudes and latitudes 
to visualize the antenna or the surface polarizations. 

To compare the alignment between antenna and surface polarizations choose 
the radii of both polarization spheres so that the spheres are tangent to the scattering 
surface at the sub“observational point as Illustrated in Figures 5.6 and 5.7. A pole of 
the surface polarization sphere will be afixed to the sub-observationol point. This 
polar axis will correspond to the z axis of Figure 4,4 The xy plane coincides with the 
equatorial plane and is parallel the surface. The antenna boreslght axis lies in the xz 
plane and points at an angle of 0^ with respect to the z axis (Figure 5.6). Now, on 
the otherhand, the equatorial plane of the antenna polarization sphere coincides with the 
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plane containing the boresight and the y axis. This plane corresponds to the x' y' 
plane in the antenna coordinate frame (Figure 5.7). The polar axis of antenna 
polarization sphere aligns with the z' coordinate. Comparison of the orientations of 
the latitudes and longitudes for any common line of sight will indicate the polar- 
ization mis*"alignment property (Compare Figures 5.6 and 5.7), Within the plane 
of observation (the xz or x'z' plane) regardless of view angles the polarizations co" 
incide. For any other line of sight there will be a difference In alignment. The mis- 
alignment is greatest in the polar regions of th*' antenna or surface polarization 
spheres. When the antenna is pointed toward the horizon the alignment is everywhere 
perfect (one must mentally rotate the sphere in Figure 5.7 so that the x* axis points 
to the horizon). When the antenna Is pointed at the sub-observational point, the 
misalignment is severe everywhere in the vicinity of the sub-observational point. 

Within the nadir region the scattering coefficients defined with respect to the surface 
as compared to those one may define with respect to the antenna differ radically. 

For example, if a significant anisotropic scattering behavior occurs at nadir, any 
finite beam scatterometer would tend to integrate this behavior. The measurement, 
as' a consequence, would be difficult to refer to the surface polarizations. The surface 
polarization character at nadir indicates that infinitesimai becmwidths must be used if 
the nadir region is to be probed and if scattering coefficients defined with respect to the 
surface are to be reported. This is clearly true if there is a difference in <jS^j > and 
<[Sh^f^> scattering properties as viewed with respect to the surface. 

As pointed out in Section 4.3 , there is an alternate method of mounting the 
antenna which will produce a different polarization character. Suppose the antenna had 
been mounted so that its horizontal polarization vector (i^i) op the boresight axis (x‘) 
aligned with the surface vertical polarization at that line of sight. The polar axis of 
the antenna polarization sphere (z*) would coincide with the -y axis of the surface co- 
ordinate system. The corresponding polarization sphere is illustrated in Figure 5.8, 
Comparison of the polarization property with that of the surface indicates that the mis- 
alignment is invariant with view angle and the polarizations do not align globally for 
any view angle. The polarizations continue to align in the plane of observation; however, 
the same mis- alignment In the nadir region remains a problem. 

Regardless of the mounting position tt is evident that for non-zero beomwidth 
antennas the discrepancy between antenna and surface polarizations prevails in the 
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nadir region. It is clear that nadir is a forbidden region when one views it from 
the point of surface polarizations. At all but small view angles the polarization 
discrepancy over the main beam of narrow beamwidth antennas is generally small 
(how small will be shown in Chapter 7) , At these angles as the beamwidth becomes 
narrower for a linearily polarized antenna, the percent antenna power occurring In 
the orthogonal surface polarization becomes smaller. As one approaches nadir 
the beamwidth must become increasingly narrower for the same degree of mis”olign- 
ment. 
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6.0 THE INVERSION OF SCATTEROMETER MEASUREMENTS FOR THE SCENE PARAMETERS 


6,1 Introduction 

The recovery of the scene scattering parameters entails an appropriate set of 
measurements and the inversion of a corresponding system of scatterometer equations 
of the type derived in Chapter 4. Within this chapter a measurement and an inversion 
technique is derived to recover a complete set of scattering coefficients. The technique 
is also specialized to the case where the scattered fields may be regarded as completely 
polarized. To assure that the technique Is both as simple and as accurate as possible^ 
certain antenna properties are specified. The consequences of not employing a suitably 
chosen antenna is illustrated in Chaper 7, The mathematical and physical aspects of 
inverting scatterometer measurements are treated in the following sections. Certain 
antenna properties which are helpful in approximating the measurements by a system 
of olgebraic. equations are identified. In this chapter the distinction between surface 
and antenna polarization Is appropriately discarded to simplify the presentation. 

The consequence of this action is treated in Chapter 7, 


6.2 Mathematical and Physical Aspects 

The inversion of scatterometer measurements falls into the same mathematical 
category as do many remote sensing problems. Typically, the observational relation” 
ship reduces to solving a Fredholm integral equation of the first kind, viz,. 


g(y) = f K(y,x) f(x)dx 


(6-1) 


where K(y,x) Is usually a continuous function over a rectangular domain attributable 
to a sensor, f(x) is the unknown sensor stimulus and g(y) is the observed sensor response. 
The scatterometer equation .s a generalization of the above expression. Since there 
are nine unknown scattering parameters, it is clear that there must be at least nine 
different kinds of measurements to retrieve all the parameters. If each kind of measure- 
ment is identified by a subscript i and if the scattering parameters are denoted by 
^j(n) where « ( 9 , tj) ) , then the system of measurements can be written as 
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where 


«J “ 


i = 1,2, ... ,9 


(6-2) 


^11 ~ 

^12 = (Shr9i,t)i/''' 

1^13 = (9vt9hr^3vr9ht^2 •'g^tSht^vr^hr^'’^ < ’ i ' 

S-4 = 2(’'9vt9ht9vrV‘=° = <®t^®r>'i/’'’‘ 

= -2(’'9vt3ht9vr9hr‘='’=' 

Ss = 2(9v/ 9vtSht9°=V9vt'' 9vr9hr'''’=^>1''’'' 

3vrShr“=®r>/''“ 

l^iS “ ^f9hr'^V^t'^°=V3ht’' 9vr9hr=°=^>i/'’' 

S-9 = -2<9hr’'”^t9°=V3ht’' 3vr9hr9“=^r>l/'"' 

(6-3) 

are the kernel functions with respect to an integration on a sphere and where 
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^ 1 
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vv * 


5a = <|Shhl^> 


?3 = <lSvhl^^ 


?4 = 


t - Im <S Sf. > 
^5 - vv hh 


E = Re <S S*. > 
^6 vv vh 


5, = I™ 


5= = 


£, = Im 


(6-4) 


the unknown scotterlng parameters. The parameters leading the integral are constants. 
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For each I one must specify receive and transmit antenna polarization states and pat” 
terns such that the resulting system of equations can be solved approximately. There are 
undoubtedly many such specifications. However, there are certain physical considera- 
tions which make the search for the appropriate kernel function {antenna polarizations) 
simpler. 


It has been shown, for example, that in the measurement of a auto”correlation 
coefficient, the kernel function can be approximated by a delta function if the antenna 
beam is sufficiently narrow to resolve the angular behavior of the coefficient (3^, The 
method assumes that the scattering parameter is constant across the significant portion 
of the kernel function. The unknown parameter is withdrawn from the integral and 
the resulting integral expression evaluated. The solution then becomes algebraic. This, 
in effect, is equivalent to assuming that the kernel is a delta function with a weight 
corresponding to the evaluation of the integral expression. The method is feasible 
since the kernel function is sharpened by a product of pattern terms as indicated in 
Equation 6-2, The two-way sharpening effect is illustrated in Figure 6,1 where both 
g and g are plotted. It should be noted that the ordinate scale has been transformed ' 
logarithmically to dB. The kernel function is consequently significant only over a 
very small domain of |(0 ,4)), 0 . 

It would be helpful if the delta function approximation could also be used 
to recover the cross-correlation scattering paramete , An examination of Equations 
6-2 and 6”3 indicates that the two-way sharpening effect is present in the gain functions. 
However, there is no guarantee that B ^ and 3^ will remain constant across the signifi- 
cant domain of the gain functions. Generally the antenna phase factors ore functions 
of (0/J0), On the other hand, if these factors are stationary on the main beam, then 
the delta function approximation can be employed for these parameters also (see 
Equation 6-6). The ability to realize the stationary condition is treated in the sub- 
sequent section. 
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PATTER!!-. AMPLITUDE (DB) 



FIGURE 6.1 INCREASED RESOLUTION OF THE KERNEL FUNCTION 
CAUSED BY THE TWO-WAY SHARPENING OF THE PATTERN 



6o3 Desirable Anirenna Properties 


!f the delta function approximation is to be employed, then ft is desirable to 
have the relative phase and constant across the main beam (See Equation (4~29)). 

This objective is equivalent to requiring that the gain and polarization be stationary 
across the main beam, Chu and Kouyoumjian [39] have derived the conditions under 
which stationary gain and polarizations can be achieved. Coincident stationarity, 
they state, can be realized by ony plonar aperture distribution which is symmetric 
with resoect to two orthogonal axes in the aperture plane. An aperture is planar 
only if the excitation lies in the aperture plane and not orthogonal to it. 

For some center fed paraboloids the above requirement can be met; however not all 
feeds result in a planar distribution even though the symmetry property is observed* 

This is illustrated for the case of a dipole feed. Although the distribution in the 
aperture plane is symmetric, it contains excitation components orthogonal to the plane. 

The orthogortal components are induced by the depolarization property of the para" 
boloid. The far field of such a dish is Illustrated in Figure 6.2. The computation was 
based on a "10 dB taper, a f/D* ratio of 0,36 and a wavelength of 2.16 cm. 

The introduction of variable cross polarized content can clearly destroy the 
stationary polarization requirement. Admittedly the cross~polarized content In the illustrated 
case is small; however, as will be shown later in Section 7,4, retrieval of the cross" 
polarized scattering coefficient can be affected by weak cross polarized pattern levels. 
Furthermore, dipole fed paraboloids with smaller f/D ratios will have a larger cross 
polarized level than illustrated here [44J . 

Recently, corrugated horns [40] (41] and dual mode horns [42] [43] with 
circularly symmetric patterns have been recognized as capable of eliminating cross- 
polarization in center fed paraboloids. The feed pattern of these horns are said to be 
balonced. Mathematically their radiation takes the form 



/cos<|)'\- 

/ sifi (|)' \_ 


( / 


\sm6 / 

\cos0' / 


exp ( -i kp) 


(6-5) 


where the z' axis is directed along the axis of the paraboloid. Chu and Turrin [59] 
have shown that center fed paraboloids with balanced illumination exhibit no cross- 
polarized content in the aperture plane. The far'fields of the above described para- 


* focal length r paraboloid diameter 
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RELATIVE PAHERN INTENSITY (dB) 



ANGLE (OEpREES) 


FIGURE 6.2 FAR FIELDS OF A PARABOLOID WITH 
A DIPOLE FEED 
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boloid wil'h a balanced feed was computed and are shown in Figure 6.3, The cross- 
polarized field was totally absent in the numerical computations. 

Balanced fed paraboloids are suitable candidates for scatterometry work when 
a complete set of scene parameters are desired. Since support struts and aperture 
blockage, in general, introduce cross-polarized radiation it Is important to minimize 
blockage in addition to choosing an appropriate feed. The Cutler type feed with 
balance radiation may be a suitable approach. 

Alternatively, an array of broadly directional radiators is also a suitable can- 
didate,. If the interaction between elements is weak, then the pattern of the array 
is the product of the array factor and the pattern of one of the elements. The polar- 
ization property in the main lobe will be dictated by the polarization property of the 
central segment of the elementary pattern. The polarization will generally be station- 
ary across a small segment of the elementary pattern; and, therefore, the array polar- 
ization will also be stationary there, 

6.4 The Inversion of Scotterometer Measurements 

When a narrow beam scotterometer antenna with a coincident stationary gain 
and polarization property is employed,the scotterometer equation may be approximated 
by* 


WfrW = ^ 2 ^^^ + [13 + 21 ^ cos (e^-Spi<ISj^> + 


i(B, +8.) 


2l4Re[e "'•'<S^Shfi>l + 2RetT5c^®^ + e + 

( 6 - 6 ) 






where 


* The degree of accuracy will be demonstrated in Chopter 7. 
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Angle (Degrees) 


FIGURE 6.3 


FAR FIELDS OF A PARABOLOID WITH 
A BALANCED FEED 
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cos edn 




^ 9 ht%r-^ 9 ,t 9 hr^ cos 9 dfl 


-//' 

= Vv^ 


do 


t^ht ^vr (6-7) 


■ /M 

= //^t V^vrShr 


= //VV^ 9hr = 

= Jj \t >/V^ 


do 


^ ^ A^W^G^G^/(4^z)^ (6-8) 

It has been assumed that observations are conducted over a planar earth so that r = z/cosG . 

It hos also been assumed that the kernel function has sufficient resolution that the scat- 
tering coefficient may be considered constant in the domain where the kernel function is 
significant'. Now suppose that the scatterometer is equipped with a dual linearily polarized 
feed or if necessary two ontenna with orthogonal linear polarizations to assure good isola- 
tion, The amplitude and phase of each feed channel is assumed controllable. Then as 
will be shown below, a series of fifteen intensity measurements with different polariza- 
tion combinations is capable of extracting a complete set of nine scattering parameters, viz,. 
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A pair of measurements is 

required to isolate the real or imaginary part of the complex valued coefficients. The 

transmit" receive polarization states are indicated for each measurement: 


1) vv = g, g^,j = g^^ = 0 ) 




(6“9a) 


2) H-H 


( 9^ = 9„r = 0. g^j = 9|,p ■ 9 ) 

hh' 


Wtp= KI <|S^jR> 


(6-9b) 


3) V-H 


^ ^vt " ^hr " 3’ ^ht " ^vr "" ^ ^ 

«tr= ''I <lSvhi^> 


(6-9c) 


4a) LC-RC ( = 9hr " ” 9* h = 3^ = 90“ ) 
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1 1 ^ 1 
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5a) Elliptical ( ^ ^ 9^^ =”9» ^t “ * ®r “ ^ 

1 11 

V= >'H7<^vI^>" 7<iShh!'>*7’”'<^v<h>] 
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6a) V- Diagonal Linear ( g^^ = g, = 0, g^^ = g 
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6b) V- Diagonal Linear ( g j. = g. gu» = D, g = ~ ' 


Mtr= 


[t <'^v 


7a.) V-RC 


i (6-91)^ 

< 9vt=9" Snt = 9. S,r = 9hr “ 7 s* ^r “ ®'’° > 


'] 


•■•tr' ''I 


[t T <l5^bl ♦ I" < v!I»>] 


1 


(6-9j) 


7b) V-LC 


t 9vt = 9. 9f,t = t), g^^ = g^^ = - g, 3^ = -90 ) 


«tr= 


[7 " 7 <’^hl ^^v4>] 

(6-9k) 


d» 

‘8a) H- Diagonal linear ( = 0, = g, g^^ = g^p = Bp = 0 ) 




1 (6-91) 

8b) H- Diagonal Linear ( = 0, gj^^ = g, g^p = g|^p = — g, == 180°) 
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(6-9m) 



9a) H-RC 


9b) H-LC 



where 


-If 


2 2 
g cos 9dfi 


(6-10) 


The transmit and receive polarizations may be interchanged without affecting the above 
equations o The abcv ; set of equations assumes that the scattering coefficients are de- 
fined with respect t; the antenna frame. As will be shown for narrow beam antenna, 
the above polarization states will retrieve scattering coefficients defined with respect 
to the surface at all but very small incident angles. The above polarization states may 
be incorporated in the equati m v^hich distinguish surface and antenna polarizations 
to develop an inversion technique based on the distinction. These equations are de- 
veloped in the succeerling chapter. 

2 2 

From the above set of equations it is noted that < | > , <| and 

< jS^|^|^> are each derived from a single measurement i.e, measurements (1), (2), and 
(3), respectively. The remaining parameters are isolated by differencing pairs of 
equations. It is clear that if a complete set of scattering parameters is desired, the 
meosurement set is over-specified. If a minima! set of equations is required, then all 
measurement pairs can be reduced to one of the members. It is advisable, however, 
to work with an over-specified set of measurements to reduce the sensitivity to measure- 
ment errors if all the coefficients are desired. It is a distinct advantage to specily 
equati,;.; pairs if a particular complex valued coefficient is to be isolated. The above 
technique does not pre-suppose that the scattered fields are completely polarized. 
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If one further assumes that the scattered fields are completely polarized, then 
the inversion problem reduces to that for non-statistical targets. Various measurement 
schemes have been reviewed for this cose by Huynen [21] . One of these schemes is 
based on field amplitude and phase measurements of a pair of orthogonal returns from 
each of the two orthogonal ilium noting polarizations. Another scheme involves 
omplitude measurements at different polarizations. The latter technique yields a set 
of target invariant parameters which must be transformed to a scattering matrix. In*" 
tensity measurements as described above will, of course, also work. The set of 
measurements may be solved subject to the constraints 


|<s s* >! 

vv hh 

= V 

1 ^ 
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II 

l^> 


1^> 


= V 

|2> 


l'> 


for additional accuracy, Norriinear regression techniques os described in reference 
[45l or \ 46 ] may be employed to solve the system of meosuremenis subject to these 
constraints. 

In retrospect one can also use correlation and cross- correlation techniques to 
isolate some of the non-coherent scattering coefficients. For example to measure 
<Svy is transmitted. During reception both and are cross- correlated 

without and with 90° phase shift injected into one of the channels to isolate the real 
and imaginary parts, respectively. 

Either correlation techniques or intensity techniques as proposed will suffer 
from poor realizations of the desired antenna properties. Since intensity measurements 
are commonly made, this investigation will restrict its attention to the intensity 
technique.' 
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7.0 PRACTICAL CONSIDERATIONS IN RETRIEVING THE SCATTERING COEFFICIENTS 


7.1 Introduction 

In attempting to retrieve the scattering coefficients by the method developed in 
Chapter 6, one Is immediately confronted with the fact that the ideal antenna polat^ 
ization states specified in each measurement are seldom achieved in practice. To de- 
termine the sensitivity of the measurement to deviations from these ideal states, computer 
simulations were conductedo Measurements were simulated on the basis of the complete 
scatterometer equation as developed in Chapter 4 and a scattering characteristic similar 
to that of the sea under low wind conditions. The scattering coefficients were expressed 
with respect to the surface polarizations; and, consequently, all simulated power returns 
involve transforming the pattern information to the surface polarization states to compute 
accurate power returns. Measurements were computed based on known deviations from 
the ideal antenna polarization requirements and were inverted on the basis of the ideal 
antenna specifications. The sensitivity in retrieving each coefficient was thus establish" 
ed, namely, by comparing the actual coefficient with the estimated coefficient. 

The computer simulation was designed not only to determine the sensitivity of 
the measurement to non- ideal antenna polarization states, but was designed to establish 
the beamwidth limitation to realize the delta function approximation for the integrand 
in the scatterometer equation. It was also designed to determine whether the distinction 
between surface polarizations and antenna polarization is important; and if so, under 
what conditions it is important. 

Within the latter portion of this chapter special consideration Is given to the 
sampling requirements when measuring an antenna pattern. The simulations described 
above were based on idealized functional representations for antenna patterns. In 
reality these ideal symmetric representations are seldom achieved (See Figures 6.2 and 
6.3, for examples of non"symmetric patterns). As a consequence, to accurately specify 
the scatterometer integrand recourse to pattern measurements is necessary. The latter 
section of this chapter develops the theory which specifies the density of points at 
which the pattern must be measured to uniquely represent the pattern o This section 
of the chapter Is important in numerically evaluating the inversion parameters In the 
scatterometer equation. 
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7*2 Description oF the Scatl'eromefer Simulation Program 

The reader will recall that the inversion technique developed in Chapter 6 was 
derived without regard to the distinction between surface and antenna polarizations. 

As a consequence to compute the return power accurately from scattering coefficients 
defined with respect to the surface polarizations, the scatterometer simulation program 
was specifically designed to compute the return power on the basis of Equation (4-50) of 
Chopter 4 rather than Equation (4-29), i.e., with the pattern transformation included, 
hor an antenna pattern and a view angle selected externally to the program, the exact 
return power is computed for all fifteen measurements described in Chapter 6, 

The inversion of the resulting measurements is performed in two ways. In the 
first method, called the approximate method, the inversion is performed without regard 
to tl»e distinction between antenna and surface polarizari^ns , It is (erroneously) 
assumed, as in Chapter 6,that the scattering coefficients are expressed in the antenna 
coordinate system. Equation (6-19) served os the inversion model. Since the return 
power was computed on the basis of the difference between surface and antenna polar- 
ization and the inversion was performed without regard to the difference, the distinction 
between surface and antenna polarizations could be evaluated. The second method, 
called the exact method, does not ignore the difference between antenna and surface 
polarizations. The Inversion is based on antenna weights that are computed by trans- 
forming the pattern polarization states to the surface polarization states for each of the 
fifteen measurements. The transformation, in general, "excites" additional scattering 
coefficients above those recognized in tho approximate method (See, for example. 
Equation (4-49),) 

A delta function approximation was also employed in the matrix inversion 
model. The model was based on an approximation of Equation (4-47) and takes the form 



74 
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vv vh 


lyCosQ^dn 


(7-1) 

where fhe Is are defined in Equations (4-50) „ The resulting fifteen equations are 
employed in .a least squares estimation technique to recover the scattering coe Ncients. 
The matrix technique was developed to test whether the fifteen measurements were 
sufficient to invert for the coefficients when the difference in polarizations is recognized. 

In addition to specifying the choice of antenna view angle, the program user 
may, through the use of the input control card, introduce cross pattern amplitude bias 
and relative phase bias into those measurements employing vertically or horizontally 
polarized transmissions or receptions. The return power is accurately computed for all 
fifteen measurements with the biases included. The inversions, both approximate and 
exact me thod*-, are performed, however, without regard to the biases, i.e., they are 
based on ideal antenna states. The sensitivity of the inversions to pattern deviations 
from ideal conditions could thus be studied. 

Cross pattern amplitude and phase biases have precise meanings for vertically 
and horizontally polarized transmissions or receptions. However, for those measure- 
ments requiring simulationeous vertically polarized and horizontally polarized patterns 
(eg,, LC, RC, linear i 45°*), it was more meaningful to conduct Monte Carlo 
studies on amplitude and phase. This technique requires many simulations to be con- 
ducted, Each simulation is based on a different set of deviations in amplitude and/or 
in phase. First and second order error statistics are accumulated from all the experi- 
ments conducted in this fashion. In the measurements requiring simultaneous cross 
potterns, it is evident from Chapter 6 that balanced patterns are required, i.e., 

* Assuming the antenna Is not simply rotated. 
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g^=g^= l/2g* So U was appropriafe to specify the pattern amplitude perturbation in the 
Monte Carlo studies as a deviation from a balanced condition. For each experiment 
the amplitude and phase are randomly perturbed within bounds specified by the user. 

The deviations are based on samples from a uniform distribution so that the perturbed 
gain and phase satisfy 


Sf, ' \ 

K ‘ 1 

$' = $ + Bpg 


(7»2) 


where A < 1/2, 0 and and are random samples from a population distrib- 
uted uniformly over ["1,1] <> °nd 6 are the ideal gain and phase requirements. 

Both approximate and exact inversions are performed for each experiment. The error 
statistics are formed independently for each. The above studies are initiated by 
qaecifylng 2A and 25 on the input control card. 

This program allows the selection of one of four symmetric antenna patterns. 

For any selection It is assumed that both dominant and cross patterns have identical 
functional forms. The relative phase between the patterns (if both exist) was assumed 
stationary. When amplitude error is introduced into an> one of the fifteen measurements, 
the deviation is applied so that the normalized gains safisfy g^(0) -tgj^(0) = 1 on the bore" 
sight axis. The specific pattern options are given by the following functions; 


Pi = 


( sin x/x) 
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Pg = ( Jj{x)/x f 

_ _ / 3 / sin X \2 

Po - ( — ^ C — z cos X ) ) 

X 

Jo(x) 


(7-3) 
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where 


X = kasinS* 
a = aperture radius 

k = 2 tt A 

The above pattern functions correspond to one~woy patterns having respective side lobe 
levels of "13.2, -17o6, -20„6 and -24.6 dB„ In addition to providing a choice in 
pattern functions, the program requires an input parameter denoted os ka to control the 
beamwidtho The beomwidth for the respective patterns are related to ka by the follow- 
ing expressions: 

= 0.88 TT /ka 


AGg = 1.02TT/ka 
AG^ = 1.15 7r/ka 


(7-4) 


AG^ - 1.27TT/ka 

For a fuller understanding of the pattern functions the reader is referred to pages 9.14- 
9o21 of reference [5] , 

The scattering characteristics on which the simulations were conducted are 
illustrated in Figure 7,1 . The coefficients except for the real and imaginary parts of 
< ^hv*^ based on theoretical results reported in reference [10] . The magnitude 

of< ^hh* ^ geometric mean of <[S^^P> and in accord- 
ance with tile results of Chapter 5, The phase characteristic of < was 

assigned to be that for small perturbation theory for a sea water temperature of 293®, 

The characteristics are similar to that of the sea under low wind conditions. In accord- 
ance with small perturbation theory the coefficients < and < are 

extremely small. For the sake of the simulations weak but identical characteristics 
were arbitrarily assigned to the real and imaginary parts of these coefficients. All 
characteristics were assumed isotropic. 

For a complete description of the scatterometer simulation program the reader 
is referred to Appendix D , 
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FIGURE 7.1 SCAHERING CHARACTERISTICS EMPLOYED IN THE SIMULATION 





7*3 Resolution Requirement 
7*3.1 General 

Angular resolution in scatterometry has been achieved either by employing a 
narrow beam or by doppler filtering or by a combination of both. Angular resolution 
is clearly required to search the scattering characteristic* It is also required to realize 
the delta function approximation in the inversion technique. It has been common pracr 
tice- to specify the resolution on the basis of some notion of the scattering character* 
istic* However, when the difference between surface and antenna polarizations is an 
important consideration, a resolution guideline can also be established to assure that 
the antenna polarization coincides with the surface polarization over the significant 
portion of the beam. An expression is developed showing the percent power incident 
on the surface in the orthogonal surface polarization for an antenna whose polarization 
is pure with respect to the antenna frame. The results can be interpreted in terms of ■ 
resolution ^eamwidth;. 

Resolution requirements are also established for the assumed scattering character* 
istics by employing the simulation program. The result expresses the measurement 
accuracy achieved by the delta function approximation with ideal antenna polarization 
specifications, 

7.3,2 Polarization Decomposition of the Incident Beam 

Suppose a scatterometer transmits a horizontally polarized wave E^i when point- 
ed in direction 0 . The total power incident on the surface is given by 
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(7-5) 


When ^ 0 t is decomposed into orthogonal surface components, the above expressions 
can be written by 

(7"6) 

p = -tA— r>r. 1 2 r. . . , 2 , I . 

0 * V* 


2Z 


/|E*'l^ [li. 


VI'] 


dn 


The percent, power appearing In the orthogonal surface polarization is given by 
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(7-7) 
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or 



(7-8) 

The latter expression is simpler to evaluate numerically. 

The above expression was evaluated as a function of view-angle for various 
beamwidthso A Jinc pattern function was employed in the computation. The results of 
the evaluation are shown in the graphs of Figure 7.2, The polarization mis-match as 
anticipated from Chapte*' 5 is greatest at nadir regardless of beamwidth. If is evident 
that small beamwidths ere able to probe closer to nadir without introducing significant 
orthogonally polarized components. The permissable level of orthogonal polarization 
will be treated in a subsequent section. Although the above results were based on a 
horizontally polarized incident wave, a similar result could have been ccxnputed for a 
vertically polarized Incident wave. 

If one chooses to avoid transforming the pattern polarization states to the sur- 
face and accurate measurements of the surface scattering coefficient are desired near 
nadir, then the graphs of Figure 7.2 are helpful in choosing the proper beamwidth. If 
the experiment requires that the cross polarized content be less than, say, -20 dB, 
then the 1% ordinant will specify how close one can probe nadir with various beamwidths. 

An alternative to the above procedure is to employ the exact inversion model 
based on the differences between antenna and surface polarizations. The delta function 
accuracy of this technique for small angles ts developed in the succeeding section, 

7,3,3 An Evaluotion of the Delta Function Approximation 

To determine the beamwidth (resolution) requirement to realize the delta 
function approximation,scatterometer simulations were conducted in the vicinity of 
nadir where angular resolution is required to search the rapidly varying scattering char- 

acferiSTicSa ffi© ability of the delta funci'ion approximation to retrieve each scattering 
coefficient was established at incident angles of 0°, 4° and 8'^, Beamwidths from 1 
degree to 12 degrees were considered. The results are illustrated in the graphs of 
Figures 7,3 and Figures 7,5 through 7.8 for both the approximate and exact methods. 

The performance of the delta function approximation at nadir is shown in 
Figure 7,3 for the approximate method. It is evident that there is little difficulty in 
retrieving >, < and Re > except for beamwidths tn 

excess of 10 degrees. The degradation at large beamwidths is, of course, the result of 
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PERCENT DEPOLARIZATION 
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FIGURE 7.2 DEPOLARIZATION OF THE INCIEDNT BEAM 
AS INDUCED BY THE DIFFERENCE BETWEEN THE ANTENNA 
AND SURFACE POLARIZATIONS 
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FIGURE 7.3 ACCURACY OF THE DELTA FUNCTION 
APPROXIMATION FOR THE APPROXIMATE INVERSION 
MODEL FOR 0 = 0° 
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the antenna beamv# dth interacting with the scattering surface ''beamwidth" . The beam- 
width interaction problem is clearly evident in the error characteristic of Im < • 

Since this coefficient has a notch character at nadir, it is impossible for any non-zero 
beamwidth antenna to retrieve this parameter, 

2 

Unusual error performances are apparent in retrieving <| > , Re < 

and Re < ^ constant 50% error occurs for <| regardless of beam- 

width; whereas a 100% error occurs for the letter two parameters. An explanation for 
the error in P> can be constructed solely on the basis of the difference between 

antenna and surface polarizations, A sirrilar explanation is thought to apply to the 
other two parameters, although no quantitative argument could be constructed. The 
error !n > can be best understood when the antenna and surface polarizations 

are projected on the surface. The surface polarizations will project as a polar grid 
whereas the antenna polarizations will project roughly as a rectangular grid as illustrated 
in Figure 7.4, From these diagrams It is understood that when a vertically polarized 
spherical wave is incident on the surface, half the power appears in the surface vertical 
polarization and the other half in the surface horizontal polarization. As shown in the 
accompanying decomposition diagram both incident components are depolarized by the 
surface and upon their return to the antenna each depolarized component is transformed 
(T) back to the ontenna polarizations. Upon transforming back to the antenna polar- 
izations, one half of each depolarized component is transformed into the antenna 
horizontally polarized state. As a result, the inversion based on the antenna polarizations 
is 50% low. This result indicates that it is futile to recover os defined with 

respect to surface polarizations with a recovery technique based on the antenna polai- 
izotions , 

It is also informative to examine the power structure for a cress-poiarized 
measurement. The third rov/ of Table 7,1 shows how the return power is distributed 
among the scattering coefficients for a cross-polarized measurement. A sizeable 
contribution arises from the polarized coefficients (columns 1 and 2); however, the sum 
of those components is fortunately cancelled by the contribution from Re 
(column 4)^ The cancellation is assured by the Isotropic character assumed for the 
surface. 

For the same incident angle no results can be reported for the exact method. 

For the nadir angle the observation matrix is singular. The singularity Is plainly evident 
in the observation matrix as shown in Table 7,2. The reader will observe that the 
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FIGURE 7.4 COMPARISION OF ANTENNA AND SURFACE POLARIZATIONS 
AT NADIR WITH DECOFiPOSITION DIAGRAM TO EXPLAIN 
CROSS POLARIZED MEASUREMENTS 
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TABLE 7.1 POWER COMPOSITION MATRIX FOR A NADIR MEASUREMENT 
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7.2 OBSERVATION MATRIX BASED 

ON SURFACE POLARIZATIONS 


following pairs of observations are identical: 1 and 2, 3 and 5, 6 and 7, 8 and 12, 9 
and 13, 10 and 14 ond 11 and 15, The rank of the matrix is consequently 8, For iso” 
tropic scenes the singularity may be removed by solving the system of measurements 
subject to the constraints: <| S^^F>=<! S^J^> , ^ 


< = y^| Syy(^ > <[ S|^j^|^>, As a result of the constraint there are only 

five independent parameters. This result was anticipated from the ^‘Gedanken Experi" 
mente" referenced in Chapter 3, 

In general, retrieval of the scattering coefficients at nadir is a difficult task, it 
is merely coincidence that the approximate method yielded as many accurate estimates 
as it did. If a scene Is anisotropic or if it has a peculiar character where t‘"^vv^hh*^/ 

^ <1 S J > <( S|^j^ p >, then there is no assurance that either method will work. 

See the fourth column of Table 7,1 where it is evident that Re"^ ^vv^hh*^ plays an 
important role in forming the polarized measurements. Careful investigations at nadir 
will require very narrow beams to search nadir asymptotically if the coefficients are to 
be reported with respect to the surface polarizations. 

The accuracies of the delta function approximation for the approximate and 
exact models at a view angle of 4^ are shown. Figures 7.5 and 7,6, respectively. From 
Figure 7.6 it is apparent that the approximate method can be employed with reasonable 
accuracy (0,5 dB) to retrieve all coefficients if the beamwidth is less than 3°. A beam- 
width as much as 10° can be tolerated for the recovery is restricted to the polarized cor" 
efficients. The exact inversion method will permit beamwidths up to 9° in retrieving 
all the scattering coefficients. Similar results are apparent in the error characteristics 
for a view angle of 8° (See Figures 7.7 and 7.8), 


7,4 Antenna Requirements for the Accurate Recovery of the Scattering Coefficients 
7,4.1 General 

A number of simulations were conducted at various incident angles with and 
without biases and also with and without random perturbations introduced into the 
measurement. These simulations served as c training set to identify the particular 
scattering coefficient or coefficients which primarily contributed to the error character- 
istic for each scottering coefficient. Invariably the best single parameter to which the 
error could be attributed was the magnitude of < > . The magnitude of this 

parameter corjveys a notion of the size of < P > and ^ >as well , These 

three coefficients generally interacted to introduce an error in the measurement when 
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the antenna transmission and reception properties deviated from the ideal state specified 
in Chapter 6, 

The training set also made it apparent that the error characteristics were primarily 
governed by the level of the cross polarized leakage for those measurements involving 
linearlly polarized hansmission or reception* The level of the leakage is expressed in 
terms of one-way depression relative to the dominant pattern. The relative phase 
between the dominant and leakage patterns was treated as an independent error parameter* 
A bias error study was, therefore, applied to the retrieval of < 

and • Although the latter two coefficients involve balanced 

cross patterns during reception, studies showed that the error performance was largely 
insensitive to smoll deviations from a balanced condition^ small deviations from the 
required phase condition in comparison to leakage appearing in the linearlly polarized 
transmission. Now in the case of > , It is more meaningful to consider Monte 

Carlo studies since both transmissions and receptions involve balanced cross patterns, 

Al! simulations were conducted for a one degree beam having a (2 (x)/x)''^ 

pattern. The resulting error characteristics apply equally as well to approximate or 
exact inversion methods. When translating the performance to small incident angles 
where the antenna and surface polarizations differ significantly across the beam, then 
one must assume that the Inversions hod been performed by the exact method , The 
graphs of Figure 7,2 serve as a guide as to when the matrix method must be used. 
Simulations with the other pattern functions yielded similar results and so are not reported. 


7.4,2 Error Characteristics 

The error characteristics for the recovery of are shown in Figures 

7,9 and 7,10. The results are shown for two phase conditions in which 0^, “ 0® 
and 3^ 90°, respectively. These two conditions result in extremal error character" 

istics in which the maximum error results from one phose condition and a minimum error 
from the other condition. The extremes ore induced by a sign change in the 
contribution from Re< ° dominant parameter. As shown by Figure 7.9 


VV III! 2 

there is no difficulty in retrieving the dominant scattering coefficient > except 

for a cross polarized pattern less than 10 dB beneath the vertical polarized pattern and 
a large separation between <1 S^J and <1 Sj^jJ ^ , The weakened dominant pattern 
results in less return power from < 1 S^J ^ , A further reduction occurs when and 

< are significantly weaker than<] , A similar result occurs when 

3 ^ =B 90° {Figure 7,10). The error is slightly larger because the coefficient 
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Re causes a "negative" power contribution, resulting in on even smaller 

resultant power. Recall that responds to the product pattern ^/"g^g^ and 

its sign is controlled by the sum /3. + jS (See Equation 4'“29). 

* ^ 2 

The retrieval performance, when attempting measurements of < l%hi 

illustrated in Figures 7,1 1 and 7.12 for pattern phase conditions corresponding to |S^ = 

0° and jSj, = = 90°, respec^'iveiy. Since > is generally less than or 

equal to< |S |^>, one can anticipate a poorer error characteristic. For the case 

where = 3^.= 0 , the ability to recover < > is shown to be strongly dependent 

>[ . Positive power contributions are made by both 

The resultant power in this case is excessive. When ( 3 ^ ~ 


on its separaJ -jn from [<S /v^hh* 

i2 


< S. 


vv' 


■> and 


0= 90^, the contributions by <5 vv^hh*> is negative and partially cancels the 

contribution. As a consequence, one may suspect that the latter phase 

condition yields a slightly better error characteristic. Comparison of Figures 7,11 and 

7,12 demonstrates that this is the case. From either graph it is observed that when 

2 

< > is 10 dB lower than , the antenna cross polarization level must 

be less than "30 dB for a error less than 0,5 dB. When the separation is 5 dB, the 

antenna cross polarization level must be better than "20 dB. The latter is probably re" 

presentative of the sea for angles of incidence up to 70°, 

The error characteristics for retrieving <|S^j^|^> for the same two relative phase 

conditions are shown in Figures 7.13 and 7.14. From Figure 7,13 it is apparent that 

the weakness of the scattering coefficient in the presence of cross"leakage makes it 

2 

very difficult to isolate. The ability to measure > is shown to depend strongly 

on its depression from the polarized scattering coefficients as conveyed parametrically 
by l<s,,s . Figure 7,13 represents a worst cose situation in which all the 
dominant coefficients to include Re < make positive contributions to the 

return power. This situation is consequently useful for formulating a criteria for accurate 
measurement of > « It has been common practice to judge the ability of an 

ontenna to measure cross"polorized coefficients by its one "way and in some cases by 
its two-way isolation in comparison to the separation between the polarized and cross" 
polarized coefficients. The graphs of Figure 7.13 sh^%w explicit the antenna require" 
ment. If a 0,5 dB accuracy is desired and if lies X dB beneath l"^^vv^hh*^l ' 

then approximately X + 16 dB one-way isolation is required. The above result indicates 
thot one must not only consider the polarized coefficients in making a judgement on 
an antenna but o complex valued coefficient must also be considered. 
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It is apparent from the graphs of Figure 7,14 that if the phase of the leakage 
pattern can be adjusted to 90° during transmission and reception, the level of 
leakage is almost immateriaL In this case contributions by and 

are almost entirely cancelled by the contribution from Re • ^hese results 

show that if the phase of the cross leakage can be adjusted for 3^ = = 90 , the 

stringent requirements on the cross pattern amplitude can be relaxed. 

The error characteristics for Re < ^yv^hh*^ "^^vv^hh* ^ shown in 

Figure 7.15c Monte Carlo studies were performed to construct this characteristic. The 

random deviations in amplitude {from balance) and in phase were uniformly distributed. 

Maximum deviations are indicated on the graphs. It is apparent that the real part of 

< > is easy to recover. Phase perturbations have little effect on the accuracy. 

The recovery of the imaginary part appears to be more difficult; but this is mainly a 

2 2 

result of its weak response in comparison to <| > and <1 > . 

The error characteristics for the cross-riorrelation coefficients < S .5. *> and 

vv nv 

are all shown in the graphs of Figure 7,16. Both extremal phase con- 
ditions are superimposed on the same plot. The graphs show that the real parts of the 
coefficients are difficult to retrieve if = = 0, Similarly, the imaginary parts are 

difficult to retrieve if ^ = 3^ = 90°. On the otherhand, the imaginary part and the 
real parts are easily recovered if 3^ = 3^- 0° and 3^ = 3^= 90°, respectively. The 
graphs also indicate that the ability to retrieve the coefficients is dependent upon the 
separation of the coefficients from the real or imaginary part of >. it is 

again evident if the correct phase property is employed that a reasonable accuracy can 
be anticipoted, 

7.4,3 Alternatives 

When the idea! antenna states as specified in Chi er 6 cannot be approximated 
reasonably and if as a consequence significant error is introduced into the measurements, 
the experimenter has recourse to specifying the complete antenna polarization states he 
is able to achieve. As long as he approaches the desired states and performs an adequate 
number of measurements, he can be reasonably assured that inversions based on the com- 
plete scatterometer equation will yield improvements in the estimates of the coefficients. 
The inversion model should be tested to determine whether his system of mecsurements 
is well conditioned. At least nine measurements must be performed unless one has prior 
knowledge that some of the coefficients are negligible. In this technique one must re- 
concile with making at least nine measurements; whereas if the beamwidth constraint 
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is met, a single scattering parameter can be recovered with at most two measurements 
if the antenna specification can be realized. 


7.5 Evaluation of the Inversion Parameters 


7o5.1 Introduction 

Essential to accurate inversion of scatterometric measurements is the knowledge 
of the actual antenna pattern. To form the integral weights for each scattering co“ 
efficient, the pattern and phase functions must be numerically integrated over the main 
beam and perhaps the first side lobes. Since pattern information is seldom available in 
functional forni, one is dependent on measurements. In measuring the pattern, the question 
arises as to what sampling density is required to adequately specify the pattern. The 
sample requirement Is derived on the basis of simple aperture theory. The results of 
the theory are applied to the SKYLAB S“193 antenna J’o illustrate the sampling require- 
ment. 


7.5,2 Derivation of the Pattern Spectrum 

It is well known that the far field E of an aperture type antenna is related to 
the aperture illumination function, A(x,y), through an inverse Fourier transform 
relationship [19] 


00 


where 


) exp [a( kj^x + )] dx dy 

(7-9) 


” (iA r) exp (- jkr) 

k = ksinO cosj^ 

X 

k. = ksinS sin^ 

k = 2^A 


(7-10) 


The relationship is considered valid for spherical polar angles 6 satisfying cos 9^0.9. 

A is assumed to be a real function* so that the main beam of the antenna is located about 
the positive z axis. Now it is convenient to rewrite the above expression in the form 

. < V I ■ «,//.< ) exp j 2 tt ( -f f^n ) dp 

{7-11) 


* A uniform phase distribution across the aperture. 
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where 


r f = sin 9 cos 0 
s . 

f = sinO slnj3 

n 


(7-12) 


c = xA 
n =yA 

K|= (] Vr) exp (-jkr) 

The far field power paftern P Is given by 

P ( ) = KgEE* (7-13) 

where K 2 Is a suitable constant. The Fourier spectrum cf P is given by 

K2;?[ee*1 (7-, 4) 


or 


= A*A 


(7-15) 


where ^ Is the autocorrelation operator. Specifically 


.-?[P] " Kg JJ*A( ) A( C+a, n+3 ) d? dn (7-16) 

and implies that the spectrum of the power pattern is proportional to the autocorrelation of 
the aperture distribution and Is therefore band limited for finite apertures. 

For circularly symmetric aperture distribution a similar theory could have been 
derived if the initial expression had been transformed to the Bessel-Fourler integral. 
However, seldom are aperture distribution circularly symmetric, as a consequence, we 
have a more general result. 

7.5,3 Sampling Requirement 

Now suppose that an aperture has maximum length x^ and maximum height y^. 
From the above result and the illustration in Figure 7,17, it is clear that the spectrum 
of P is restricted to the product domain (xy"X, X **7(/^* ^7^^® sampling 

theorem, the pattern can be specified uniquely if samples are taken at 


( ) = ( 




nx 


2Xo * 


2yo 


(7-17) 
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where m, n e .... ... . The above result can be written in terms of 

6 and 0 by means of (Equation 7’“14). Specifically 


and 


sin e cos <i> - 


sin 6 sin 4> - 


fT]\ 


2x, 


n\ 

2yo 


(7-18) 


(7-19) 


As can be easily shown the above relationships require that the antenna pattern be 
sampled at points 0^^ satisfying 

.2 2 


»nm = 


. -1 r 

in 


( 


m 


X, 


V )^] 


(7-20) 


'^'mn = 


-1 _n_ 
m 


(7-21) 


In the principal planes the above sampling requirements reduce to 




m\ 

2Xo 


(7“22) 


in the “x" plane and 




(7-23) 


in the "y‘* plane. Between the planes in the pattern must be sampled in accord with 
Equations (7-22) and (7-23). 


7,5,4 Illustration 

To develop an understanding of the sample requirement. Equations (7-2Q and 
(7-21) were evaluated for an aperture having a maximum dimension of 1 ,12 meters in 
the X as well as the y dimension and Illuminated at 13.9 GHz, The sampling points 
for one quadrant out to approximately seven degrees in theta is illustrated in Figure 
7,18, Sampling in the remaining quadrants is performed in an identical fashion. 

It is noted that the sampling array forms a square matrix in polar coordinates where 
theta is represented as the polar radius and phi as the polar angle. 
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FIGURE 7.18 SAMPLING POINTS FOR A SQUARE APERTURE 
1.1 METERS BY 1.1 METERS OPERATING AT A WAVELENGTH 
OF 2 .16 CENTIMETERS 
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The above results are representative of the sampling requirement for the S"193 
antenna [12]. Although, since the physicci aperture was underbill ummated, the 
above result is a very conservative sampling density. A description of a program 
which computes the sampling points when given the aperture dimensions appears in 
Appendix E. 


8.0 CONCLUSIONS AND RECOMMENDATIONS 


8.1 General 

The scatterometer equation was derived for scenes whose mean plane is flat 
and for an antenna having an arbitrary polarization. Ten scattering coefficients were 
identified for scenes not satisfying reciprocity and six were identified for scenes 
satisfying reciprocity. Some of the scattering coefficients were demonstrated to be 
complex valued and were shown to impart a relative phase between the vertically 
and horizontally polarized scattered components. As a result of the complex valued 
coefficients, the definition of a scattering coefficient had to be generalized. A new 
descriptive notation for the coefficients was suggested. 

As a consequence of linearily scanning the scene to obtain a spatial average, 
it was demonstrated that the scattering coefficients must satisfy Schwartz’ inequality 

l<S..S*^>r< <|S..r><|5^,r> (8-1) 

where i, j, k or I = v or h. This naturally implies that equality is assured for the 
polarized and cross-polarized coefficients. However, equality should, in general, 
not be anticipated for the cross- correlation scattering coefficients. It is this feature 
which distinguishes coherent and non-coherent scattering coefficients. As a result of 
this inequality, scatterometer returns can be partially polarized. Furthermore, the 
inequality also implies that one cannot employ the properties of the (coherent) scattering 
matrix to describe non-coherent measurements. For a coherent target five independent 
parameters (from the scattering matrix) are required to describe its scattering coefficients. 
However, a non-coherent scene requires as many as nine independent parameters. 

The scatterometer equation under the reciprocity assumption was extended to 
account for the difference between antenna and surface polarizations. It was 
illustrated that the difference in polarizations was significant only ck nail view angles 
for narrow beam radars. The effect of misalignment can be minimized by reducing the 
beamwidth os one approaches nadir as illustrated in Figure (7.2), Minimizing the mis- 
alignment is important if an experimenter wishes to compare his measurements with 
theoretical predictions which are invariably reported with respect to the surface 
polarizations. It is shown, for example, that a cross polarized measurement at nadir 
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cannot be interpreted as an attempt to retrieve defined with respect to 

the surface polarizations. In view of the difficulty in interpreting measurements at 
small angles with respect to the surface polarizations, it is recommended that the 
nadir region be probed in an asymptotic sense with a very narrow beam antenna when 
scattering parameters are to be reported with respect to the surface polarizations. 

When a scene has an anisotropic behavior at small incident angles, it is particularly 
advantageous to report parameters in this fashion. 

A measurement and inversion technique was proposed to measure all nine 
scattering parameters. The technique was formulated without regard to the distinction 
between antenna and surface polarizations. Since the difference between the polar* 
izations is negligible for narrow beam radars at Cil but the small view angles, the 
formulation without alteration is valid there. In addition, it was shown that the 
system of measurements (antenna polarization states) is sufficient to retrieve all the 
parameters qt small incident angles under an isotropic surface assumption if the inversion 
is based on the extended formulation, i.e,, accounting for the difference between 
antenna and surface polarizations. 

. The computer simulations based on the above technique demonstrated that the 
dominant scattering parameters could be recovered with modest realizations of the antenna 
polarization requirements. However, retrieval of the weaker scattering parameters, as 
diown by the simulations, requires more careful adherence to the antenna polarization 
requirements. Cross polarized leakage in the case of linearily polarized transmissions 
or receptions causes the antenna to couple to the dominant scattering parameters. The 
leakage results not only in coupling to the real valued coefficients but also to the 
complex valued coefficients. The degree of coupling depends strongly on the relative 
phase of the leakage as well as on its amplitude. For scattering characteristics similar 
to that of the sea {where < ^yy^hy*^ "^^vh^hh*^ considered weak), strong 
undesirable contributions can be anticipated from < I , and 

, as demonstrated by the simulations. For a scene having randomly oriented 
linear re~radiators such as vegetation one can anticipate not only strong contributions 
from the above coefficients but also from 

<Sj^v^hh*^* scatter ng coefficients have been cited to emphasize that the 

.scattering processes are differe it although under the reciprocity assumption there ore 

■ 

only two independent coefficients. 

It is evident from the simulations that when only the amplitude of the ortho** 
gonot leakage is known and not its phase, stringent specifications on the amplitude 
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are required fo achieve, say, an accuracy of 0*5 dB for the weaker coefficients. This 
is illustrated when< [ S^|^| > is to be recovered from scenes having weak cross- correlation 
coefficients and < . When it is suspected that< j is X 

dB beneath the geometric mean of < [ 5^^ \ and < the permissible level 

in the cross leakage is -(X + 16) dB. On the otherhand if the phase of the leakage can 
be adjusted so that It is at or near 90° (or it is known to be near 90°) during transmission 
and reception, then the amplitude specification can be relaxed os demonstrated by 
Figure 7.14. For scenes in which Re and Re are dominant, 

this same phase condition can minimize contributions by these terms in the case of 
cross-poiarized measurements. This may be concluded by an examination of Equation 
(4»29). 

Although the assumed scattering characteristics reflected a wide latitude of 

conditions, the error characteristics generated here are by no means exhaustive. The 

retrieval accuracy to some degree Is dependent on the assumed scattering characteristics. 

For example, S, . * > and <S were assumed weak and their error character- 

> ' vv vv vn hn 

istics reflected this weakness. It is recommended that simulations similar to those report- 
ed in this effort be conducted whenever significantly different scattering behaviors are 
encountered. The antenna specifications can be established on the basis of these 
simulations. One may employ the program described in Appendix D in which case sub- 
routine SIGMA must be replaced with a subroutine that will compute the scattering 
parameters of interest. 

It is evident from these efforts that there is a fine opportunity to extend the 
three standard measurements to nine measurements. When the distinction between 
polarizations is not important,* any combinations of measurements can be selected to 
isolate particular coefficients. It is intriguing to consider certain combinations of 
measurements to observe soil moisture, crop maturity, etc. From small perturbotion 
theory there is evidence that the cross- correlation coefficients may contain additional 
information on the dielectric property of the scene when compared with the auto- 
correlation coefficients. The comparison of like and cross- correlated coefficients may 
be the key to distinguishing dielectric effects, say in agrarian scenes, from volume 
roughness effects. 


* This can always be acti:?ved at all view angles except nadir if the beam is 
sufficiently small , 

no 


The above results also have an impact on emission theory and radiometer measure" 
ments. It is clear tnat the backscatter coefficients employed within this effort can be 
extended to the bi"static case. As a consequence, we may address emission theory 
from the aspects of bi-static coefficients as Peake [241 did. Generalizing Kirchoff's 
radiation law, Peake has shown that the definition of emissivity, when assigned 
standard surface polarizations, may be related to the bistatic differential scattering 
coefficients in the following way 


e 


P 


1 - 



p^q 


(8-2) 


where the integration is performed over all incident angles. The corresponding bright- 
ness temperatures were given as 

T^ - T^ (8-3) 

P P s 

where T is the physical temperature of the emitting surface. Peake's formulation for 
s 

brightness ignores the possibility of correlation between emitted components. When 
correlation exists between the components, the concept of brightness temperature must 
be extended as shown by Ko [ 47 \ . Ko had shown that an emission of total brightness 
(intensity) and with a normalized coherency matrix p can be , regarded as a unique 
superposition of two coherent oppositely polarized emissions, i.e.. 


where 


Pll 

P12 


Pll 

P12 

" 82 

Pll 

P12 

P21 

— t 

P 22 

1 

P21 

P22 

P21 

P22 


(8-4) 


Pll - P22 

Pl2 = -P2i 
P21 = "Pl2 
P22 - Pll 
Pll + P22 = 1 

+ P^^ = 1 


(8-5) 


III 


Temperatures are assigned to the brightness according to the Rayleigh" Jeans law 

= B. X^/k (8-6) 

Arbitrary measurement of this emitted field, say, with any two orthogonal polarizations 

will not necessarily result in any unique temperatures. The correlation between emitted 

components plays an important role in defining the brightness temperatures „ Within the 

context of bi-static theory, the cross-correlafion coefficie.'its <S <S S , 

f* vv hh ' vv vh 

^^hv^hh* ^vv^hv*^ and < S^h ^ establish this correlation. For some surfaces, 

the first three coefficients are not important unless the emissions within the radiating 

body are correlated. Under this circumstance the correlation is governed by < S^Sj^y*> 

and f I*®*/ by the correlations which the emitting surface induces. These 

cross- correlations for the sea are assumed negligibly small . The brightness temperatures 

are, therefore, given by the vertically and horizontally polarized emissions and the 

corresponding decomposition into coherency matrices is given by 

Pll Pl2 0 

= kT a' * kT-Zx"" 

Pai P22 0 0 

For an agrarian scene the above simple decomposition may not occur at all, since 
multiple reflections are likely to induce correlations info the emissions internal to the 
radiating boundary and because nr\v*^ and < •^ol' negligible. 

Therefore the brightness temperature concept must be altered for agrarian scenes. 

8,2 Final Remarks 

The above observations as well as the developments in the earlier chapters 
indicate the importance of having derived the complete scofts^-cmeter equation and in 
particular having derived it in the framework of coherency theory. The interaction of 
the transmitted fields with the scattering surface was expressed as a transformation of a 
coherence matrix (Equation (4-41)). Interpretation of the scattered fields from its 
coherence matrix imported meaning to the cross-correlation scattering coefficients. 

The complete scattering action of the surface, when interpreted in the context of 
emission theory also enables one to interpret the coherency properties of microwave 
emissions. The reception of scattered or emitted fields is also expressed as the product 
of two coherence matrices. 
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Alfhoudh a fuller interprefaHonal basis lies in coherency theory, practical 
application of the theory has led to a technique for measuring all six scattering 
coefficients. The measuring technique was evaluated for practicol antennas. As a 
result of this evaluation it becomes apparent that measurement standards or standardized 
reporting procedures or both should be instituted. There is also a clear need to dis- 
tinguish scattering parameters reported with respect to the antenna polarizations from 
those reported with respect to the surface polarizations. The measurement of weak 
scattering coefficients requires stringent realization of the antenna polarization require- 
ments, The documentation of the antenna transmission and reception property must be 
complete, to include amplitude and phase properties, to valldote a measurement. Until 
such a procedure is followed it could be erroneous to report, for example, "cross- 
polorized" measurements as cross-polarized scattering coefficients. To assist the ex- 
perimenter it is also clear that a program should be initiated to develop a scatterometer 
antenna whigh ts capable of meeting the antenna specifications for most if not all the 
scattering parameters for a variety of scenes. 
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APPENDIX A 

CorrelaHon ond Cross~Correla1-ion 
Products From KIrchhoff Theory 


1.0 INTRODUCTION 

The scattering and coherency properties of a finitely conducting random 
surface satisfying the Kirchh off approximation are investigated within this appendix. 
The expressions for the polarized scattered fields in the plane of incidence ore 
specifically derived for both vertically polarized and horizontally polarized incident 
plane waves. The scattered fields are derived under the assumption that the surface 
slopes are small . {Only zero order and first order slope terms are retained within 
the derivation). The resulting expression is specialized to the backscatter case to 
derive the self~correlation and cross"correlation scattering coefficients. The angular 
coherency of the scattered fields about the backscatter direction is also considered. 


2.0 THEORY 


2.1 General 

For a plane wave incident with direction n. on a gently undulating finitely 
conducting bounded surface Fung [36] has shown that the far field scattered in direction 
rig is given by 


E 


s 



_ _ 1 jk(n^-n.)«p 

' (hxH ) j e ' ^ ^ dS 


(A-1) 


, -jkR 

where -jke 

K = 

4-FrR 


(A-2) 


R = distance from the surface to the far field point 

^ " position vector from an origin local to the surface to a point on that surface 
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■nxT 

= [ (l+R„)Ca-T. )(nxT.)- 

(1-Ry) (n* n . ) (a-d . )t^.j I'EqI (A-3) 

rTxH 

= [-(1+R^) ) (TTxT. )- 

(1 - Rj^ } ( a'.T. }tj \ Eq\ (A-4) 

E 

_ j(£Ut- k-*p) 
= E a e 

(A-5) 

kj 

= kn"^. 


k 

= 2tt/X 



ff. xrT 


s- 

iTi.xnl 



= TT^xt. 



7j = inirinsic; impedance 

R = Fresnel coefficient for horizontal (vertical) polarization 
v,h 

The geometry associated with the scattering problem is illustrated in Figure A-l, 
[ti.> ft", forms an orthogonal triad of vectors at each point of the random surface. 
The plane of incidence coincides with the y - z plane. Now expanding the terms 

Figure A-1 , Scattering Geometry 
z 



within the integrand we have 

n^x(nxl) = |Eq1 j (1+R|^ ) (a*T. ) - 

{1-Ry ) (TT •’n^, ) (T •? . ) (n^xt . ) | (A-6) 
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and 

-'j?ngX{?gX(nxH')) = IEqI j (1+R^) (a-3‘. ) [(n^- t.)ngXn - {ns*")ngXt^] - 

(n-F.)a-Rh)(i't.)t.| (A-7) 

where a radial component was dropped under the far"fie!d assumption. Now restrict 
observations to the plane of incidence so that 

rTg = sinGg sin<j5g Ty + cose^ (A-S) 

with 0^ = Tf/2 or ~^/2, depending whether the forward or back scatter quadrant is, re- 
spectively, chosen. Denote 

n*. = sine, iy - cose^.T^ (A -9) 

where 0j is the incident angle. Now it is easily shown that [36] 

t. = [(sine.-cose.ZyH^ + cose.Zj^Ty + sin9,Zj^T^]/Dj (A-10) 


and 


where 


[^x^x cos 0^.(cose.Zy-sin e.)iy + sin a. (cos e.Zy-sinB^. )i J/n^ 

(A-n) 


of = + (sine^-cosa.Zy)' 

Z ? Z(x,y) 


(A-12) 


Z - 
x,y 3x,y 


Then 


n • t j = • (sin8 sin<|j cos e. + cose sine. )Z /D. 
n*n^ = “(sinS.Zy + cose^)/D2 

iig-n = (cos0g - sin0gSin0gZy)/D2 (A-13) 

" [(sin0g3ln<J)gSine.-cos0gCos0.)Zj^Tj^ + cos e^Csin 0.-cos e-Zy)Ty 
-3 1 n 9gS i n ( s 1 n 0. - cos 9^ Zy )T^ j/D 

n^xn = [(sina^sineg + cas6gZy)ij^ - cosa^Z^Ty + 
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v/here 


4 = «y 


(A -14) 


2,2 Horizontally Polarized Incident Wave 

Suppose a = 1 , i.e., the incident wave is horizontally polarized. Then 


and 


a-t. 


a * d. 


-(sine^-cose^Z^)/ 




(A-I5) 


(A-16) 


The horizontally polarized field scattered in the plane of incidence may be shown to be 
given by - 


s (p 


= K 


where 


// 


I^l^exp [ik(n^-n^. )•/>] dxdy 


(A-17) 


(A-18) 


^hh 


= r-(l+Rj^)(cos 0 g-sinegSin(j>gZ ) + (l-Rj^)(cose.+sin0.Zy)] |Ej 

■ (A-19) 

Only terms to first order in Z and Z have been retained. To the same order it may be 

X y 

shown that the depolarized component is zero, i.e.. 




(A-20) 


2.3 Vertically Polarized Incident Wave 

Suppose a = -cos6|T ^ ” sin 6.1 
on the X - y plane. Then 


i.e,, a vertically polarized wave is incident 


. a-t, = -Z/D, 

and 

a*d. = {sin9^.-cos0.Z )/D, 

12 ! 


(A-21) 

(A“22) 




The vertically. polarized scattered field may be shown to be given by 

^s* ’0s 

where 

ieg = cos0gSin4)gTy - sine^T^ 


(A“23) 


(A-24) 


" sin<!)g[-(l+Ry)(cos9^-sinegSin(f)gZy} + (1-R^)(cos0. + sin9.Z )] lE^I 

, (A-25). 

Only terms to first order in Z and Z have been retained in I . To the same order it 

y w 

may be shown that the depolarized component is zero, i<.e,. 


f • t , = 0 

S (|) 


(A“26) 


2,4 Linear Approximations for the Reflection Coefficients 

It is necessary to understand that the reflection coefficients are functions of the 

local incident angle and are therefore functions of the local slopes, Z and Z . For 

X y 


small slopes we may approximate R^and linearily by 


az, » azy ^ 


(A-27) 


Now the relation between the local incident angle 0’ and the local slopes is 


cose = -nj*n 


or 


= [ZySine. + cos6.]/D2 


(A-28) 

(A-29) 


The derivatives within the linearily approximated reflection coefficients can then by 
the chain rule be written 


0R, 


v,h 


a.%., a»’ 


az. 




be’ 


dz. 


(A-30) 


x.y 


An evaluation of the derivatives yields 


as' 


Z =Z =Q 
X y 


V«j,-sin^0^. 
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he’ 




hi. 


2 =Z =0 
X y 




Z =Z =0 
X y 


-26pSine^R^(9^) 

ye^-sin^0.( e^cos^e.-sin^e^ ) 


(A-31) 


= 0 


X -1 


It hds been recognized in the above expressions that 0' = 0. when Z - 2 - 0. We 

I X y 

finally have the approximate expressions 


Rv(Zjj,Z^) = R^{0.) + 


2e^s1ne.R^(9.)Zy 


2 2 2 
e^-sin e^(e^cos 0^.-sin 9.) 


(A-32) 


and 




2stn0^Rj,(0^-) Z^ 


(A-33) 


2,5 Partial Evaluation of the Field Integrals 

The evaluation of the polarized field expressions requires that integrals of the type 

“ ^^Zyexp[jk(Wg-n.)-P]dxdy (A“34) 

be considered. By employing an integration by parts technique, specifically by letting 


and 


we get 


dv = exp[jk(coseg+cos0^)Z]Zydy 
u = exp[jk(sin6gSin<i>g-sin0^)y] 


(A-35) 

(A-36) 


/■ 


Intg = / vu( dx 

30undary 


// 


^ e>«p[jk(ng-n^)-pl 


(A-37) 
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where 


sin0,sin<j),- sine. 

S O & 1 

C COSe+COS 0 . 

5 1 


(A-38) 


The first term is identified as the edge effect and may be neglected. 

When the above results are incorporated In the field integrals we can write a 
unifying expression 




(A-39) 


where 


I =* 


^exp[jk{ng-n^ 


)-p] dxdy 




+ (l-Rj)sine^.3-~ 


■ tl-Rj+-|--^)cos0. - [(l+Rj)sin6gSin 




s s 
(A-40) 


H.(j=v or h) denotes the like polarized field component when a jth polarized wave of 
amplitude illuminates the surface. The reflection coefficient and its derivative are 
evaluated at the incident angle. 


3,0 THE BACKSCATTER COEFFICIENTS 

Now specialize to the backscatter case. Specifically, let 0 = 0. and 0 = - ^2,, 

» I 3 

Then 

^ (A-41) 

where 

(A-42) 

The differential scattering coefficient employed in this effort is given by 

* 2 
<E.E,>R*^ 

; — (A-43) 

|E^| Acos0^ 
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where A is the illuminated area. Consider the ratio of the intensities 

<E-E?> 9 

<II*> 

. . 2 a K 


(A-44) 


v/here 


JJJJ exoH2ky(y^-y2)]<exp[j2k2{Zj-Z2)]>dXjdx2dx2dy2 


V = ^sine. (A -45) 

• h “ {A-46) 

2 

Suppose that 2-j and Z 2 are joint gaussian variables with zero mean, variance ^ and 
correlation A(x^, X 2 , y-j, • Then it is easily shown from the characteristic function 
properties of gaussian variables that 


<exp[j2k,(Zj-Z2)]> ■ exp[-4k ^cr^(l-A)] 


(A-47) 


Now transform the resulting integral to the center of mass coordinates. Let 

U = Xj^-Xg 

V = y^-yg 
Xg =* Xg 

yg- V2 

The integral then can be written as 


<U*> = 


Further transform the integral to cylindrical coordinates where 

u = pcos^ 

V = psin^ 

Xg “ P'cos^ 
yg = psin £ 


(A-48) 


JJJJ exp[-j2kyV] exp [-4k^^ a^(l-A)] dudvdXgdyg (A“49) 


(A-50) 
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If it is further assumed that the surface is statistically stationary and isotropic, then 


<ri*> = 


JJJJ exp [-jak^sin^pJexD o-^{l-A(p)l P dpd^p'dpa^ 


(A-51) 


where G(p) is a gate function describing the limits of the illuminated area. Specifically 


G( p) 


1 if p£ A/tt 

0 Pi: 


(A-52) 


where A is the area of illumination in the mean plane of the surface. Now recall that 


rat 

I exp(±j«si 


jasin^)d|= 2Tr0o(“) 


The integral- will then reduce to 


2 2 

An asymptotic evaluation of the above integral for large yields 

-tan^072m^ 
irAe ^ 

<II*> - g”g 2 
2k m COS 6^ 

2 .2 

where o / p“(0) has bsen identified as the slope variance m , 

Combining the above results it is clear that 

* B-.bJ 

= — 

where j, k = v or h , 


(A-53) 


<II*> = ZTiPiJ 


(A -55) 


(A-56) 
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4„0 THE ANGULAR COHERENCY OF THE SCATTERED FIELDS 


At this point consider the mutual coherence function 


9p<?3)E*{e2,«3)> 


{A-57) 


where i(j) = v or h. The coherence function denotes the cross-correlation between two 
field components scattered in the plane of incidence at scattering angles 0^ and 02 with 
a common range R, The expectation is an ensemble average over all random surfaces 
satisfying the Kirchhoff approximation. Now in view of (A- 39) the coherence function 
can be written in the form 




where 


'jk 




(A -58) 
(A -59) 


Now transform the center of mass coordinate system where 

“ " (A-60) 

V = y^-yg 


Also let 0^ = 0. and 02 = 6. + A0 where 40 is a small deviation fi;om the backscatter 
direction. We have 


(nj-n^)*Pl 


-2(s1ne,yi + cose^Zj) 

(-2sine. + Aecos0^)y2 + (2cos 0^- Aesin 0-)Z2 


Now for a gaussian random surface having a surface height characteristic with zero mean, 

O 

variance tr and correlation function A, the expectation within the integral becomes 

< > = expf-jZkyV+jk^Aeyg I exp [-ky^ar^A6^/2 ] exp[-4kj,^o-^{l-A)l 

(A-62) 

The integral then can be written as 


Intg 


exp[-k^VA0^/2l 


//// 


exp[-j2kyV + jk^AGyg]. 

exp dudvdXgdyg 


(A-63) 


I27 


Now fransform the integral expression to cylindrical coordinates by letting 

U = PCQS^ 

V - PS inf (A ”64) 

Xg = p'cosf 

yg = psinf 


Then 


^jii(Oi,02) = I E^lfixp(-k20^AQV2)Intg 


where 


(A ”65) 


Intg 


■ ////‘ (p)G(p')exp [-jEkySinfpJexpfjk^AQsinCP'l 


exp[2k ^o-2(i-A(p))]pdpdfP'dpd{ (A-66) 


and where it has been assumed that z is stationary and isotropic, G(p) is a gate function 
defining the region of uniform illumination on the mean plane. Using Equation A“53 
twice, we can write the integral as 


where 


Intg = (27r)^I^Ig 

(A”67) 



y 6(p) (2k^p) exp E4k^^or^(l-A)]p dp 

(A-68) 


h - 

/g(p00o (kjj Aep^p'dp' 



For a circularly illuminated area of radius R the latter integral can be evaluated to 
get 


where Jinc(x) = 


we have 



1 2 = R^Jinc(k^R 60) (A”69) 

Jj{x)/x. , As a consequence, for o circular region of .area A 

= 4 A lKi='C^j^exp{-kj0^A6?/2)Jinx(k^R^A0)I 


(A”70) 
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Now normalize the mutual coherence function in the following fashion 


r 




4TrR / ( I c’ I ^ Acos 
o ' o' 



(A-71) 


Now recognize that for small j10 

a 2exp(-kyVAS^/2)Jinc(x)47T<S^jS*j,> 

where x = equations A-43 and A- 54) . The degree of coherence or 

correlation is consequently related to the mutual coherency function by 


D = > (A-72) 

or 

= 2 exp{“k5a^A0^/2)Jinc(k R A6) 

Consider the character of the degree of coherency. Except for extremely large 
ko- -values the exponential term contributes negligibly to D at small incident angles. 
The decorrelation is consequently largely governed by the Bessel function for small 
incident angles, D vanishes at the zeroes of Jinc. The first zero occurs where 

AQkR^cosO = 3.832 ‘ (A-73) 

The corresponding angular separation is given by 


A0 = 3.832/kR cos0 
0 


(A-74) 


Suppose k = 291 (f = 13,9 GHz), = 10 meters and 9 = 25 *^; then decorrelation occurs 
when A9 = 0.00146 radians or at .,08 degrees. It is concluded that radar returns de" 
correlate rapidly with changes in view angles. 

At large incident angle (grazing angles) the exponential factor will pre- 
dominate, This result is physically reasonable since the surface roughness predominates 
the view; whereas at small angles the area of illumination as conveyed in Jinc is the 
dominant factor. 
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APPENDIX B 

The Scatterometer Equation Within 
the Context* of a Scat't'ering Theory 


loO INTRODUCTION 


The scatterometer equation Is once again derived within the context of a 
specific scattering theory. The structure otid meaning of the formulation, as a 
consequence, readily becomes apparent. Specif'cally, the angular correlation 
assumption is shown to be equivalent to the non-ccherent property of scattering; the 
relation between scattering operator and the scattering coefficient is clarified also. 


1.1 Derivation and Discussion 


Silver [19] has shown that the far field radiated in the direction n^ from a 


bounded closed surface S having surface excitation E and H is given by 

„-jkR = 


E = 
s 


47tR 


iig x: jy hxE-niigxCnxH) 


jkp*n 


s ds 
(B-1) 


where 

R = distance from the surface to the far field point 

P~ position vector from an origin local to the surface to a point on the surface 
n - surface unit normal 
k = 27 t/?l 

TJ = intrinsic impedance 

The geometrical entries of the above expression are illustrated in Figure B~1 . 

Suppose that the surface is smoothly u ndulating and perfectly conducting. Then under 
the Kirchhoff approximation the tangent surface fields are given by 


n X e: == 0 (B-2) 

and 

hxfi = 2nxH^ (S-3) 
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where Hj. Is transmitfed field incident on the surface. The scattered field, therefore, 
sirhplifies to 




jkn 

2ttR 


hg >« JJ ( n ds 


(B-4)- 


Now when specializing to the backscatter case we can write 


Hg ^ ( fig ^ ( n X )) = - { hg • h ) hg X 


so that 


-jke 


-JkR 


2ttR 




) ds 


where for the backscatter case 

= n X H 


(B-5) 


(B-6) 


(B-7) 


s ” "t 

In beginning with the far field expression we have in effect anticipated the use of a non- 
coherent assumption since in scatterometry one would not necessarily be in the far field 
if the entire scattering aperture were coherent. 

Now the spherical incident field is denoted as 


where 


= -a K 0>t - kR ) 

K = i^4n'R 

ij^ = antenno input current 

complex effective height vector of the antenna 


(B-8) 

(B-9) 


The range R is measured from the antenna as illustrated in Figure B-1 . The spherical 
wave can be approximated by segments of plane waves each illuminating a patch of the 
surface. In addition, the incident field amplitude components as conveyed by may 
be considered constant on a given patch. Suppose the entire illuminated surface is 
segmented, into N patches. Then the incident field on the mth patch may be approximated 
by 


E = _i K r pj( tut - kR + k p *n ) 


(B-10) 
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where 


R range to the centroid of the mth patch 

= position vector from the centroid to a surface point In the mth patch 

n = unit vector in the backscatter direction for the mth patch 

Now from Equation (B-6) we note that the backscattered field for the mth patch is 
given by 




.-j 2kR 


"sm 


2ttRi 


‘m » 
L 


‘m 


tm 


//' "'-'sn. 


\ -j 2k p*h j 
; sm ds 

(B-n) 

The integration Is performed on the surface within patch m. 

Now define the scattered field per differential steradian subtended about the 


m 


antenna as • 


t rI 

sm m 


sm 


Sm 


(B-12) 


where /iA is the area of the mth patch in the mean plane of the surface and 0 is the 
m ' m 

incident angle on the mth patch. The scattering operator employed within Chapter 4 
may be identified with the above expression, viz,. 


S ( 6 ,(J) )= liin 

vv ' m’^m ^ 


sm 




-j ’ vttS 


4-itR 


(B-13) 


where 


hh ^ ^ f sm ‘ jim 


’vtnT '-t ■ ’em 




(B-14) 


(|)m 


The inner products above isolate the vertically and horizontally polarized components as 
defined with respect to surface. It is appropriate to denote^^ os an operator since it 

Hr 

must recognize the phase of the incident field relative to the centroid of the patch . 
Whether these differential scattering operators exist is not important to the development 
within the appendix. However, within the main text they are assumed to exist at least 
in approximate form. 
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Within the context of this theory “-^hv simply a statement of a well 

known result that a smoothly undulating perfectly conducting surface does not depolarize 
the incident field. In general, the latter operators are not zero. The above increment- 
al field is defined so that the total field at the antenna is given by 


or in the continuum limit 


N 


E, = E Afl. 


Ea = 


tn=T^ sm m 

ds 


■/ 


(B-15) 

(B-16) 


It must be recognized that the antenna does not respond to the total field, a 
quantity often computed by the theorist. Instead the antenna responds so that the open 
circuit voltage induced into the antenna terminals by the mth patch is given by 

~ '^stn * ^rm (B**17) 

where is the complex effective height vector in the direction of the mth patch 
during reception. The total induced voltage is clearly approximated by 

The average power available under matched conditions at the antenna terminals is given 
by the ensemble average 

^ SRy. ^ (B-20) 

where is the radiation resistance of the antenna during reception. 

Now the scattered fields and effective heights can be decomposed into polarization 
components coinciding with the surface polarizations, viz,, 


as^d 


E = E 

sm smv 





(B-21) 


L 

rm 


= 1 


vrm 



+ 1 


hrm 



(B-22) 


134 


The expecfal'ton can then be written in the form 


< £ • 
SID 


or 


— * 

L £ " 

rm sm 


> = < ( £ 
rn ' vsm 


= tr M M 
rmn smn 


\rtti "** ^hsn ^hrn ^ ^ ^vsn \rn ^ ^hsn ^hsn ^ ^ 

(B-23) 

(B-24) 


where 




^vsm 

£ vsn ^ 

< £ 

vsm 

^hsn ^ 


M 

smn 

/ ^hsm 

£ vsn =■ 

< e 

hsm 


and 


r ^vrm 

* 

Vrn 

1 

vrm 

* 

^hrn 

1 


(B“25) 


M 


,rmn 


^hrm ^vrn 


^hrm ^hrn J 


(B-26) 


are identified as mutual coherence matrices. The mutual coherence matrix for the 
scattered field is composed of elements correlating fields arriving from patches m arid 
n or equivalently from different angular directions 0 ^ and (0^, 0^). Now within 
Appendix A it is shown it is reasonable to assume, on pragmatic grounds, that returns 
arriving at different view angles from the same patch are uncorrelated. The assumption 
is exact in the geometric- optics limit. It is even more reasonable to assume here that 
fields arriving from different angles are uncorrelafed since they arise from different 
patches. As a consequence the mutual coherence matrices have the special property 


M 


sjk 


= N .. 6 ., 
sjj jk 


(B“27) 


for every jk. is the Kronecker delta and 


^sjj 




SJ 




'hsj. 


<£ 


hsj 


£ * . > 
^ vsj 


i^hsi > 


(B-28) 


The non-coherent assumption, consequently, allows us to write the received power in the 

form N N t 

tr E ( S Kiu '^ik ^^k ^ ^si.i (B-29) 


W,= 


8Rv 
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or 


8R* 




where 


>^J rjj 


(B-30) 


(B-31) 


is the coherency matrix for the receiving antenna in direction (6,, 0 .). Now expanding 

* .A ■ 


N aH.,, we have 
sjj J" 


4 k^jK/R? 


N„*. M2. “ ^ 

(4Tr)^ AA- cos 0. 


'jvv 

L^jhv 


B- u 

jvh 


B 


jhh 


(B-32) 


where 

B * “ 1 nf 1 ^ 

Jpq 


qtj • "sj ^ ■ ” "sj^> <iS dS 




(B-33) 


o 

From the obove expression we note that the. incident field complex amplitude components 
and I have been separated from the scattering integrals. The relative phase 
between incident amplitude components is retained in the products ^'he 

ncxi-coherent differential scattering coefficients per unit steradian for the jth patch is 
defined as 

< § > = 

pp qq 


4k' 


(4tt)' 


jj< (^*nsj)(n'-ngj) 




sj “ > dS dS' 

(B-34) 


(See Equation (A“43), Appendix A). As a consequence we can write the coherency 
matrix for the scattered fields based on intensities per unit steradian in notation similar 
to that of Equation (4-16) of Chapter 4, We have 


or 


'^sj ° ”sjj “j 


= |K/ 




vv ^ vtj 1 htj 


(B-35) 


^vv ^hh ^ ^ vt j ^ htj 




(B-36) 
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We have shown above thai- the elements of the coherency matrix change units 
after the non-coherent assumption has been applied to the double summation. The 
return power can now be written as 

H - ^ E tr c^. cj. An. (B-37) 

or upon taking the limit of the sum as we get an integral approximation 

” ° "sljr / tr dS! (B-33) 

where the above equation has a form identical to that of Equation (4" 13) of Chapter 4 
when^^j^ = 0. 


APPENDIX C 

Correlation oncf Cross-Correlation 


Scattering Properties oF a Slightly Rough Surface 


UO THEORETICAL DEVELOPMENT 


For a plane wave incident with angle 0. in the x - z plane on slightly rough 
surface, satisfying the requirement 

(kcycos9)^«i (C-1) 

2 

where cr = surface height variance and k = 2ir/x, the solution for the scattered fields is 
expressed in a perturbation expansion of spatial Fourier components [32] [33] [34] [35] , 
The Foi'Tier components are interpreted as an angular spectrum of plane waves. Suppose, 
the spectral components are denoted as A'^ (k , k ) where 


X 

= ksine^cos?^ 


ky = ksin0^s1n4>g 


(C-2) 


p^ x,y,z 

Figure C-1, Scattering Geometry 


z 
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when a verHcatly polarized plane wave illuminafes the surface and as A (x,y,z) when 

a horizontally polaiized plane illuminal'es the surface, (See Figure C-1). For either 

case the electric field E at point (x,y,z) is given by an inverse Fourier transform re- 

P 

latDonship 


(27t)‘ 


A (k , k ) e 

D X y 


+ kyY - k^z) 


(C-3) 


o 5^ j? 2 

where k = k - k - k . The superscripts denote the incident polarization whereas 
the subscripts denote the scattered cartesian components (E , E , E ), If we suppose 
that G highly directional antenna points in the direction {0^/ 0j with beam^volume 
A0A^^ then only certain spectral components near 


= ksin6 COS0 

A S S 

ky = ksin6^sirnj)g 

k, = kcos0 
2 s 


(C-4) 


will be observed. The angular spectral space surrounding (0^,^^) with angular volume 
sin9^AQ^ is given by 




J I A0A<f> 


(C-5) 


where 


J = 


3kj 

d0( 

Sk. 






2 

k sin0^ cos0^ 
s s 


(C-6) 


Now the cross” correlation of electric field components observed within the 
antenna beamwidth is given by 

<ae;ae^^% — ^ 


(2^) 



A > 
p q 


AxA 


J(k^x + kyy) + /y) 


dk dk d«: dfc 
X y X y 

(C-7) 


where- 
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AxA = (Ak ) X (^i< Ak ) denotes the cartesian domain of intcgrationo The solution 

X y ^ y 

for the spectral components are given by the form [35] 

"" A„p’^(k,k')Z(kj^+ks1n0o, ky) (C-8) 

for the first order solutions and by the forms 


= bI 


x,y X' y 


°x,y 

r 

^x,y 

r 

'x„v 

'x.y 



''x,y Z(iV“.y 6)dadS + 





f, (l<, k) Z( a+ ksine , 0 ) Z(k - a, k - 0 )da ds + 

i X y 

fp(k, k') Z( a+ ksin 0 , B ) Z(k„- a , k - B )da dg -J- 
t- X y 

^y»x ’ V ^ 


and 


" “* + — B ^ 

^ ' k, 


(C-9) 


(C-10) 


for the second order fields, ky) is a random vari able describing the Fourier spectral 

heights of the rough surface and =\/k^ " - ky^o The coefficients ^ 

Bn and Bn ^ are determinstic functions depending on the propagation constants in 

the upper and lower medio . For the purposes of this appendix it is sufficient to observe 

that the first order fields ore dependent on Z(k , k ) whereas the second order fields 

X y 

are dependent on Zfa+ksin9^. B) l< "8).. The functional form of the A and 

B coefficients is at this point immaterial ^ 


2.0 CORRELATION PRODUCTS FROM THE FIRST ORDER FIELDS - CASE I 


Now consider a typical correlation product from the first order solotlonso V/e 
have from the integrand of (C~5) 

= I = A^„(k J') A® * («,«') 


p' X y' p' X y' 


<Z(k^,-fksin0o, k^,)Z“(fCj^+ksin0o, > 

(C-11) 
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The expectaHon can be writfen as 

<Z(u ,v )Z (u‘,v'}> = /y//^ 2 (x ,y )z(x\y’)>. (C- 12 ) 

J J J . -j(ux + vy) j(u'x’+v’y') 

e e dx dy dx'dy' 

where z and Z are Fourier transform pairs. If the random process is stationary, then 


<z{x ,y ) z(x',v')> = R(x-x‘, v-y') 


(C- 13 ) 


where R is the auto- correlation function of the surface heights. Now transform the 
variable to the cen 
and y' = y' . Then 


variable to the center of moss coordinates, i.e,, leti = x " x', t - y ' y'/ x' - x 

X y 


<Z{u ,v )Z (u',v')> = 

j[(u‘-u}x' + (v'-v)y'J 

2 dx ' dv ' 


or 


- fful \ -j[(u-u')x' + (v-v')y'] 

" J J e *- J dx'dy' 

(C- 14 ) 

where W{u,v) is the spectral density of the surface heights. Finally we recognize that 


<Z(u ,v )Z (u‘,v')> = (Ztt)^ W(u ,v ) S(u-u ' ,v-v' ) /p,. 


(C- 15 ) 


When the above results is substituted into the correlation integral we get 


< ah’" AE*^ 
p q 


/ j % (C-16) 


2 / f ’p °q 


When A is sufficiently small, the Incremental complex intensity may be approximated by 


<ae!'ae*®> = 

p q 


A*" 

«op ^oq 


W(kj^+ksin0^, k M^Ak (c-17) 


(Ztt)" ^ - y x y 

The cross'correlation per unit steradlan is therefore given by 


<ae|!’ AE*^> 

p q 

s1n0gA6^A4>g 


2 v» *c 

k^cos0^A' 
s °p 

(2x)^ 


W{k^+ ksine. ky) 


I 4 I 



where (C-5) has been used. The total cross’'correlation in direction (0 , ) is given by 


^n> ^ 

p q 


k^cos^Gg A^p Af* W(kj^ + ksinGo, k ) A A 


(2tt)‘ 


R 


2 (C-19) 


where AA is the illuminated area in the x “ y plane and R is the range to the element of 
area. 

To consider the horizontally and vertically polarized backscattered components 
we must realize that the horizontal polarized component is related to the cartesian 
components by 


AjJ(kjj,ky) = -aJJ (kj^,ky) + A^ coscji^ (c-20) 


and the vertical component by 

^''x’ '‘y' " ^x ^’‘x’^y' cosG^cos*^ + aJ (kj^, k^) cose^sln^ij - 

(ky.ky) sin9j (C-21) 

When we specialize to the backscatter direction (6^ “ 0^, '^)f we see that 


and 


''h ('x’",’ 


(C-22) 

a; (kjj.ky) 


(C-23) 


with k = ksin0 y k = 0, k = kcos0 , It is well known that the first order backscatter 
X o y z o 

fields do not involve depolarized components [34] . From reference [35] we now identify 
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''h ' -j2kcose„Rj^ Z(2ksinS„,0) (C-24) 


and 


where 


(''x’l'y) = -j2kcose,T^[e^(l+ sin\) - sin^e„]Z(2k5ln0o.O) 

(C-25) 

= Fresnel reflection coefficient for horizontal polarization 


T - Fresnel transmission coefficient for vertical polarization 
V 

= relative (complex) dielectric constant. 

and where a unit amplitude incident wave has been assumed. Therefore the correlation 
components in the backscatter direction become 

|2 AA 


V V 


E*^*> 


4 4 

4k cos 6 


/i. 4 4 

4k cos 0. 


> 1 r ? 9 1 ^ 

jTy[€^(l+ sin ) - sin VI(2ksin8^. 0)- 


(C-26) 


{27t)' 


\ e^{l+ sin\) - Sin\ R^j W(2ksin9,, 0) 


AA 


<E^ E^*> 


4k^cos^e „ AA 

— |Rj^r W(2ksine„, 0)— 

( 2 tt)^ ^ 


(C-27) 


(C-28) 


The corresponding generalized differential scattering coefficient per unit intensity per 
steradian is given by 

,2 

(C-29) 


<S*^ S*^> = <E>" pS* 

pp qq ^ tp tq >■ 




AAcos 0, 


(See the generalized definition of the scattering coefficients in Chapter 4). So 

I 

4 


<|S, 

I vv 


k 2 

cos^e„ |T„[€^(1+ sin\)- sin\]| W(2ksin0o> 0) 

(C-30) 


Tf- 


O' V r 
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k^cos^e 




[Ty[ sin\)- sin\]jR* W(2ks1n0^, 0) 

(C-31) 


4 3 

o k cos 0 o 

<|SJ> = iRj^r W(2ksin0,. 0) 

tf2 


(C-32) 


3.0 CORRELATION BETWEEN FIRST AND SECOND ORDER FIELDS - CASE II 


To develop the correlations between first and second order fields it is sufficient 
to note that the correlations involve expectations of the type < Z(k^,ky) Z* (a, ^ Z* (K^ 
'’and of type Z(kj^,k^) Z*(a+ksin0,3) Z*(K^-a , , These 

expectations involve independent gausslan random variables with zero mean and con- 
sequently vanish. The first and second order fields are therefore uncorrelated „ It is 
concluded that = 0 and = 0 at the lowest order. 


144 


APPENDIX D 

5ca^^■erometer StmuIaHon Program (SCATSIM) 


1 ,0 INTRODUCTION 


The theory and operation of the scatterometer simulation program is described 
within this appendix. The following section shows how the scatterometer equation of 
Chapter 4 was implemented with idea! and non” ideal antenna parameters. The compu” 
toHon of the inversion models with and without recognition of the difference between 
surface and antenna polarizations is also described. Finally the operation of the pro- 
gram Is treated by means of a flow chart. A source listing and a sample output is also 
presented. Additional program documentation is provided by comments within the pro- 
gram . 


2.0 THEORY 


2.1 Simulation of the Scatterometer Equation 

Since the scattering characteristics were based on surface polarizations. Equations 
(6-47), (6“49) and (6-50a) through (6-50i) were implemented for use on the computer. The 
equation was simulated using identical functional forms for the verticolly and horizon- 
tally polarized patterns (if they are both present during a transmission or reception). 

Recall that the normalized patterns are given by 


'v»h 




* |1|,(0,0) 


P-1) 


where (0‘, 0 *) = (0, 0) is the boresight point. As a consequence, we require 

+ g^(0,0) = 1 (D-2) 
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Now if g denotes the functional form for the pattern and has the property g(0,0) = 1 
and if gj^ is assignee' the value ag where a <1 , then we require that g^ = (1 - a)g„ The 
scatterometer equation under the above assumption can be written as 


W(0o) 


(4 IT 2 ) 


/ 


i^p(g cose)'^ do. 


where 


(D-3) 






21? + 2IaRe<S^^S*^> - 2IgIn,<S^^S*^> 

(D-4) 


where 


= ( 1-a^) ( l-a^)cos^(i^ + a^a^sin% 

(l-a|)a|. + 4c^c^) sin^i|icos^ijj 

^2 "" a^a^cos% + (l-a^)(l-a^)sin^ 

(l-a^)a^ + 4c^c^)sin^fJJCos^ijj 

I3 - ( (l-a^)a^ + (l-a^)a^ -2c^c^)(cos% + sin%) + 

2((l-a^)(l-a^) + a^a^ - ^ (2a^-l) {2a^»l) ) * 

s1n^i|jcos^i{/ + 2s^s^ (D”?) 


+ (d-a )ax + 

(D-5) 


^ + (d-a^)a^ + 

(D-6) 
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= C^c (oos^il^ + Sin^ip) + ( (2a^-l) (2a^-l) -2c^c^)“ 

2 2 
sin ipcos ^ 


I 5 ' 

( ‘=r=t ‘ ' 

(D“9) 

'6 = 

(d-a^)c^ + 'T-a^)c^) cos'll!) - (a^c^ + a^c^) 

Sin% + 3 ((2a„-l)c {2a^-l)c,l C 0 S^ 4 r 

(D-10) 

^7 

((l-a^)s^ + (l-a^)s^) cos^ilJ + (a^s^ + 

0 

) sin ijj 

p-11) 

II 

00 
t— t. 

(^tC^ + cos% - ((l-a^)c^ + (l-a^)c^) 

sin^ij) - 3 ((2a^-l)c^ + (2a^-l)c^) cos^ij; sin 

|)-12) 

'9 = 

^ (d-a^Js^ + (l-a^)s^) 

- 2, 
sin i|j 

(D-13) 


■ = /a\j.U - a^) cosB^ 

- a^) cosg^ 

(D-14) 


Sp = »^ay.U - a^J sinB^ 


The infegration is performed in the surface coordinate system (See Figure 4.1), In the 
above expression all odd powers in sinip have been pragmatically dropped. These factors 
are odd functions of ^ and will not contribute to the integral (See Equation 4~33), 
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Now since the scattering characteristic was assumed isotropic, the return 


power can be aporoximated by 

W(0„) = tr y 


(4ttz ) 


^vvl ^^(0 J cosef dsi 
2 ^ 


_ f “u 

I ? I 2 

/ (O 

IgCg coi:0f d 


+ < I s 


- 2 Im<| 


* r / ^ 

vv^hh^u cose) 

i) 


+ 2 Re <s s 
vv 


hv^w / ^6^ 

*/o. . 


g cose) dn 


- 2 Im f ijig cosef dS! 


+ 2 Re <Sy^s 


- 2 Im <s . s. . 

vh hh 


* 2 

hh^u / 

^Oi 2 

hh^oi j ^9^3 "“SJ 

»/r\ 


p-iS) 


and where ^ ^ = 1 ,2, o . .,N ^ Is a set of half degree annuli centered about the sub- 

observation point. < denotes an evaluation of the scattering coefficient on the 

J2 annulus. Function subroutine SIGAAA contains the functional representations for all 
(0 

nine scattering coefficients. The integrations are performed, of course, only over those 
annuli where the pattern function g is significant. The above approximation reduces 
the computation to integrals of the following kinds: 


2 2 4 

g cos e cos !/; dS^ 


■I 




f 2 9 4 

= I g cos e sin dQ ^ 

J n 


Jy,{m) = 


2 2-2 2 

g cos 0 sin if; cos t|j dJ2 


cos^e 


J ( w) 


L 


2 2 2 
g cos 0 cos ij; dn 




■L 


2 2 
g cos e dfi 


w ''to (D“16) 

These integrals are evaluated in subroutine DINTEG using a two“dimensiona! Gaussian- 
Legendre quadrature technique [48] . For a selected antenna pattern and a selected 
view angle the return power is computed in accord with the above expression. The antenna 
gains G^ and G are formed in a separate computation. These factors are based on an 
evaluation of the expiession 

2 

(D-17) 




J*g sin0* d0‘ 


The evaluation of the above integral is performed in subroutine SOLID, which employs 
a single dimension Gaussian-Legendre quadrature. The numerical evaluation of the 
pattern functions is provided by subroutine LAMBDA ^ All of the above integrations are 
executed from the mainline of SCATSIM, 

Since the relative phases 3^ and 3^ were assumed stationary over the main beam 
and first side lobe, the return power could be evaluated for various combinations of a^, 
a^., 3^ and 3^ without re“eva!uating the double integrals. As a consequence, an arbitrary 
pattern condition within the above constraints (g^= (l"a)g and g^^ = ag) could be esta- 
blished. The combination of relative amplitudes and phases for the fifteen prescribed 
measurements are shown in Table D*1 . Subroutine ANTENNA, when addressed with 
zero arguments, generates those prescribed values. When amplitude and phase biases 
and/or perturbations are entered as arguments, subroutine ANTENNA will apply biases 
of the prescribed value to all measurements In which a^ or a^ is zero or unity. Random 
perturbations are applied to the remaining coses if the perturbation arguments ore non- 
zero, In this fashion either measurements based on 15 ideal or 15 deviated antenna con- 
ditions can be generated. The actual coefficients required in the integrand factors 

‘ 9|are computed in subroutine COEF. COEF fills a 15 x 9 x 6 array with the 
appropriate values so that the return power can be computed for each of the 15 measure- 
ments. Let C.jj^ denote the array. The i subscript designates the measurement number, 
the j subscript identifies one of the nine scattering coefficients within the integrand, 
and the k subscript identifies one of the six kinds of integrands (Jj^(ui)}. See Table D.2 
for the entries in Let Yp i = 1,9 denote the nine scattering coefficients and 


f 


149 


TABLE D.l 


MEASUREMENT 

NO. 

COEF 


St 

a 

r 

0r 

1 


0 

— 

0 

— 

2 

<'=hh''> 

1 

“ 

1 

— 

3 


0 

— 

1 

— 

4 

Re<S S,,*> 
w hh 

0.5 

-90° 

0.5 

90° 

5 

Re <S S, , *> 
w hh 

0.5 

0° 

0,5 

180° 

6 


0.5 

45° 

0.5 

-135° 

7 


0,5 

-45° 

QZ5 

135° 

8 ■ 

Re<S St *> 
w hv 

0 

— 

0.5 

0° 

‘ 9 

Re<S S, *> 
w hv 

0 

— 

0.5 

180° 

10 

Im<S St *> 
w hv 

0 

— 

0.5 

90° 

n 

Im <S S, *> 
w hv 

0 

_ 

0.5 

-90° 

12 


1 

— 

0.5 

0° 

13 


1 

— 

0.5 

180° 

14 


1 

— 

0.5 

90° 

15 

V> 

1 

— 

0,5 

-90° 
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T.^RLE D-2 

THE COEFFICIENT MATRIX C,., 

Ijk 


ROW/COL, 1 '* 2 3 4 5 6 


1 

I 

A 

°r°t 

(1-a )o + (1-a }a 
. r t t r 

+4c.c 
t r 

0 

0 

O 

0 

2 

Vt 

a-Of) (1-a^) 

(1-a,) Of +(1-0^0^ 

+4c,c^ 

0 

0 

0 

3 

{l-a^)a^+ (l-a^)o^ 

“2c.e 
t r 

-2c, 

(1-ap {i“o^)+ aj.a^ + 
(2o^-l)(2a^-l)-6c^C|. 

0 

0 


4 

2c,c 
t r 

2c.c 
t r 

2 [ {2o^-1) (2a^-l)-2c^c^J 

0 

0 

"Vr 

— 5 

Cn ^ 

0 

0 

0 

-2(V,+',V 

2(CrS,+ c,5,) 

0 

6 

21(1-o,)c,+ (1-0,)=, 1 

-2[o^Cj^+a^c^l 

6[(2a^-l)cj+(2a^-l)c^] 

0 

0 

0 

7 

0 

0 

0 

-2 [ (l-a^)sj + 

(l-a^s^ 1 

-2 [0,5,+ 0,5,1 

0 

8 

2 [a,c + a c ) 
t r r t 

-2 t{l-aj.)c^+ (l-aj,)c^l 

-6l{2a^-l)c^+(2a^-l)c^l 

0 

0 

0 

9 

0 

0 

0 

-2 (a.s H- a s ) 
' t r r t 

-2I(1”o^)5^ + 

(l-a^)s,l 

0 


NOTEj The different planes of C,.j^ I = 1,2, 15, are 
formed by substituting values for a^, and 0^ 

from Table D.2, 
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let- Y.( w) denote its evaluation on the annulus. Furthermore let 

N 

“ S Y](w) J[^(a)) 

u=l 


where denotes the evaluation of the J^^annulus. Then the return power is 
given by 


W^-(0o) 


(47rz ) ' 


tr 


6 

Lk=l 


jk 


(D-19) 




Subroutine EXACT performs the above computation for i - i,2, ,o. .15, For each 
measurement (i) the contribution by each Y. is isolated by EXACT and stored in its second 
argument. The power matrix is clearly given by 



X^G^G 
t r 

(4itz)^ 


6 


(D-20) 


The structure of the return can thus be examined , 

Another routine, called IDEAL, also estimates the return power but without regard 
to the distinction between surface and antenna polarizations. The computation follows the 
above scheme, however, it recognizes that = 0, In this case, all the necessary inform™ 
ation is carried in Again IDEAL isolates the power contributions by each scatter" 

Ing coefficient and consequently forms a power matrix also. 


2,'' The Inversion Models 

The above formulation of the return power was designed so that the exact or 
approximate inversion model parameters could be isolated from intermediate steps. 
Either model assumes that the measurement can be approximated by 



where M.. is a 15 x 9 matrix. Each row of M corresponds to one of the fifteen measure” 
ments. When the distinction between surface and emtenna polarizations is required, the 
elements of M.. are simply filled by forming 

6 N 




k=l 


j k 




w=l 


(D-22) 


M.. is constructed in subroutine EXACT. Once M.. has been constructed the inverse 
'J U 

model is computed by forming the normal equations M.. M.j^and computing its inverse 

OdEMINV in the mainline). The inversion for this mode! is performed In subroutine 

MATRIX. Within MATRIX a least squares solution is executed, viz,. 


y;(6o) = w,(9j 


(D-23) 


This subroutine also accumulates the first and second order error statistics during a Monte 
Carlo study, A call to a secondary entry MATSHOW will display the statistical results. 
Similarly when the distinction between polarization frames is not required, an 
approximate inversion model may be formed from Equation D-22 above by simply setting 
= 0 in each {See Equation D-16). Symbolically we have 

6 N 


13 


k=l u)=l 


(D-24) 


The integrals need not be re-evaluated since all the desired information is contained in 

Particularly tjj = 0) = ( <4 ijj = 0) = d^(w). The remaining J are identically 

zero. These special properties were recognized and accordingly a routine (IDEAL) was 
prepared to evaluate the elements of M.j for this case. The inversion of this model is 
performed as suggested in Chapter 6. Recall that the ^ 

<{ ^vh^^^ each computed from a single observation (a row of M). The remaining 
coefficients are computed by differencing pairs of equations (rows). The inversion for 
this model is performed in subroutine DIFFER, Again first and second order error statistics 
are accumulated. They are displayed by calling the secondary entry DIFSHOW, 

The reader should note that routines EXACT and IDEAL play dual roles. Either 
can form their respective inversion models or they can compute the return powers for 
the fifteen kinds of measurements. Only EXACT computes the exact return power since 
the scattering coefficients are defined with respect to the surface polarizations. The 
option to use IDEAL to compute return power exists to compare the two polarization frames. 
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2,3 Documentation for SCATSIM 

A mocro“flow chart of program SCATSIM is shown in Figures D,l through D,4, 

The program is organized into roughly four functions. Each figure covers one of the 
program functions. Figure D.l documents the part of the program which reads the in- 
struction card and initializes various parameters for use In the actual simulations. This 
portion of the program, once validating the instruction, prepares various descriptive 
antenna parameters such as beamwidfh and gain. These parameters and the input para- 
meters are displayed to document the case study. The second zero in the pattern function 
is employed to establish the domain of integration. The domain is broken into N half 
degree annuli . On each annulus k = and is computed . 

Once are formed, dj^(^) and E i|J = 0) are computed and scaled 

for the antenna gain effect. ^ ^ 

Following the above initialization, the program simulates the fifteen measure- 
ments under ideal antenna specifications. The structure of this program portion is 
illustrated in Figure D.2. Subroutine ANTENNA is called with zero arguments to pre- 
pare an ideal antenna parameter set {a^ a^, i = 1 ,2, . , . .,15, Subroutine 

IDEAL forms the approximate (ideal) inversion model and a power matrix which ignores 
the distinction between surface and antenna polarizations. Both the inversion model 
and the power matrix are displayed. The inversion mode! is stored for subsequent use 
by DIFFER. Subroutine COEF forms C.^j^ from the ideal antenna parameter set. In turn, 
exact then uses C.j|^ to compute the exact inversion model, M... It is again called to 
form the exact power matrix. The normal equations are prepared* from the model and 
then is inverted by HEMINV, If the system is singular, a flag (ISING) is set true. All 
matrix inversions are subsequently by-passed by an appropriate test. The exact inversion 
model and the inverse of the norma! equations are stored for subsequent use by MATRIX, 
Once the return power is computed, both difference and matrix inversions are performed 
and the statistical (accuracies) results are shown for the ideal antenna. 

If the bias parameter ABIAS is non-zero, a bias error study Is performed. This 
portion of the program is illustrated in Figure D.3. The processing follows, for the most 
part, that performed in characterization of the ideal antenna; however, the two bias 
parameters ABIAS and BBIAS are employed in the arguments of ANTENNA to introduce 
pattern deviations from the ideal case. The inversions are perfornjed using the ideal 
antenna models. | 

If the perturbation parameter AAMX is non-zero a Monte-Carlo study is perfotmed, 
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INPUT PARAMETERS 


ANTENNA TYPE - ITYPE = 1,2,3, or 4 
BEAMWIDTH PARAMETER - ka 
VIEW ANGLE - tnot {9^) 

AMPLITUDE BIAS - ABiAS 
PHASE BIAS - BBIAS 

MAXIMUM AMPLITUDE PERTURBATION (AMAX) 
MAXIMUM PHASE PERTURBATION (BMAX) 



Figure D Ja “ MACRO-FLOW CHART FOR SCATSiM - 
PROGRAM iNJTIALIZATION AND ANTENNA PARAMETERIZATIONS 
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Figure D . 1 b - MACRO-FLOW CHART FOR SCATSIM - 
PROGRAM INITIALIZATION AND ANTENNA PARAMETERIZATIONS 


!56 






Figure D,2a - MACRO-FLOW CHART OF SCATSIM- COMPUTATION 
OF THE IDEAL ANTENNA RESPONSE 
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Figure E).2b - MACRO-FLOW CHART OF SCATSIM - COMPUTATION 
OF THE IDEAL ANTENNA RESPONSE 




Figure D.3 - MACRO FLOW CHART FOR SCATSIM - 
COMPUTATION OF THE BIAS RESPONSE 
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Figure D.4 - AAACRO FLOW CHART FOR SCATSIM - 
COMPUTATION OF THE MONTE CARLO STUDY 


l<50 







The documenl-aHon for this portion of the program is illustrated in Figure D.4, The 
course of this portion of the program is identical to the bias study except that many cases 
are examined. The number of cases is specified on the instruction card. For each case 
the antenna parameter set is perturbed randomly within subroutine ANTENNA, 

2.4 Program Listing and Sample Output 

The source listing of SCATSIM is shown in Figures D.5 through D.18. Sufficient 
comments have been inserted to identify variables with the theory and to track the 
operation of the program, A sample output is shown in Figures D,19a through D,19f, 
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SCAT SIHULATIQN POOGRAH 


LAS EL HAIM PAGE 


SCAT SIHOUATION PfiOCRAH 

THIS PPOGRAH EflAOLtC TilE USER TO STUDY THE PERFOR- 
HADCE OF HIS SCATTEROHETER AHTESMA hHEH IT POSSESSES 
leakage problems FPOH the orthogonal PQLARliATlQH OR 

HHEH ITS POLtPIZATlCM PROPEHTitS AREN’T KNOUN HUH 
CFTTAIMTV, the PROCPAH PPE-SUFPrjSES THAT THE USER 
HILL ATTEMPT ANYONE OF THE FIFTEEN MEASUREMENTS AS 
CESCfilfiEO IN THE HEFORI QY J.P. CLAA3SEH EHIITLEO 

(CR I. THE USER HAS FCLP CHOICES OF AHTFNNA 

PATTEPNS AS EPECIFIEJ OY ITYPE =liZ<3 OR «t. THE 
types COrtXEjTOfO TO THE LAHCCA PATIEhNE OF TYPE 
H=-l/2.Qtl/2,l AS CESCPIQEO IN THE AAOtP HANPOODK. 

ChP. g . THE BFAH HIDTH CF THE CHOSEN PATTERN IS 
GOYERNOREO By ThL INPUT PARAHtltN KA. Th£ QEAMHIDTH 
IS PtUATLS to KA IN SIATEHCNT KC. OF TKf PROGRAM 

THE VIEM angle AT HPTCH THE USER WISHES To CONDUCT 
hIS STUDY IS SPECIFIEO IN TNOT. THE OUTCOME OF THE 
simulation is ilASEO OT) A SCATTERING CPARACTCR- 
ISTIC SIMILAR TO THAT OF THE SEA, RY RF-OLACIHG 
SUaPCUTIriE SICHAi tPE USER MAY ImtHCOOCE AnoTHER 
CHAFACTE"ISTIC. no IE that THE hOUTlNE MUST COMPUTE 
THE SCATTEPlNC COEE PfR UNIT STESADIAN. DIAS 
LEAKAGE DY THE OMThCGONAL POL AK IZA T ION IS INTROOOCED 
eV THE USE OF THE INPUT PARAHEIER ADIAS. 

THE PHASE OF This leakage is defined relative 

TO VERTISAL POLARIZATION ANC IS CONTROLLED WITH 
INPUT PARAMLTES DDIAS, TO CONCUCT PONTE-CARLO STUDIES 
CF THE OUTCCPE OF THE SCAT MEASUREMENTS WHEN SHALL UN- 
CERTAINTIES IN The AMPLITUDE ANC PHASE PROPERTIES OF 
THE antenna EKIST, INFUT PARMtTERS AHA* AND BMAX 
PAY OE SPECIFIEO TO OE OTHER Than ZERO. MMEH AHAX=0 
and DMAYsOt IT IS ASSUHEC THAT NO SLSH STUDY IS OE- 
EIREti. the CONTRAINTS OH TPF DIAS AND RANOOH PARA- 
PEtFRS Are OESCRIDIO in SUDROUTINE ANTENNA. 

NHEH BIASFS ARE NON-ZERO TriE MONTE CARLO 
STUDIES ARE CONDUCT FD WITH BIASES IhSERTFO. 

P0E = INTEGRAL(PATTERN*C0S(TMETAn**2 

PSCm = INTEGRALt IPATT£BN<*C0SITHETAI I**2*SIGHA«I»» 

Pnas = ICFAL oosep.vation matrix 

pact = actual observation PATRIX WITH PERTURBATIONS 
PINV = INVERSE OF THE NCRMAL EOUATIQNS FORMED 
FROM POBS 

SC = actual scattering COEFFICIENTS AT'tnOT 

TNOT = View angle 

ITYPE = ANTENNA TYPE 

KA * ANTENNA N0RPALI2E0 RADIUS 

AT = RELATIVE GAIN OF HORIZOMAL PATTERN 

CURING XMtSSIOH 

AR = RELATIVE gain Of HORIZONTAL PATTERN 


0000007D 

onaaooaD 

00000090 
DOOOOICO 
00000110 
Qooa 0120 
00000130 
OOaODlNO 
oooaoiGO 
OOOOOlGo 
00000170 
aoaoaioD 
oaoooloo 
ojooozaa 
ooocozia 
ooooozzo 

00000730 

DAOOOZAO 

00000200 

00000200 

00000270 

Doooozao 

OOOOQZ'ID 
000003(10 
00000310 
0000 (|320 
00000330 
000003NO 
O0OQD35O 
00000300 
00000370 
00000300 
O0G0O390 
OOODONfO 
OOOCQNID 
D0D0CN20 
OOOOOr.30 
OOOOONI.D 
OOODONGO 
OCOCGNGC 
OOOOONTO 
00000400 
00000490 
OQOCaSQO 
00000510 
00000520 
OODO0S3O 
00000540 
00000550 
OOOD05GO 
D00D0570 
00000500 


Figure D ,5a - FORTRAN LISTING FOR THE 
SCATTEROMETER SIMULATION PROGRAM - MAINLINE 


Taifc»» 01 0e-15“V5 23.069 


scar SIHLLftTION PPOGRAH 


LAQEL PAGE 


Z 



o 

OJ 


53 

54 

55 

56 
67 
58 
55 

to 

61 

62 

63 

64 

65 

66 

67 

68 
65 

70 

71 

72 
71 

74 

75 

76 

77 

78 
75 
80 
81 
82 

83 

84 
86 
at 

87 

aa 

ag 

go 

91 

52 

53 

54 

55 

56 

57 

58 

59 
100 
101 
102 

103 

104 


C 

c 

c 

c 

c 

c 

c 

c 


c 


c 


DUPitiG i<£C£r>Tian 

Q7 = PHASF PFLailVe TO V PATTEfiH OUHTHC T«AHS 
OB = PHflSF PEiaTIvE TO V PATTERN DURING RECP 
ADiaS = PATTL'm AMPLITUDE 01*$ 

PilIA5 n PATTFPN Pl(A5r Q1A5 
AHAK = HAX PATTffilJ PFRTUBBATIOH 
DHAX = MAX PHASE PERTUBOATION 

COHKOF) /ONE/ Htisli P0QSHT5t9»« POQSEClS.git PI«VI9i5Jt SCT3) » 
1 LABELUOI. AT(15>, ARUS). QT(15). GRllSI. SINTNt CDSTN. 
t XAi ITVPf 

EIHEH5IC6 THrTA{15a>. SNODIA). ANTt4). HI0TH(4>. PACTilStlOlt 
t P5CEt6.5»iPSCT«5I.CI9.F.,151 ,P8E 16 ,5 1 « PD 1 15> .PAT161 .Q(6T 
EXTEPNAL LAMBDA 
PEAl ka.lahdca,kao 
logical JUHPt,JunP?.ISING 
CATA STATEHCNTS 

CATA XAC.TIIOTO.I TVPEQ /2*-1.0«-l/ 

3ATA HICTH /0,«8, 1.02. 1.15. 1.27/ 

CATA ANT /-C. 5,1. Oil. 5, 2.0/ 

DATA SNCa/6.283,7.ai6,7.725.S.416/ 

CATA Pi/ 3, 14155265/, PELTA/4. 36332313E-0 3/. OEG/D. 0174532525/ 

CALL SETCIH(Pim/,5.91 


C 

C 

C 

C 

C 

C 

10 

L 

1000 


A 

2000 

L 


INPUT ANTENNA TVPE, NORHAUIZEC RADIUS. 
viEH angle, ahpl bias, phase bias, max AHPL 

FFPTIJflATIQN. HflX PPA5F PERTUeBAIIOH. 
sample size in KONTE carlo STUDY. 

PEAO (5. 1000.ENO = 230I IT YPE ,KA, TNDT ,AOI AS .BPNA5E, 

AHAX.FPMASG, JCASES 

fOPMAT(I2«6rt.5.I5) 

IFCITYPE .GT. 4 .OR. ITYPE .UT, II CALL ABORT I2HTVI 
IFl.NOT. IAkAx .LT. 0.0 .OR. ABIA5 .LT. 0.0 .OR, 
U6IAS.-AMAXI .GT. l.OII CO TO 20 
HPITE 14, 20001 ADIAS.AHAX 

FORHATtmi, IX, •PERTURBATIONS TOO LARGES5X. 
•AerAS=’F10.e,3X,*AMAX=*F10.Sl 


CO TO 10 


c 

C ROUND ANT LOOK ANGLE TO NEAREST HALF DEGREE. 

C 

20 JUHPl 3 KAO .EO. KA .AND. I TYPED ,EO* ITYPE 
TNOT » IFIyldNOr ♦ 0.251 * 2.01 / 2.0 
UUHP2 = TIIOTO .EQ. TNOT 
ISIHG 3 .FALSE. 

IFtJUMP2i GO TO 25 
JtJOTP n TNOT * DEC 
COSTH = COS I TIIOTR I 
CINTN ■ SIM I THP»'F< 1 
25 WRITE (6.3000) TNr , 


00000590 

DDOOOT-00 

00000610 

00000620 

08000630 

OOOnOF4Q 

0000DE5D 

00000860 

00080671) 

OODOOebO 
00000850 
00000700 
00 000710 
00000720 
00000730 
03000740 
00050750 
00000760 
00080770 
00000780 
00000750 
00000800 
ODOCOSlQ 
00000820 
OOOROil.lO 
O00OPA40 
00010R50 
DDOCCnGO 
00000870 
OODQOauO 
OOOGQ890 
00000500 
00000510 
00 00 0520 
00000530 
00000940 
00000950 
09000980 
09000970 
0ODO09H0 
00090990 
OODOIOCD 
00001010 
00001020 
00001030 
00001040 
O00D1050 
QD001060 
00801070 
00001080 
00DD105O 
00001100 
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01 02-iS-rS 23.CB9 


SCAT StHULATION PROGRAH 


LABEL MAIN 


C»A&E 


3 



105 

3000 

FORHATClHl//' SCAT STUOV FOR VIEM AKGLE OF*, 

00001110 


106 

t 

F5.1.* 0FG«£ES*//» 

00001120 


107 

c 


00001130 


lOA 

c 

COMPUTE THETA (HAK ) 

00001140 


ItB 

c 


OODO 1150 


110 


IFIJUHPU GO TO 30 

D00D1I60 


III 


GIHTM » SNCOlITTPEl / «A 

QD0Q1170 


liz 


COSTH » SQPTC 1,0 - SIHTH * SIHTH 1 

00001160 


113 


THAX 3 ATAIK SlHTM,COSTH l/OEG 

00001190 


IIA 


TMAX = IFIX{(THAXtQ,2GI * 2.0 1 / 2,0 

OOQO 1200 


115 

c 

• 

00001210 


lib 

c 

CCiPUTC A'lTESdA paraheters. 

00001220 


117 

c 


00081230 


114 


call SOlIOI LAHB[)A,ITrPE,KA,l.Q,COSTHiia,<i,S> 

00001240 


ll"» 


GAIN £ 2,0 / S 

- 00001260 


t?C 


reCTOF = GallieGAIH / t 4.0 * PI • KA »**2 

00001260 


121 


GDf) = ID.O • ALOGloe GAIN 1 

00001270 


122 


he«H = HIOTH(l7VPe>*PI/IKA*CEG) 

00001200 


12J 

30 

00109 = l!l.*AUJG10tAniAS*l.l)E-2s) 

00001290 


12A 


r»ANO = in.*ALOGHItAHAx+l.(lE-7Bt 

00001300 


125 


MOITE (6,4000) ANT(ITVPE) ♦«A,QEAM,GDe.CaiAS,ePMASE.CRAHD,RPHASE 

00001310 


12b 

4000 

FO»flAT(//,20X, •ANTENNA PAR4i*et£4S*/ 

0 0001 ’20 


127 

4 

IX, ‘ITPE'TX.'KA'.GX, •HI0TM*,5X,*CAIH*,6X, •cross*. 6X,*RELA*/ 

00001330 


124 

1 

20)f,* (0EGJ*,5X,* (OIJJ '.ex.'ICOl’.eX.’PI-ASE’.SX.'AKAX'.rX.'BHAX*/ 

00001340 


129 

1 

F4.1,1X,7FID .2//) 

00001350 


110 


PRIA? = rpMASF*OFC 

00001360 


131 


QHAX = ftPPASE’DEo 

00001370 


132 

c 


00001380 


133 

c 

CHECX MHETHER sake AnTEMHA and ANGLE 

00001390 


1 JA 

G 


00001400 


U5 


IFUumPI .and, JUHP2) GO TO 195 

00001410 


136 


KRITFtb.SOOO) 

00001420 


137 

500C 

FORMAT (//IX, ‘STRIP*, 2X,* theta *,33X, ■HEIGHTS*. 25X, 

00001430 


138 

t 

•PRECISICH*//) 

00001440 


139 

c 


00001450 


lUC 

c 

COMPUTE NUMBER OF SAMPLING ANNULI* 

00001460 


t«-l 

c 


0000 14?0 


1A2 

c 


00001400 


1A3 

c 

AJ RIGHT OF EORESIGHT. 

0000 1490 


lAA 

c 


00001500 


145 


IHCR = 2.0 * THAX 

OQOOlSlw' 


146 

c 


00001520 


147 

c 

B1 LEFT OF eORESIGHT. 

00001530 


148 

c 


00001540 


149 


IF 1 TKQT - THAX !35,40,40 

00001550 


150 

35 

JljCL ■ 2.0 * THOT 

000G156Q 

f' 

151 


GO TO 50 

00001570 


152 

40 

IHCL » INCR 

09001580 


153 

C 


00001593 


154 

c 

TOTAL NUMOER OF ANNULI* 

00001600 


155 

c 


00001610 


156 

50 

I TOTAL * IHCR + IHCL + 1 

00001620 
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iriitimmia 1 1 II iiinmWHtot 



T81«.ii 01 02-15-75 


23.0S9 


SCAT SIHULATtort PPOGRAH 


LAB£t HAIH 


PAGE 



» » 


<> 


157 

C 



00001630 

isa 

c 


COHPttTE HIQPOINTS CF SAKPlING STRIPS, 

000D1640 

159 

c 



00001650 

leo 



IP 1 I«CL *eO. 0 1 GO TO 70 

Q0001660 

161 



ITEHP = IhCL ♦ 1 

00001670 

162 



CO 60 la l.IKCl. 

DDODlOao 

163 



THETAtl) - TMOT - 1 ITEHP -11/ 2.0 

D0QQ1690 

16t. 


6(> 

CONTIKUE 

00001700 

165 


70 

DO ao la lilMCRtl 

00001710 

166 



THETAUNCL *11 a TROT + t I “ 1 1 / 2,0 

00001720 

167 


60 

CONTINUE 

00001730 

168 

c 


- 

0000171,0 

169 

c 


CLEAR ACCUMULATORS 

00001750 

170 

c 



00001760 

171 



DO S5 la 1,9 

00001770 

172 



FSCKII a 0.0 

00001760 

173 



PBIin a 0.0 

DDQ0179D 

17<. 



DO 65 Jai,6 

OOOOlflOO 

175 



pQEfjtnaQ.a 

0000 IHlO 

176 



FSCEtJ. 11=0.0 

00001820 

177 

85 

COHTIUUE 

0000 1030 

179 



CO 150 11 = IflTOTAL 

00001040 

179 



I = II 

00001050 

160 

c 



00001660 

161 

c 


A1 LIMITS ON COSdHETAl 

00001B70 

162 

c 



0000 laaQ 

163 



THETAR = TMETA(t>»0EG 

D00D1690 

164 



X2 a COSITHETAR-DElTAI 

0000 [900 

165 



IF (ThETAR .LT. 0,00011 X2 a 1,0 

00001910 

166 



XI = COS(^':£TAR*oeLTAl 

00001920 

167 

c 



OP001930 

186 

c 


01 LIMITS ON PHI. 

00001940 

189 

c 



00001950 

no 



CEKOH = SIH(TKETAR)»5TNTH 

00001960 

191 



IF (CENOH ,LT, 0.000015 CO TO 90 

00001970 

192 



CD3PH = ( COSTM - COStTHETAfiMCOSTH l/DENOH 

00001960 

193 



IF 1 COSPHI .GT, -l.D 5 GO TO 100 

0090 1990 

194 


90 

PHI = PI 

oaaosiioo 

195 



GO TO 110 

00002010 

196 


100 

FhI = ATAk 2( S0RT( ABSti,0-CaSPHI*COSPHI) 1 • 

00002020 

197 


X 

COSPHI 1 

0009 2C3D 

196 

c 



00002040 

199 

c 


Cl SET NO. OF IKTECRATION DCMAlNS 

ODDC205D 

200 

c 



00002060 

201 


110 

INCV a PHI / OEG * 1 

. 00002070 

^ 00C020R0 

202 



IFCIHCV ,GT. 31) INCY ■ 31 

203 

c 


• 

00002090 

204 

c 


C) INITIALIZE CONVERGENCE TESTING PARAMETERS. 

000O2100 

205 

c 



00002110 

206 



00 114 J a 1,6 

OOOQ2120 

207 



OIJI a Q.D 

00002130 

200 

114 

CONTINUE 

00002140 
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01 

03-15-75 

29.C49 scat sIHOlATIOH PPOGRAI* 


label 

309 

C 



00003150 

310 

C 

El ItlTEGRATE WITH PRECISIOH J, 


00003160 

211 

c 



0D0O2170 

313 ■ 


DO 120 J ■ 4iB 


00007100 

2H 


Jj a J 


O00C3190 

314 

c 

> 


OOOD33DO 

316 

c 

F) COMPUTE MEt^HT CN ITH AKNULUS* 


00003310 

316 

c 



00D0333D 

317 


CALL FXOPTI67,l,t,01 


0D0D233Q 

31Q 


CALL OIKTEGI P AT, X3 iXl ,PH I.-PHl tI»CV* JJ« JJ ) 


0000334D 

319 


CALL F)tOPrC67tl,0,C> 


00003350 

330 

c 



00003260 

331 

c 

G) TEST FOR convergence. 


00002270 

333 

c 

* 


00063300 

?33 


CO 116 X * 1,6 


00003390 

334 


IPIpATCKI .LT. l.OE-281 60 TO 116 


00003300 

335 


IF < A[jS(PAT(KI-QlKn / PATCKJ ,LT. l.E-5 1 GO TO 

116 

00002310 

336 


GO TO 119 


00007330 

337 

116 

CONTINUE 


D00C2330 

336 


GO TO 130 


00003340 

339 

114 

DO 120 K 3 1,6 


00003350 

310 


next = PATtKJ 


00003360 

3 n 

120 

CONTItlUE 


00003370 

333 

C 



00003390 

331 

C 

M convergence aohitteo 


00003390 

3 54 

C 



00003400 

335 

110 

hPiTE [6,60Q0> 1 1 .TPETA tl 1 ,PAT , JJ) 


OE003410 

316 

bOOC 

FComatI I4,F5.1i6E11,3»I5I 


00(103430 

237 

c 

II 5YMIHESI7F parameters TO COMPUTE EXACT RETURN 

AND 

00003430 

316 

c 

CeSFPVATIOH MATRIX 


00003440 

3 39 


CO 150 Jal,9 


O0QD3460 

340 


Jjaj 


000034GO 

341 


5CATC = StGMAUJiTPETAlI) > 


00003470 

343 


CO 140 K?l,6 


00003400 

343 


FSCMK,J)=P5CE1K,J) «PAT (K 1 « SC ATC 


00 003490 

344 


r0£(X,JI =POE lK,J>+PATtKI 


00003500 

346 

140 

CONTINUE 


00002510 

346 


PCCHJI = PSCHJl VPAT (6>*SCATC 


00003520 

347 


P6HJI = PBl IJJ+PATI61 


00003530 

349 

150 

COtlTinUE 


00003640 

349 

c 

HSITE16, 51001 FACTOR, Pal,PSCl 


D0B03550 

360 

C 6100 

FOPHATUH,IOE12.4/9E12,4//1 


00003560 

351 

c 

call HAT0UT«PDF,6,9,6,9rl«PeE,3HRAHI 


00093570 

353 

c 

CALL HATCUT IP5CE,6, 9,6,9, 4HPSCE,3HRAH» 


00002540 

361 

c 

1) COMPUTE PEFERENCE SCAT. CCeF, 


00003590 

254 


30 HO I = 1,9 


00003600 

356 


II = 1 


00003610 

256 


SCtlJ = SIGHAJII.THCTl 


00002620 

357 

c 

JJ SCALE FOB antenna GAIK , 


00003630 

354 


FQIUI a POI tll*FACTOR 


00002640 

259 


P5DH1J » PSCni)*f ACTOR 


00003650 

260 


DO lED J » 1,6 


00002660 
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T014I, Ot 02-15-75 23,D«5 SCAT SIHULATIOM PROGRAJ1 LABEL MAIM 


261 


PSCECJfll ■> P5CE( J. :> “PACTOB 


O0OD2670 

262 


P5ElJiI» = PBEIJ,I)*FACT0R 


OOOD2GBO 

263 

160 

C0»TIHUE 


00002690 

SGA 

c 

HPlT£t6*510DI FACTOfitPUI.PSCl 


001IG2700 

265 

c 

HALL KATOUTtFflr *6*9,6. 9i3HPCE,3hSCA1 ' . 


00002710 

266 

c 

CALL HATCUr (PSCE*6,9,6,9,4HPSCE*9HSCA» 


00002720 

267 

c 

EO^H IDEAL PENCIL QEAH HEIGHTS AND POKER HATRi;( 


D000273D 

266 


WRITet6*7Q00> 


000027411 

265 

7000 

FOPMATIlHl, 'lOEAL ANTEriHA HEIGHTS AND POWER HATRIX* > 


00002750 

270 


call AHTFRHAtO. 0,0.0 ,0.0, 0.0) 


00002760 

271 


CALL IDFAL(PCI,PO0Sri 


00DD2770 

272 


CALL IOC AUP9CI, PACT) 

• 

000O27B0 

273 

c 

APPFHO POWER VECTOR TO PACT FOR OlSPLAY 


00002790 

27A 


C’O 170 I = 1,15 • • 


Q00G2900 

275 


PACTtr.lC) = HID 


00002910 

276 

170 

COMTIliUE 


00002920 

277 


CALL HATOWT(PD031,lS,9,15,9,ShDEUTA,6HHEICHTI 


00002930 

27B 


call hAtCUT {PACT,15,10,19,10,5HPOHER,6HHATRIXI 


QDD021I40 

279 

c 

FORK Exact pencil heights and poher hatrix 


00302050 

2B0 


HPITf «6,7S00) 

, 

00002060 

291 

7500 

FOPHAT JlHl, 'EXACT AHTEHRA HEIGHTS AhO POHER RETURNS* 8 

, 

03002S70 

2B2 

c 

HPITF (6,71QGI (AT II) ,QTI 11 . ARtll ,BR(1I ,1=1,15) 


OOD02UBQ 

283 

C 7100 

FOPKATI//15I4E12.A/I ) 


QD002'’9D 

296 


CALL COFFtC) 


00002900 

295 

c 

WRITF{6,61QCI f ( 1C 1 1 , J ,K) , J=i ,61 , 1=1.9) ,K=1,15) 


00002910 

296 

C 6100 

F0P«AT(//9I6E12.4/)I 


030D2920 

287 


CALL FXACTIP9E,PaOSE,C,H) 


00002930 

296 


CALL HAlCUTCP0aSF.l5,9,15,9,5hDELTA,61-HEIGHT) 


00002940 

289 

c 

prepare normal eqks. 


00002950 

296 


DO IBO 1 = 1,9 


00002960 

291 


DO 190 J = l,I 


00002970 

292 


PINVII.J) = 0.0 


C00029BD 

293 


no 175 K = 1,16 


00002990 

294 


PIMVCI.J) = PINUtI, J) ♦P0BSE«,I)*P0BSE<«,J) 


DDCO3D00 

295 

175 

COHTIliUE 


QC0Q3C1Q 

296 


piHvcj.n = PiHVd, J) 


DD003'>20 

297 

180 

continue 


00003030 

299 

c 

COMPUTE AND DISPLAY INVERSION HATRIX 


00003040 

299 


CALL HEMlrJV(PIHV,9,PBI,$260) 


OODDItlSB 

300 

c 

FORM EXACT POHER hatrix 


00003060 

301 

105 

CALL EXACT CPSCe, PACT, C,H) 


00003070 

302 

C 

APPEND POHER VECTOR TO PACT FOR DISPLAY 


OJ0G3DB0 

303 


CO 19G I = i,lS 


00003090 

304 


PACm.iO) = Hit) 


000031 DO 

305 

190 

comtihue 


00003110 

306 


CALL HATCUTtPACT,15,10,15,ia,6HPOHCR,6HKATRIX) 


00003120 

307 

C 

COMPUTE SCATTERING COEFFICIENTS 


00003130 

300 

c 

A) DV THE DIFFERENCE HETHOD 


00003140 

309 


HRITEI6,e500) 


00003150 

310 

6500 

FaRHATIlMl) 


00003160 

311 


call differ 


00003170 

312 


CALL CIFshOH 


DD0Q31S0 
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02-15-75 

23.CA9 SCAT SIMULATION PROGRAH 

label 

313 

C 

El ev the matrix methoo 

00003190 

3lti 


ir<ISlriGt GO TO 195 

00003200 

315 


CALL matrix 

00003210 

316 


CALL HATSHOh 

00003220 

3ir 

195 

TIIOTO = tNOT 

00003230 

3lfl 


KAO s KA 

00001240 

315 


ITYPEO = ITYPE 

00003250 


C 

CHECK IP CASE UITM DIAS IS DESIRED 

00003260 

12l 


IFf.hOT. IA9IAS .GT. 0 . 01 ) Go To 210 

D0ni!3?70 

322 


HR1TE(6,6D0C> CQIAS.nPHASE 

0000 3250 

323 

oooo 

FORMAT IJH,* AUTERNA HITM DIASES ONLVV/ 

00003250 

321. 

. t 

lX.'AHPL nIA5=*.F6, 1. • OS'. SX, 'Phase BIAS=.'»F6.2.* D£G*//B 

, 00003300 

325 


call ANTENNAtO.OiAOIAS.O.O.QQIAS) 

00003310 

326 


call CQEIFICI 

00003320 

322 


Call cxact(psce,pact.c.«i 

00003330 

32H 

c 

append POHER vector to pact for DISPLAV 

■ OPOOiSAO 

329 


to 200 I = 1.15 

00003350 

330 


PACTtl.lOl = mil 

00003360 

''1 

200 

CONTIM'JF 

O00C3370 

’,2 


CALL MATQUT (PACT. 15 . 1 D 1 15 . 10. SHPOuE.9. 6MHATRIXI 

03003300 

333 

c 

Al BY THE OIFFEREnCe METHOD 

000C339D 

331. 


CALL differ 

00003400 

3 35 


CALL CIFOHOM 

0000 3410 

336 

c 

6) ev The hatr'k hethoo 

00003420 

337 


IFdSiNGl Go TO 210 

00003430 

3 3fl 


CALL MATRIX 

00003440 

339 


CALL MATGHOM 

000O3450 

340 

c 

CH5CK IF MONTE CARLO STUDY DESIr^o 

QUOC 3450 

-141 

210 

IF (.hot. (AHAX .GT. 0.01) GO TO 10 

00003470 

342 


HRITR I5,900QI CDIAS .CPhASE.CRAhO. RPMASE 

00003480 

343 

9000 

FORMAT (IPl//* MONTE CARLO STUDY*// 

0000 1490 

344 

1 

IX.'AhPU DIA5=’ .F5.1,' DB'.5X. ‘PHASE EIAS=*' .F6. 1. ' DEG', 

ODQD 35C0 

145 

5 

MAACQH aMPL = * .F6.1.* Oe'.SX.'RAMOOH PHASE»'.5X. 

003035 10 

346 

t 

F6.1,* t,rGV/l 

0000 3520 

347 

C 

perform MONTE CAR! s, . ' V 

0000353Q 

340 


ro 220 I-l.ICASES 

OCOJ 3540 

349 


call ANTFMNA (AMAX.ABIAS.SHL, -.IAS) 

00003650 

350 


call COEFCC) 

00003550 

351 


CALL EXACT(PSCE.PACT.C.H) 

0Q0C3570 

35? 


CALL 1 IFFER 

00003560 

353 


IFtlClnC) GC TO 220 

00003550 

354 


call matrix 

00003503 

355 

220 

COIITIliUE 

00003510 

355 

C 

SHOW RESULTS OF STUDY 

00003520 

357 


call cifshoh 

00003530 

350 


IFtlSING) GO TO 10 

00003540 

359 


CALL MAlSHOM 

DOllCaEOO 

360 


r.D TO 10 

00003660 

351 

230 

STOP 

D00S357Q 

362 

240 

HRITF(6t9500l 

00003660 

363 

9500 

FORHATI/lX.lO I1H*I , 'MATRIX SIHGULAR'.IQ (lM«l//> 

00003550 

354 


ISInG a .TRUE. 

00003700 

355 


GO TO 155 

00003710 

366 


END 

00003720 


PACE ? 


Figure D,5g - FORTRAN LISTING FOR THE 
SCATTEROMETER SIMULATION PROGRAM - MAINLINE 


01 02-15-75 23.052 A SUSPOUTInE AHTEtJNA 


LABEL AnT2»H 


1 

2 

3 

k 

c 

6 

7 

a 

5 

ID 

U 

12 

13 

1«* 

15 

16 
17 
la 

19 

20 
21 
22 
23 
29 

25 

26 
27 
20 

— 29 

O 30 

vO 31 

32 

33 
39 

35 

36 

37 
3a 
39 

90 

91 

92 

93 
99 

95 

96 

97 
90 
99 
5D 

51 

52 


cahtehwa suopqutihe antemma 

SUBUOUTINE ANT£HNA1AHAX,ABIAS,BHAX,BBIAS» 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS SODPOUTIME SVHTHESIZES THE RELATIVE AMPLITUDES 
AMO PHASES OF THE TMO ORTMOGOMAILV F0LAR,I2E0 ANTEHHA 
PATTERNS DURING TRAMS HISS SION T ANO RECEPTION R FOR 
THE FIFTEEN STANOAFD HEASUREHENTS. THE INPUT AR- 
CUPFNTS PERHIT THE PRECISE ANTENNA REOUIBEHEKTS FOR 
EACH OF the FIFTEEN HEASUREHENT CONOITIOHS TO BE 
PERTURBED with BIASES EItMEK FIxEO OR RAmOOH OR BOTH 

AhAXsHAXIMUH RAHOOH PEtiTURBATICN INTROOUCEO IHTO THE” 
CORZONTAlLV POLARIZED PATTERN 

A6IAS=BIAS INTBOOLICED INTO THE HORIZONTALLY POLA- 
RIZED pattern 

CHAYsHAXIHUH RAnJOK perturbation 1 TROOUCED INTO THE 
RELATIVE phase BETWEEN ORTHOGCNALLV POLARIZED 
PATTEpmS 

eeifiS=EiAs ihtroouceo into the relative phase 

BETHEEN ORTHOGONALLY POLARIZED PATTERN 

IF THE INPUT ARGUMENTS ARE SET TO ZERO, PRECISE 
AHTEHNA REQUIREMENTS ARE ESTABLISHEC IN ThE OUTPUT 
VECTORS: 

AT<R!=RELATIV£ amplitude of TPE horizontally POLA- 
RIZED PATTtHH DURING TRAIISH lES ICN (RECEPTION) 

CT(P)=PELATIV£ PHASE EETREEH OWTHOZONAL POLARIZATION 
CURING TPAIiSHtSSION (RECEPTICKl 

OTHfRMiSE, PFHrUPgATIOHS ARE INTRODUCED IN ACCORD 
NITH ALGORITHMS BELOH. BIASES ARs EFFECTIVE ONLY 
HHFN polarized TRANSMISSIONS OR RECEPTIONS ARE HADE 
IN IHE FIFTEEN H EASUREHEH TS. IT IS ASSUHEO THAT 
LEAKAGE OR DE-POLARIZATION IS ThE CAUSE OF THE 
elASES. the user hlST AQ3ERVE THAT 

11AHAX .GT. 0. 

21AOIAS .GT, D. 

3) AHAX'i'ABlAS .LE. 1.0 
91-PI .GT. BDIAS .LE. PI 
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POT (1,1.1 =SK»C05(dR( It *BT ( I J >*PIM(AI 
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0l926«73 16,309 SUSDOUTINE SOLID 

CSOLID SUOFfOUTlN^ SOLIO 

SUBHOUTIije SOLID «FCT.K*,WNA,X8P,KlP,N>t*iNHP,S|JH) 

C 

C THIS HlUTINe EKPLOVS A OUASStAN UEOENDUe OUADnA?UHE 

C (INTEOHatIOH) pPOCEDUREi INTEGRATION OVER }i SE0HENT3 

e ARE PEriFoRHED AFTER TRANSLATION TO l-ltllXtcilti) , 

e KEPiXiri a UPPER AND LOUER UH|TS ON X 

C fiNX « "EtJNENTS IN X 

C NNP 9 lEtiHEE OP PRECISION IN K 

C 

OIHGNSlO J SAHPLEtOtB|,OOEr(8,8),5x(01,CI8} 

data ((SARPLEU I I J*w6)* 1 92(6 1/-0. 577350269,0, 577350269» 

8 0i0,0,0,0i0,0i0(0,0i0,0( 4 

8 "0 .774596669, 0, D, 0, 774596669, 0,0.O« 0,0. 0,0 1 0, 0,0* 

6 •0,06u363l2, -0.339931044,0,339981044,0,861136312, 

8 0,0, 0,0, 0,0, 0,0, 

8 -O, 905179846, -j, 53846931010. 0,0, 530469310, 0.906179846, 

A 0,0i0,1,Q,Q, 

A »Q,932.16 75i4,-.rj,65l2o9306,-'Oi2386i9l86, 0,230519 166,0,661209306, 

A 0.932449514,0, 1,0,0, 

A »0, 949 107912, •), 741531 t»5, -0,4 05045151, 0,0. 0.405845151, 

A 0,741531185 ,0.949107912,0,0, 

,1 -fl, 96 0289856, «n. 795666477, -0,525532410, *0,183434642, 0,183434642, 
A 0,5255.32410, 0,796665477,0; 960209056/ 

data uoErii. )),j>i,Bi,ip2ja) /i;o. i,o, o,d,o,q,o.Oi o.o,o.o, 

A 0,0, 

A 0,555555556 ,0,880809089 ,0,555555556 , 0 , 0, 0 i 0 , 0.0, 0,0, 0 ,0 , 

A 0.347854851,0,652145155 ,0,6S2l45l55 , 0 ,347854851, 0,0, 0 , 0 i 0. 0, 

A 0,8, 

A 0,236926,385,0 1 470628 670,0; 566688009, 0,4786286701 
A 0,236926835,0, ,0,0, o'.Di 

A 0,i7i324<92, 0.360761573,0 .467913935, 0,467913933, 0,360761573, 

A 0fl7l334492,o,a,o.0, 

A 0,l29434966,0.2797oS39l,o;38i83o051,0,4l79S9l8j, 

A O; 301030051, 0,279705391, o', 129404966. 0,0, 

A 0 1 101228536 1 0, 222381034,0,313706646, 0,362633783 i0,362683783t 

A 0;3137q6646, 0,222881034,0,101220536/ • 

C 

c clear SUHhInQ variable 

e 

SUH 4 OiO 

c 

c re«ass|On Input arguments 

c 

rN « RUA 
H » HA 
NX ■ NNX 
NP a NMP 
C 

c cohpute length or cell sides 

c 

DEuX 1 tX2P-XlPJ/FL0AT|MXJ 


label solid 

00006800 

D000689Q 

00006900 

0000691Q 

00006920 

Q00069JQ 

00006940 

00006950 

00006960 

00005970 

00006988 

00006990 
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OOOO7Q10 

00007020 

00007030 
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00007068 

00007678 

00007080 

00007090 

00007100 

OOOOTllO 

00007120 

00007130 

00007140 

00007150 

00007160 

00007170 

00807180 

OD007l9o 

00007200 

ODOQ7210 

00007220 

00007238 

00007740 

00007250' 

00007260 

00007270 

00007200 

00007290 

00007300 

00007310 

00007328 

Q0U07330 

00007340 

00007350 

00007360 

00007370 

00007380 

Q0007390 


Figure D. 15a - FORTRAN LISTING FOR THE 
SCATTEROMETER SIMULATION PROGRAM - SUBROUTINE SOLID 
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5UBR00Y1KE SOLID 


LABEL SOLID PASS 


3 


93 


HDELX a OELXaO,9 

00007400 

54 

C 


0(3007410 

95 

C 

rORH SAKPLS rAOTOR FOR X 

00007420 

96 

c 


00007430 

97 


DO 2(3 t a liHP 

OOQQ744Q 

98 


SX{13 a sample {NP f I JtHCELX 

00007490 

39 


C<I> = COEr(NP.t) 

00007460 

6D 

30 

COHTtNUE 

00007470 


c 


DQQ07400 

62 

c 

integrate in strips or oelx 

00007490 

63 

c 


(30007500 

64 


XH a XlP * MSEtX 

• QOQ07510 

69 


DO 40 n l,tiX 

00007520 

66 


XM a XH * DElX 

Q0QQ7530 

67 

c 


09007540 

66 

c 

transfohh to call 

0000755(3 

69 

c 


00007560 

70 


DO 30 1 3 tiNP 

00007570 

71 


X a SX(1)*XK 

Q000756Q 

72 


XHSIN a RHaSoRTIl.a-X*X7 

00007590 

73 

c 


00007600 

74 

e 

PORH PARTIAL SUH3 

00007610 

73 

c 


0(3007620 

76 


SUH a SUK*rCT(M,UNS1N}aC(l) 

00007630 

77 

30 

coNTiwje 

(30007640 

78 

40 

CONTlNjE 

00007650 

79 


SUH a Qi5»DELX«SUH 

00007660 

8(3 


RETURN 

00007670 

Bl 


END 

00007680 


00 

JSk. 


Figure D .15b •*“ FORTRAN LISTING FOR THE 
SCATTEROMETER SIMULAHON PROGRAM - SUBROUTINE SOLID 
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1 

CI.AK80A SUBI^OUTINE I,AHB0* 


00007600 

2 


REAU rUNGTlON 


Q0007700 

3‘ 

c 



00007710 

A 

c 

H IS A7 IKTERGER WITH VALUS le2(3| 

OR ^ 

00007720 

3 

e 

DEPENDItJQ ON antenna TVPEi 


00007730 

6 

c 



00007740 

7 

c 



00007790 

9 


GO TO ( i00i200,300i4QO 


00007760 

9 

c 



00007770 

10 

c 

H 4 1 ItEANS antenna TYPE ■ -1/2, 


0001)7700 

11 

c 


* 

DQ007790 

12 

100 

ir ( ADS(U) ,UT. l.E-27 1 GO TO 600 


00007000 

13 


LANBDA a ( SlNtU) / U iasE 


00007810 

H 


RETURN 


00007620 

15 

c 



00007030 

16 

c 

H 8 3 (ieans amtEnna type 8 li 


0000)040 

17 

c 



00007860 

18 

200 

U*hBDA * 12,0 » BJiXOXfOJ >**2 


□0007060 

19 


return 

. 

00037070 

20 

c 


. 

00007880 

21 

c 

H <• 3 MEANS antenna TYPE f 3/8, 


00007890 

22 

e 



00007900 

23 

3Q0 

IF t AaS{0) ;ut; 1,E-2T j go to 600 


OOOQ7910 

2* 


LANBDA 8 ( 3,0 / ( u • u > * (( SINIOJ/ U 1 

1 « GOSIU} 

00007920 

25 


RETURN 


00007930 

2b 

C 



OQ007940 

27 

c 

N 8 4 ‘IEANS antenna TYPE s 2t 


00007960 

29 

c 



00007960 

29 

400 

IF ( AQS<U> ,LT, i,E«27 > GO TO 600 


QDDQ7970 

30 


LAHBPA a { S,0/(U*U) • C 2,0 4 BJIXOX(U) <• 


00007980 

31 


RETURN 


D0007990 

32 

600 

lakbda H 1,0 


00006000 

33 


RETURN 


ooDoeoio 

34 


end 


ooooaoao 


Figure DJ6- FORTRAN LISTING FOR THE 
5CATTEROMETER SIMULATION PROGRAM - SUBROUTINE LAMBDA 
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CBJ2ER0 soaHouriNE Bjzeno 

heal function trJZEROCX) 

C t 

c cohputes the bessei. function of Integer qrdqr zero, 

C uses a polynomial HETHOO, 

c 


RFAL X.T1.T2 
C 

C CHECK TO SEE HNiCH APPROXIHaT tON IS NEEDSO, 

C 0<x<3 

IF ( X ,QT, 3|0 1 QO TO 100 
C 

n B (333333333 • X 

n ■ Tt»Ti 

EJZERO 4 (t((f ,0002100 • Tl > ,Q0a9444)*Tl • |0444<I79)4T1 
6 -'.3163366>*T1 • 1,2656200>*T1 ■ 2; 2499997 »»Tt • 1(0 
RETURN 


C 

c X > 3,0 

c 

100 71 » 3,0 / X 

BJZERO a t(C({,00014476«Tt 00072005 >*Tl * . D0137237JOT1 . 

i ,t>0O09gi2)»Tl ■(, 00552740) aTl - , 0000007?)«T1 • ,79788456 
C 

T2 » ttm(iOO0l35SS*Tl-,0Ofl29333)aTl-,0OO54lg3}aTia,00263S73> 
& aTi-,OOOo3954)*Tl», 00003954 >»Tl*,Q 4 l 66397 )»Tt-, 7 a 339816 *X 
BJZERO 9 OJZERO » CD 5 (T 2 ) / SQRTIX} 

RETURN 

END 
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Figure D.17 - FORTRAN LISTING FOR THE 
SCATTEROMETER SIMULATION PROGRAM — SUBROUTINE BJZERO 
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19 
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17 
10 
19 
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23 


CaJlHOV SUBROUTifjE BjiuOl; 

runCTjoN ejixoKOo 

)r<2 ,ar. 3,o> go to ib 

c 

c 0 |GE| K .UE. 3 

r 

1 B ,33333333303 
¥ 0 ¥“V 

RJlKO¥ c <(<({, B(}OOUCi20VoO,OCia3176l}OVa,OOfl4g3l9SoVo 
6 ,O395‘l289)B¥*,21(f93S73J«¥o,9O249905K>T6i9 
ftfTURB 
C 

C X ,GT, 3 

C 

ID V D 3.3/S 

RjjXOl B Ut£tB,aOB2QD33oVo,(!OH36S3)OVo,OQa495jj}*V<, 

fl , 0 0017105 )•¥«, 0165^7667 (ov., 0 3000156 ►or*, 79703490 
BJlXQX n BJlXQS/()!oS3R7 (KU 

Z a («<<£■•, 0B029l66o¥i;00079a24JoV4,0a074348)sra 
S ,0 0637079 }«¥* I nB0O56501®¥o. 12499612)0 7,2,356194390); 
BJlXOX 5 BJlSOXoCOStZ) 

RETURN 

END 
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Figure D.18- FORTRAN LISTENG FOR THE 
SCATTEROMETER SIMULATION PROGRAM — SUBROUTINE 8J1XOX 
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Figure D,19a ~ SAMPLF i^UTPUT FOR 
SCATTEROMETER SIMULATION PROGRAM 
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Figure D.l 9b - SAMPLE OUTPUT FOR 
SCATTEROMETER SIMULATION PROGRAM 
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0.1R73E-01. 

O.Zt72E-Ol 

0. 

0. 

0. 43242-01 

0. 

0.112DE-03 

0 . 



q 

Q. 2l7aF'01 

0. l«7’r-ij(. 

0.2172F-01 

0. 

0. 

-0.4324E-01 

-0.06112-10 

-0.1120E-DJ 

-0.7435E-13 



10 


0.lli73F-0£. 

0,2l72t-Dl 

0. 

0 . 

0.4269E-10 

0.4339E-01 

O.llOSt-12 

0.3747E-04 



11 

0 , 2 i?at*oi 

0.I973E-C'. 

0.2L72F-IJ1 

0 . 

0 . 

Q.43L2E-10 

-0.4339E-01 

0.11172-12 

-0. 37472-04 



12 

O.H71f-Ot. 

D.2170E-01 

0.2172F-01 

0. 

0. 

0.1120E-03 

C . 

0.4324E-01 

0. 



13 

0 • 1 S 7 3£* O't 

0.2170E-01 

0.7172F-01 

0. 

0. 

-O.U2DE-03 

-0 .7435E-13 

-0 .4324E-C1 

-O.D&llE-lO 




O.IH/IE-OU 

D-2170E-01 

D.2172E-01 

0. 

0. 

0.1105E-12 

5 .3747E-04 

0 .4265E.10 

0 .4339E-01 



15 

Q.1073E*ai* 

D.2170E-01 

0 .2172E-(U 

0. 

0. 

0.1117E-12 

-0.3747E-04 

D.4312E-10 

-0.43I9E-01 






POhEH HATRIX 








NEaS/CCEF vv 

MB 

VH 

VVHHR 

UVMHl 

UUVHR 

l/VVMI 

hvmhr 

HVHHI 

PO«ER 


1 

0.«31££ 00 

0.2113E-05 

0 .OG77F-05 

0.143fC-02 

0. 

0. 

0. 

0. 

0. 

Q.8321E 

00 

2 

Q.2UJE-0S 

0.P316E to 

0.96772-05 

0.1436E-02 

0. 

0. 

0. 

0. 

a. 

O.0331E 

00 

3 

0,7l7iJt-03 

0.7179C-Q3 

0.2513E-32 

-0.1436E-02 

0. 

0. 

0. 

0. 

0. 

0. ’4132 -02 

J* 

0.20 'HE 00 

a.2C03E 00 

-0 .43662-13 

0.4165E 00 

o.ieoaE-13 

0.4223E-12 

0. 

0.4223E-12 

0. 

0.0331E 

DO 

5 

0.2075E 00 

0.2076E DO 

0.06772-05 

-0.4151E 00 

-0.3746E-11 

0. 

-0.423OE-12 

0. 

-0.4230E-IH 

0.129 5£. 

-04 

6 

O.ZOT-jE 00 

D.?07«jE 00 

0.4 TT‘lf-(J5 

Q.7178E-03 

D.ieaaE-Q2 

0. 

0, 

0. 

Q. 

0 . 4 1 8 4 F 

00 

7 

0.2079c 00 

O.FST'F OD 

0.43362-05 

0.7l7aE-03 

-o.ieooE-oa 

0. 

0, 

0 . 

0. 

0 . 4142E 

O'! 

9 

0,t.lE2E 00 

0.3599E-03 

0.1261F-02 

0. 

0. 

0 .4245E-03 

D . 

0 .ia99£-05 

0. 

0.4152C 

00 

5 

Q.hlf. 2f 00 

D,359';r-0J 

0. 12612-02 

0. 

0. 

-0. 4245E-03 

-D,a453E-12 

-0.10992-05 

-Q.7299E-15 

0.4i74r 

00 

10 

O.aitOE DO 

tJ.J599t-03 

0. 1261F-02 

0. 

0. 

0.4150E-12 

0.4259E-03 

D.1D85E-14 

0.367dE-06 

0.4 18 22 

00 

11 

D,<.lf)2E 00 

0.I59SE-03 

O.1261E-02 

0. 

0. 

0.4233E-12 

-0.4259E-03 

0.10962-14 

-0 .367QE-06 

0. 41742 

00 

12 

0. 1599E-03 

0.41622 00 

0.12612-02 

0 . 

0. 

0 .1D99E-05 

0. 

0.4245E-03 

0. 

0. 4102E 

03 

13 

0.3S3SE-03 

0.4152E 00 

0.1261E-52 

0. 

0. 

-0.1099E-05 

-0 .7299E-15 

-0.4245E-03 

-0.S453E-12 

0.4174E 

00 

li> 

a .15‘j9E-03 

0.4162E 00 

0.12612-02 

0. 

0. 

0.10S5E-14 

0 .367BE-06 

0. 41902-12 

0.42592-03 

0.4182E 

00 

IS 

0.3S99E-03 

a,<>l62E 00 

0 .1261E-02 

0. 

0. 

Q.1096E-14 

-0 .3670E-06 

0.4233E-12 

-0. 42592-03 

0.4174E 

00 


Figure D.-19c — SAMPLE OUTPUT FOR 
SCATTEROMETER SIMULATION PROGRAM 


STATISTICS Ton THE OIFPESEHCc H£THf>3 BASET, OH 


t OBSERViltCMS. 


SCAT COEF 

OAIOE 

pfah 

PHS 

KHEAN 

'-SRI'S 

liV 

O.lOSE 02 

0 .276E-02 

a.7a2F-05 

0.D14 

o; 

t-M 

A. 1‘iit 02 

0 .27EC-02 

C . 762E-05 

0.01<i 

0. 

WH 

C.sntc-fll 

-0. P07F-D3 

a 4 A?<AF'07 

-0.357 

0. 

VUHHR 

D,-32£ J2 

0.23BiE-fl2 

0,e&A£-05 

0.C12 

0, 

WHKl 

0 <OT2‘^-01 

•B.2i6£-03 

0-*i6/3E-07 

-0.250 

0. 

SiUVMP 

0 .SB2E’02 

-5.173E-fltf 

0.2'9‘J£-D9 

-0.176 

1)4 

UWVHt 

a.9f2E‘02 

* 0 , 3GB£-af. 

Q.956F-13 

-0.003 

0. 


0,B02E-a2 

-[»4 17JE-0". 

0.2B9£-O9 

-C.l ?6 

0. 

fVMHl 

OflBZf-O? 

-0.30'JF-a6 

D.asaF-ia 

-0 .003 

0. 


STATISTICS FO« Thf MATRIX ('eTMOI) OASE& CN t OBSERVA TICKS 


SCAT COEF 

UALUE 

pfah 

H li 

SHE an 

*RFS 

UV 

0.IB2E 02 

0 .256F-02 

O+e5fe£-05 

0.013 

0. 

PH 

0.157E 02 

O.266L-02 

0-fc57£-Q5 

a. 013 

0. 

VH 

0.5fllE-Ql 

-0.7O6F-05 

O.F18E-10 

'0,015 

0. 

VVHH» 

D.l'iJC 02 

Q.266E-D2 

C.£5h£-1)6 

0.013 

0. 

yVHHl 

Q.a72£-01 

-0.£62E-05 

0 .5jae-O0 

-0.076 

0. 

VVt/HP 

o.'iBtt-az 

-O.33lie-06 

0.J05E-12 

-0.003 

0. 

VUVHl 

0.g«2£-02 

-0.2S4E-06 

0.B65E-13 

-0,003 

0. 

pvkhR 

Q.°e2E-Q2 

-0.359E-C6 

0.12qE-12 

-0.005 

0. 

MVMMi 

0.9S2E-02 

-0.317E-06 

O.lOtE-12 

-0 .003 

0. 


Figure Dol9d - SAMPLE OUTPUT FOR 
SCATTEROMETER SIMULATION PROGRAM 


monte CA»lO 5T<jDV 

tPPL Bias= -i^O.O Da PHASE BIAS* 0 * DEC f?fttiDOH ftHPL” -iO.D BS ft&HOCH P(lflSE« 


statistics fOfl TME DIFFEBEHCE METHOD DASED 0« 150 OBsefitf ATIOHS 


SCAT CCEF 

VALCE 

(-FAH 

FHS 

2MEAN 

i:rks 

Vlf 

0*19ZE 0? 

0,317E-02 

0.113E-05 

0.017 

O.DOO 

HH 

. Q.1S2E 02 

D.310E-02 

0.119£-05 

0.017 

0.000 

WH 

0.*eiE-01 

0.7B5E-02 

D> 3O6E*0s 

13.530 

0.005 

VVHMR 

O.ISZE 02 

“D.e75F-0t 

D.55.IE-01 

-0.352 

0.2UR 

VVHHI 

0.fl72E-Ql 

-0.235F-02 

0, ULai-OX 

-2.701 

51,910 

WtfVKO 

0.Sft?E-02 

0.17tE 00 

0.362E 00 

1739. 70A 

3691.119 


0.')fl?E*02 

0.3AAF-01 

0.395E DO 

350 ."17 

9062.620 

HVHMR : 

C.Sii2E-02 

D.1E2F 00 

0.J77E 00 

1653.317 

3039.095 

FVHMl 

0 .sa?"-Q 2 

-0.e2>lE-02 

0.359E 00 

-Oa.299 

3659.711 

statistics 

FQO THE MATRIX HFthOD 0ASEO ON 150 

OBSERVATIONS 


SCAT CCEF 

value 

KEAN 

FhS 

■/.mean 

ZRKS 

Utf 

0. 152E 02 

-0. 102E-01 

0.2Z2E 00 

-O.IOQ 

1,155 


0.I52E 02 

0. 15EE-0I 

0.222E 00 

0.D61 

1.156 

Urt 

0 . 5*ir-oi 

-0.et5F-02 

Q.22CE DO 

-19.889 

369.720 

VWHHfi 

0.1S2E 02 

-0.E75E-01 

0. E5**E-01 

-0.352 

0 .209 

WVMHl 

O.e72£-Ot 

•0.221E-02 

0 4 AA9E*01 

-2,533 

51.506 

WUHR 

O.SB?E»02 

0.171E 00 

0.365E 00 

1793.109 

..709.0/0 

VVVHl 

Q.5fl2e-02 

0,3A5E-01 

0.395c 00 

351.193 

9066.099 

FUHhftf 

0.5S2E-02 

0.163E 00 

D.370E 00 

1656.123 

3053.521 

HVMMI 

0.5ezE-02 

“0.e3tE-02 

0.359E GO 

-09.675 

3657.770 


Figure D.19e “ SAMPLE OUTPUT FOR 
SCATTEROMETER SIMULATION PROGRAM 


flMTtWfi HITK eiaSE* CM.V 


APPL Bias= oa Phase eias^ c= oeg 


'“lEAS/CCEF VV 

HH 

P£H<t'R HATRIM 
VH VUHHR 

VVHH 

VVVMlJ 

UVVMl 

HUhHR 

liVnhE 

>'oucn 


1 

o.H^ise 00 

0*7SSlE-05 

0 .‘T670E-0f 

0,t6Die-02 

0. 

O.169OE-0b 

0. 

0.b56SE-D7 

0. 

G.OAllE 

00 

Z 

O.ESSIE-05 

0.&31SE 00 

0.967Bf-(l5 

O.lf D1E*02 

0. 

0.bS&BE-07 

0. 

0. l69Hr-nb 

0. 

0 .93T1E 

05 

3 

o.aoiiE-os 

o.eoiu-03 

0.2513E-D? 

-0.1?eS£-02 

0. 

O.BSllE'OS 

0. 

D.DSlIE-OS 

0. 

D.23P3E-02 

b 

Q.^o‘^^e 00 

g.2QD3E 00 

”0 .436&E- LO 

O.bieSE DO 

O.lOOOE-13 

0.ii2Z3E-l2 

b. 

0.i.223E*12 

3. 

0.3 THE 

00 

5 

B.2<irst 00 

fl.JorsE 00 

O.Hb77r-35 

-0.41S1F 00 

-o..37bee-u 

0. 

-O.btSOE-l? 

0. 

~0.b230e-iZ 

0.l295E-9b 

6 

0.2C7-JC CO 

0.?D7S£ DO 

c .bSJRF-r; 

0.71FAE-03 

O.tHOeE-02 

0. 

0. 

0. 

D. 

D.blAbt 

05 

7 

0.20T1E 00 

0.?07®P 00 

C .413'“.F“0» 

0 . 7 J TnS-Q 3 

-a.iseec -02 

0. 

0. 

0 . 

0 . 

0.4l4fcf 

DO 

D 

S< (<l^tE DO 

0.41SqE“03 

D .levOr-D? 

0 .D JOIE-OZ 

0. 

O.b207E-O3 

0 . 

0.S3-17C-D5 

0. 

0 .42FtE 

00 


3.4tbl£ 09 

0.1K,72E-D1 

D.103t.r-0f 

-0 .03010-02 

-0.7I.S2E-13 

-0.bJD2£-D3 

-0 .flbsze-12 

D.3114E-QS 

“Q.aibse-io 

0 .-lOBCE 

00 

fC 

C.b’^lE 00 

O.4015E-03 

O.lOblE-SO 

0 .019SE-11 

D.3776E-0b 

0.b2’SSE-D6 

D.beS9E-Q3 

0.42*.6£-05 

G .b| D 4E-Cb 

0.41S5E 

03 

IT 

0. '.I « If 03 

0.40 lSE-03 

E.lPhtE'Op 

0.8?7flE-ll 

-0.3776e»Db 

0.426&E-D5 

-0 .b259E-B5 

3.b2Sbt-0S 

-O.bl[J4E-0e 

0.4J73E 

OD 

12 

Q.41s9p-03 

0.41&1E 00 

0 .;2«6f-02 

G .U301E-02 

0 . 

0,53g7£-0S 

0 . 

0.4Za7E-03 

0. 

0 .42b6E 

90 

IJ 

0. lD7?E-aj 

O.lilblE 00 

0.123CE-B2 

-o-a3oir-02 

-0.7ii93e-13 

o.3nbe-os 

-D.aib3E-l5 

-0.47D2E-03 

-E.flb52E-12 

S . 4C90E 

00 

lb 

0.43 isr-oi 

Q.blblF 00 

0,t2blE-Q2 

D.OiOSE-li 

O.377«i£-Ob 

0.t,25BE-06 

O.blflbF-OB 

O.b?*jbE-0S 

0 .b25SE-03 

0.41B3E 

00 

IS 

0.40JSE-0J 

0.41'jIF 00 

O.lJblF-O? 

0. A27HE-.il 

-0.377£E-0b 

D.bZSBE-OS 

*9.bl0bE-0b 

B,4?'jbE-05 

-0»b259E-03 

0,'.173E 

00 


SHTIStlCS Po** T«e DIFFFREKCE method BfiSeO QH t OSSERWATIOhS 


SCAT COEP 

value 

hEAtj 

PH5 

2HEAU 

ZRI'S 

VV 

0.1V2E 0? 

3.317F-02 

lOlE-Ob 

0.017 

0. 

HH 

0,192£ 02 

B.JIOE-OZ 

9lE-0b 

0.017 

0. 

Vh 

0. 5V5E-01 

1. 7SbF'‘D2 

D.rl7E-D4 

13.S30 

0. 

VVHHP 

C.1G7E 02 

0,?3*.£-D2 

0, 633E-06 

D.D12 

0. 

VVhkI 

0. 8720-91 

-O. 214E-03 

0. 466E-07 

-0.2bB 

0. 

vvvhp 

0,17'2£-q2 

0.1 S2F JO 

0.37DE-0i 

1959 .587 

D. 

VVWMI 

B.5flrE-02 

0.66 E-D3 

0.755E-06 

8. 852 

0. 

HiTHMA 

O.OHPE-OE 

a.tSFE OQ 

0.370E-D1 

1959. bOb 

0. 

hvkhi 

0. 5fl2E-92 

Q.fi65E-03 

0 * 7SSE-D6 

8.852 

Q. 


statistics PCS The fiATRI« (-EtHOO BASED Oh t OBSERVATIONS 


SCAT COEP 

VALUE 

mean 

PHi 

XhEAM 

isPhS 

VV 

B.J52E 02 

a.lIIE-02 

0.199E-DE 

0.007 

0. 

HH 

G.lb2E 02 

Q.lJaE-02 

Q.i9DE-05 

0.0 D7 

0. 

VM 

0.5ftlE-01 

0.374E-D2 

0.V40E-C4 

6.422 

0. 

VVHHR 

a.l92E 02 

0.255F-D2 

D . 6&4E-*-G5 

0.013 

0. 

vvhhI 

D.872E-01 

-0.662E-[lb 

B , AI8e-(j8 

-Q.076 

0. 

VVVMP 

D .982E-02 

0.193F 08 

Q.371E-01 

1963.140 

0. 

VVVMl 

0 .982E-02 

D,C69F-03 

a.755E-0b 

8.852 

0. 

HVHHR 

0.982E-02 

D.MIE 00 

D.37t£-01 

1963.16b 

0. 

wvmhi 

O.9fl2E-02 

Q,869f-93 

0.755E-0& 

8.852 

0. 


Figure D. I9f “ SAMPLE OUTPUT FOR 
SCATTEROMETER SIMULATION PROGRAM 


APPENDIX E 
RouHne WHERE 


KO PROGRAM DESCRIPTION 


Fortran program WHERE was developed to compute the sampling points for apertures 

having maximum dimensions and across the x and y axis, respectively. The program 

will compute and list (0 .0 ) for m^O and 0 out to values of m and n restricted by 

‘ mn mn 

0,9 £ COS 0 £ 1.0 (E“l) 

If the vol j 3 of m or n exceeds 48 the value is restricted to 48 to limit the storage and 
printpd output to a reasonable amount. The listing of the program is shown in Figure E*”i . 


2.0 EXAMPLE RUN 


The maximum aperture dimensions (x^, y^) and operating wavelength (X) form the 
program input requirement. These must be dimensionally in the same units. An input 
data card containing these parameters must be prepared in accord with the read statement 
and its accompanying format statement. 

An example output of the program is illustrated in Table E.l for an aperture having 
a maximum dimension of 1 ,1760 meters and illuminated at a .02158 meter wavelength. 


194 


t 

CHMtP£ 


2 

C 


3 

ft 

ROUtSHP HHtaE 

h 

C 



c 

iHifbHft iijtpuifj'; POINTS spscipi?*' pqs 

b 

c 

a RtCTAHOttLiC iPSi,Til-F rp HIMH ^ lOT ar.0 L£*J'‘-rH VS'ST. 


c 

THP TP'Ooy IS =»i6tS Ofi COtS TM 25U-1 J. P* SLlASSEIi 

9 

c 

MftVFLtNGTH tLa^-iCil -NO aPt9TUP£ aiH£N0l3fiS >“JST 62 

9 

c 

SPECIPIfcO IN Sfl“2 i,WTTS, 

to 

c 


tt 


DIhE*<stOH KKi<te». THH*a,.4av, PHKoa.oej 

ts 


REat 

t3 


06T* T«I‘4X, OEO/e.fifiCg, i.CmST29?6/ 

l«i 

c 


IS 

c 

SPECIFIEO aPEPTJte DINEN5IQNS ANO HiVEoENTH 

16 

£1 

«DlnE!IST'?‘JS COPC'ED IN PlGaT'-H 0N9E0 COOROlNATE 

t7 

c 

SrSTEV A-T-7J 

18 

e 


IS 


PESO |5,t0aSJ XNOT.tNeT.LAPITOfi 

20 

lOtO 

fOPMfly 

St 


MCITE fC5*150C» TCMO''. VU'^T. LA'IBDa 

ss 

isoe 

FOKHPT J3>,*’C=»Fl.£.i*,3y,*vn*FiC.t,3>(,*|.Ah3ais»F10,a» 

S3 

c 


St) 

c 

ESTAPLISH vauio Obnalh Sat'PLING 

2S 

c 

NOT TO EXCEED A A«*46 natRIX 

26 

c 


27 


HHAX P 2.0*‘SlNtTa4iX»*Xt!JT/tAM9DS -- 

SO 


IF CF«AX .ST* “H*K=A‘l 

29 


HHAX = a-O^SINITMAAl^rioT/LANaPA 

30 • 


ir INHaX .GT, N8J 

31 

c 


32 

c 

DETtfl«IH£ SAhP Lifts POINTS 

33 

ft 


3<» 


HLAH - LAH30A/2.D 

35 


DO 1C 

36 


9P a 1-1 

3? 


DO to Jalt'WAX 

38 


RH a 3-1 

35 

q 


<*C 

c 

FORH SIN(ThETAI:,JI) 


c 


hZ 


-SINTH = HLA-'t“30PT{(tR‘t/:V*t0TI*®2+£RN/Vt0T}»»2):» . 



IF tSINTH .SEt l.CI SO TO 15 



= ATArltSlN*H/3<:^Ttt,0-STl.TH»ll:tTHU/OEC 

ftS 


3F a ,EQ, 1 .AND. J .20. i» GO TC 5 



PPUltJV = ATAN2(at.*«M0T.FH»YH0T)/0EG 

*i7 


GO TO 10 

00 

'5 

PHItltJJ -.i «,0 

OS 

to 

COftTIHliE ■ 

S3 

c 


51 

ft 

DISPLaV SfiPPLIMS FuIftTS 

52 

ft 


S3 


UO 2P tsl.NMX.n 

^Zi. . 


« = X+T ■ 

55 


if IN oCT, Nrtay) N=Ni<ax 

S6 


00 15 K = X,N 

57 


KXtXJ = K-1 

56 

Ifi 

CONTINUE ■ 

69 


HPITE |6.2DC5> VNOT, YNOT, L4KB0A. t«tKt,X=T*Nl 

65 

2000 

FClOHAT llPl . as*, •SAMPLING HAToiy trHETA, 

fit ■ 


7a8X,’'^OP. RFEPT'JP£*,F11.5,* QY * ,F ll* 6 , • AND MAVELENCTN* , 

62 

■* 

F10.6,/pyi t7I1.6//' 

63 


OO 2d JaliNHAX 

66 


H a J-l - 

66 


M»m tb.TCti-'l H, (TN tJ,K» .PHKJ.X} ,K=r,*;t 

6S 

.3050 

FPOMAT fr5,6r'Afl‘Ht.F4.1,2M, tFA,l. 1^)11 

6? 

20 

CONTINUE 

66 


STOP 

69 ’ 


END 


FIGURE E-1. FORTRAN LISTING OF PROGRAM WHERE, 

: 155 

OBKMAU t AGB 38. 


